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Quíck Start 



If you need help installing the software, see Chapter 2, "Installing MATLAB & 
Simulink Student Versión." 

At the heart of MATLAB is a programming language you must learn before 
you can fully exploit its power. You can learn the basics of MATLAB quickly, 
and mastery comes shortly after. You will be rewarded with high-productivity, 
high-creativity computing power that will change the way you work. 

If you are new to MATLAB, you should start by reading Chapter 4, "Matrices 
and Arrays." The most important things to learn are how to enter matrices, 
how to use the : (colon) operator, and how to invoke functions. After you master 
the basics, you should read the rest of the MATLAB chapters in this book and 
run the demos: 

• Chapter 3 , "Introduction" 

Introduces MATLAB and the MATLAB desktop. 

• Chapter 4, "Matrices and Arrays" 

Introduces matrices and arrays, how to enter and genérate them, how to 
opérate on them, and how to control Command Window input and output. 

• Chapter 5, "Graphics" 

Introduces MATLAB graphic capabilities and the tools that let you 
customize your graphs to suit your specific needs. 

• Chapter 6, "Programming" 

Describes how to use the MATLAB language to créate scripts and functions, 
and manipúlate data structures, such as cell arrays and multidimensional 
arrays. 
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About the Student Versión 



MATLAB and Simulink are the premier software packages for technical 
computing in education and industry. MATLAB & Simulink Student Versión 
provides all of the features of professional MATLAB, with no limitations, and 
the full functionality of professional Simulink, with model sizes up to 1000 
blocks. The Student Versión gives you immediate access to high-performance 
numeric computing, modeling, and simulation power. 

MATLAB allows you to focus on your course work and applications rather than 
on programming details. It enables you to solve many numerical problems in a 
fraction of the time it would take you to write a program in a lower-level 
language such as C, C++, or Fortran. MATLAB helps you better understand 
and apply concepts in applications ranging from engineering and mathematics 
to chemistry, biology, and economics. 

Simulink is an interactive tool for modeling, simulating, and analyzing 
dynamic systems, including controls, signal processing, Communications, and 
other complex systems. 

The Symbolic Math Toolbox, also included with the Student Versión, is based 
on the Maple 8 symbolic math engine and lets you perform symbolic 
computations and variable-precisión arithmetic. 

MATLAB producís are used in a broad range of industries, including 
automotive, aerospace, electronics, environmental, telecommunications, 
computer peripherals, finance, and medicine. More than one million technical 
professionals at the world's most innovative technology companies, 
government research labs, financial institutions, and at more than 3500 
universities, rely on MATLAB and Simulink as the fundamental tools for their 
engineering and scientific work. 

Student Use Polícy 

This MATLAB & Simulink Student Versión License is for use in conjunction 
with courses offered at degree-granting institutions. The MathWorks offers 
this license as a special service to the student community and asks your help 
in seeing that its terms are not abused. 

To use this Student License, you must be a student either enrolled in a 
degree-granting institution or participating in a continuing education program 
at a degree-granting educational university. 
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You may not use this Student License at a company or government lab. Also, 
you may not use it if you are an instructor at a university, or for research, 
commercial, or industrial purposes. In these cases, you can acquire the 
appropriate professional or academic license by contacting The MathWorks. 

Student Versión Actívatíon 

Activation is a secure process that verifies licensed student users. This process 
validates the serial number and ensures that it is not used on more systems 
than allowed by the MathWorks End User License Agreement. 

The activation technology is designed to provide an easier and more effective 
way for students to authenticate and use their product than prior releases of 
MATLAB & Simulink Student Versión. 

The quickest way to actívate your software is to use the activation program 
that starts following product installation. The activation program will guide 
you through the activation process. Alternatively, you can actívate your 
software on the mathworks.com Web site. 

Activation requires completion of three activities: 

• Provide registration information by creating a MathWorks Account. 

• Provide your serial number and the Machine ID for the computer you will be 
installing the software on. 

• For those students who did not provide proof of student status at the time of 
purchase, submission and verification of proof of student status is needed. 

For more information on activation, see 

www.mathworks .com/academia/student versión /activation. html. 
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Obtaining Additional MathWorks Products 

Many college courses recommend MATLAB and Simulink as standard 
instructional software. In some cases, the courses may require particular 
toolboxes, blocksets, or other producís. Toolboxes and blocksets are add-on 
producís that extend MATLAB and Simulink with domain-specific 
capabilities. Many of these producís are available for MATLAB & Simulink 
Student Versión. You may purchase and download these additional producís at 
special síudení prices from íhe MaíhWorks Store ai 
www.mathworks.com/store. 

Some of íhe producís you can purchase include 

• Bioinformaíics Toolbox 

• Communications Blockseí 

• Conírol System Toolbox 

• Fixed-Poiní Toolbox 

• Fuzzy Logic Toolbox 

• Image Processing Toolbox 

• Neural Neíwork Toolbox 

• Opíimizaíion Toolbox 

• Signal Processing Toolbox 

• Síaíisíics Toolbox 

• Síaíeflow® (A demo versión of Síaíeflow is included wiíh your MATLAB & 
Simulink Síudení Versión.) 

For an up-ío-daíe lisí of available producís and íheir producí dependencies, 
visií íhe MaíhWorks Store. 



Note The toolboxes and blockseís íhaí are available for MATLAB & 
Simulink Síudení Versión have íhe same funcíionaliíy as íhe professional 
versions. However, íhe síudení versions of íhe íoolboxes and blockseís will 
work only wiíh íhe Síudení Versión. Likewise, íhe professional versions of íhe 
íoolboxes and blockseís will noí work wiíh íhe Síudení Versión. 
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Getting Started wíth MATLAB 



What I Want 



What I Should Do 



I need to install and actívate 
MATLAB. 

I want to start MATLAB. 



I'm new to MATLAB and 
want to learn it quickly. 



I want to look at some 
samples of what you can do 
with MATLAB. 



See Chapter 2, "Installing MATLAB & Simulink Student 
Versión." 



(Microsoft Windows) Double-click the MATLAB icón on your 
desktop. 

(Macintosh OS X) Double-click the MATLAB icón on your desktop. 

(Linux) Enter the matlab command at the command prompt. 

Start by reading Chapter 3, "Introduction," through Chapter 6, 
"Programming," in this book. The most important things to learn 
are how to enter matrices, how to use the : (colon) operator, and 
how to invoke functions. You will also get a brief overview of 
graphics and programming in MATLAB. After you master the 
basics, you can access the rest of the documentation through the 
online help facility (Help). 

There are numerous demonstrations and video tutorials 
included with MATLAB. You can see these by clicking Demos in 
the Help browser or selecting Demos from the Help menú. 
There are demonstrations of mathematics, graphics, 
programming, and much more. You also will find a large 
selection of demos at www.mathworks .com/demos. 
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Fíndíng Reference Information 



What I Want 



What I Should Do 



I want to know how to use a 
specific function. 



I want to find a function for 
a specific purpose, but I 
don't know the function 
ñame. 



I want to learn about a 
specific topic such as sparse 
matrices, ordinary 
differential equations, or cell 
arrays. 

I want to know what 
functions are available in a 
general área. 

I want to learn about the 
Symbolic Math Toolbox. 



Use the online help facility (Help). To access Help, use the 
command helpbrowser or use the Help menú. "MATLAB 
Functions: Volume 1 (A-E), Volume 2 (F-O), and Volume 3 (P-Z)" 
are also available in PDF format from Printing the 
Documentation Set on the MATLAB product page. 



There are several choices: 

• From Help, browse the MATLAB functions by selecting 
Functions — Categorical List or Functions — Alphabetical 
List. 

• Use lookf or (e.g., lookf or inverse) from the command line. 

• Use Index or Search from Help. 

Use Help to lócate the appropriate sections in the MATLAB 
documentation, for example, MATLAB -> Mathematics -> 
Sparse Matrices. 



Use Help to view Functions — Categorical List under 
MATLAB. Help provides access to the reference pages for the 
hundreds of functions included with MATLAB. 

See Chapter 9, "Introducing the Symbolic Math Toolbox," and 
Chapter 10, "Using the Symbolic Math Toolbox," in this book. 
For complete descriptions of the Symbolic Math Toolbox 
functions, use Help and select Functions — Categorical List or 
Functions — Alphabetical List from the Symbolic Math 
Toolbox documentation. 
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Troubleshootíng 



What I Want What I Should Do 



I have a specific MATLAB From Help, select Support and Web Services and then choose 

problem I want help with. Technical Support. 

I want to report a bug or Use Help or send e-mail to bugsOmathworks . com or 

make a suggestion. suggest@mathworks . com. 
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Other Resources 



Documentation 

When you install MATLAB & Simulink Student Versión on your computer, you 
automatically install the complete online documentation for these producís. 
Access this documentation set from Help. 



Note References to UNIX in the documentation include both Linux and 
Mac OS X. 



Web-Based Documentation 

Documentation for all MathWorks producís is online and available from the 
Support área of the MathWorks Web site. In addition to tutorials and function 
reference pages, you can find PDF versions of all the manuals. 

MathWorks Web Site 

At www . mathworks . com, you'll find information about MathWorks producís and 
how they are used in education and industry, product demos, and MATLAB 
and Simulink based books. 

MathWorks Academia Web Site 

At www. mathworks . com/academia, you'll find resources for students and 
instructors for courses in engineering, mathematics, and science. 

MATLAB and Simulink Based Books 

At www.mathworks.com/support/books, you'll find books in many disciplines 
that use MATLAB and Simulink for examples and assignments. 

MathWorks Store 

At www. mathworks .com/store, you can purchase add-on producís and 
documentation. 
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MATLAB Central — File Exchange/Newsgroup 
Access 

At www.mathworks.com/matlabcentral, you can access the MATLAB Usenet 
newsgroup (comp . sof t - sys . matlab) as well as an extensive library of 
user-contributed files called the MATLAB Central File Exchange. MATLAB 
Central is also home to the Link Exchange where you can share your favorite 
links to various educational, personal, and commercial MATLAB Web sites. 

The comp. sof t- sys .matlab newsgroup is for professionals and students who 
use MATLAB and have questions or comments about it and its associated 
software. This is an important resource for posing questions and answering 
queries from other MATLAB users. MathWorks staff also participates actively 
in this newsgroup. 

Technícal Support 

At www.mathworks.com/support, you can get technical support. 

Telephone and e-mail access to our technical support staff is not available for 
students running MATLAB & Simulink Student Versión unless you are 
experiencing difficulty installing or downloading MATLAB or related producís. 
There are numerous other vehicles of technical support that you can use. The 
"Resources and Support" section in the CD holder identifies ways to obtain 
additional help. 

After checking the available MathWorks sources for help, if you still cannot 
resolve your problem, please contact your instructor. Your instructor should be 
able to help you. If your instructor needs assistance doing so, telephone and 
e-mail technical support is available to registered instructors who have 
adopted MATLAB & Simulink Student Versión in their courses. 
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Dífferences Between the Student and Professional Versions 
MATLAB 

The Student Versión includes all of the features of the professional versión of 
MATLAB. 

MATLAB Differences 

There are a few small differences between the Student Versión and the 
professional versión of MATLAB: 

• The MATLAB prompt in the Student Versión is 

EDU» 

• The window title bars include the words 
<Student Version> 

• Printouts contain the footer 
Student Versión of MATLAB 

• The Check for Updates menú option on the Help menú is not available in 
the Student Versión. 

Simulink 

The Student Versión contains the complete Simulink product, which is used 
with MATLAB to model, simúlate, and analyze dynamic systems. 

Simulink Differences 

• Models are limited to 1000 blocks. 

Note You may encounter some demos that use more than 1000 blocks. In 
these cases, a dialog will display stating that the block limit has been 
exceeded and the demo will not run. 



• The window title bars include the words 
<Student Version> 
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Printouts contain the footer 
Student Versión of MATLAB 



Note The Using Simulink documentation, which is accessible from the Help 
browser, contains all of the information in the Learning Simulink book plus 
additional advanced information. 



Symbolic Math Toolbox 

The Symbolic Math Toolbox included with this Student Versión lets you access 
all of the functions in the professional versión of the Symbolic Math Toolbox 
except maple, mapleinit, and mhelp. For more information about the Symbolic 
Math Toolbox, see its documentation. 
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This chapter describes how to install and actívate MATLAB & Simulink Student Versión. 

Installing on Windows (p. 2-2) Steps to install and activate the Student Versión on 

Microsoft Windows machines 

Installing on Mac OS X (p. 2-7) Steps to install and activate the Student Versión on 

Macintosh OS X machines 

Installing on Linux (p. 2-16) Steps to install and activate the Student Versión on 

Linux machines 



Á Installinq MATLAB & Simulink Student Versión 



Installing on Windows 



This section describes the system requirements necessary to install MATLAB 
& Simulink Student Versión on a Windows computer. It also provides 
step-by-step instructions for installing and activating the software. 



System Requirements 



Note For more information on system requirements, visit 
www.mathworks . com/ academia /student_version/ requirements. html. 



MATLAB and Simulink 

• Windows XP® (Service Pack 1 or 2) or Windows 2000® (Service Pack 3 or 4) 

• Intel Pentium III processor or later, Celeron processor, or Intel Xeon 
processor family; AMD Athlon / Duron / Opteron ; or compatible 
processor 

• 256 MB RAM (512 MB RAM recommended) 

• 500 MB disk space 

• 16-, 24-, or 32-bit OpenGL® capable graphics adaptor 

• CD-ROM or DVD drive (for installation) 

• E-mail (required), internet access (recommended) for product activation 

MEX-Files 

MEX-files are dynamically linked subroutines that MATLAB can 
automatically load and execute. They provide a mechanism by which you can 
cali your own C and Fortran subroutines from MATLAB as if they were built-in 
functions. 



For More Information "External Interfaces" in the MATLAB 
documentation provides information on how to write MEX-files. "External 
Interfaces Reference" in the MATLAB documentation describes the collection 
of these functions. Both of these are available from Help. 
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If you plan to build your own MEX-files, you will need a supported compiler. 
For the most up-to-date information about compilers, see the support área at 
the Math Works Web site (www.mathworks . com) . 

Installing and Activating the Student Versión 

This list summarizes the steps in the standard installation procedure. You can 
perform the installation and activation by simply following the instructions in 
the dialog boxes presented by the installation program; it walks you through 
this process: 

1 Exit any existing copies of MATLAB you have running. 

2 Insert the MATLAB & Simulink Student Versión CD in your CD-ROM 
drive. To start the installation program, run setup . exe from the CD. 

3 Read the Welcome screen, then click Next. 

4 Enter your ñame and school ñame, then click Next. 

5 Review the software licensing agreement and, if you agree with the terms, 
select Yes and click Next. 

6 Review the Student Use Policy and, if you satisfy the terms, select Yes and 
click Next. 

7 Choose your installation type. Typical installation installs all producís; 
custom installation gives you control over the installation such as selecting 
which producís to install. Select the installation type and click Next. 

8 The Folder Selection dialog box lets you specify the ñame of the folder into 
which you want to install MATLAB. You can accept the default destination 
folder or specify the ñame of a different installation folder. If the folder 
doesn't exist, the installer creates it. To continué with the installation, click 

Next. 

9 The Confirmation dialog box lets you confirm your installation options. To 
change a setting, click the Back button. To proceed with the installation, 
click Install. 
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Note The installation process installs the online documentation for each 
product you install. This does not include documentation in PDF format, 
which is available only at the Math Works Web site. 



10 When the installation successfully completes, the activation process begins 
by displaying the Activation Overview, which describes the three steps in 
the process. During activation, you will: 

- Enter your serial number and e-mail address. 

- Provide registration information by creating a MathWorks account. 

- Provide proof of student status, unless previously provided. 

11 As the activation process proceeds, read the screens and enter the 
corresponding information. At the completion of the activation process, you 
will be able to use your Student Versión of MATLAB & Simulink. In certain 
cases, your software will be temporarily activated for 30 days until your 
proof of student status is verified. In these cases, you will be reminded that 
your activation is temporary and that you need to complete the activation 
process. Once your proof of student status is verified, your activation is 
complete. 



Note If you encounter a problem during the activation process, check 
www.mathworks . com / academia /student_version/ activation. html for more 
information. 



12 To start MATLAB, double-click the MATLAB icón that the installer creates 
on your desktop. 

13 Customize any MATLAB environment options, if desired. For example, to 
specify welcome messages, default definitions, or any MATLAB expressions 
that you want executed every time MATLAB is invoked, créate a file named 
startup.m in the $MATLAB\toolbox\local folder, where $MATLAB is the 
ñame of your MATLAB installation folder. Every time you start MATLAB, 
it executes the commands in the startup . m file. 
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14 Perform any additional configuration by typing the appropriate command at 
the MATLAB command prompt. For example, to configure the MATLAB 
Notebook, type notebook - setup. To configure a compiler to work with the 
MATLAB External Interface, type mex -setup. 



For More Information The Installation Guide for Windows documentation 
provides additional installation information. This manual is available from 
Help. 



Installing Additional Toolboxes 

To purchase additional toolboxes, visit the MathWorks Store at 
www.mathworks.com/store. Once you purchase a toolbox, the product and its 
online documentation are downloaded to your computer. 

When you download a toolbox, you receive an installation program for the 
toolbox. To install the toolbox and its documentation, run the installation 
program by double-clicking the installer icón. After you successfully install the 
toolbox, all of its functionality and documentation will be available to you when 
you start MATLAB. 
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Accessing the Online Documentaron (Help) 

To access the online documentation (Help), select Full Product Family Help 
from the Help menú in the MATLAB Command Window. You can also type 
helpbrowser at the MATLAB prompt. The Help browser appears. 



i jiíw 



File Edit View Go Favo rites Desktop 
Help Navigator X_ 

Conlents índex Searcri Demos ' 

B~0 Rolease Notes 
B-0 Irstallaticín 

+ -> MATLAB 

y Symholir rviaíh Toolbúx 
É--# Simulink 
£-#Stateflow 

@ Support and Web Services 



Window Help 

+■ -► c I m l #4 



Tfflfe Reléase 1 4 wrtl i Service Pack 3: Begin Hele - Student Versión 



~B 
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MATLAB & Simulink Student Versión Reléase 14 with Service Pack 3 



IfYou Are Using MATLAB for the First Time... 

At the heart of MATLAB ¡s a new lana.uage íhaí you musí learn í;efcre you can fuüy eíploit iis power. 
Thls ¡snt as hsrd as it might sound; you can learn the basies of MATLAB ver y quickly. Yon will be 
rewarderi with h¡ g h- pro d u ct ivity . high-creativlty compuíing power that wiil charas tris w/ay you work. 

The best way to get started is to read 1 1 ■ r i I _ - ■clions (Chspters 3 through 8) in the 

Learning MATLAB tutorial with MATLAB open so you car follow along The Learning MATLAB book is 
¡ncluded with üie ¡v'ATlAB Síudsnt Versión. Yon can siso access iris . from the 

Help browser, 

If you dont want to take (he time lo rsad it thoruughly, here are links to the most important sectians: 

■ What Is MATLAB? 

An overvie ■ otti ■ i e. ofMATLAB 

■ Matrices and Arrays 

How to enter matrices, use the colon (:) operater, and invoke functions 

■ Graphics 

How to plot data, annotate graphs, arel work with ¡mages 

■ Programminq 

How io use flow control siatemems and créate scripts and íunciions 

You should also watch the 7-minute to learn howto 

opérate the MATLAB desktop iníerface. Thie le not a substituto for reading the above sectione about 
the MATLAB langiiaga, but wü! orlent yon to The cievelopment environment. There are rnany other 
you might iinJ helpfu; as wal!. 
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Installing on Mac OS X 



This section describes the system requirements necessary to install MATLAB 
& Simulink Student Versión on a Macintosh computer. It also provides 
step-by-step instructions for installing and activating the software. 



System Requirements 



Note For more information on system requirements, visit 
www.mathworks.com/academia/student_version/requirements.html. 



MATLAB and Simulink 

• Mac OS X 10.3.8, 10.3.9 (Panther™), or 10.4 (Tiger™) 

• PowerPC G4 or G5 processor 

• 256 MB RAM (512 MB RAM recommended) 

• 500 MB disk space 

• 16-bit graphics or higher adaptor and display (24-bit recommended) 

• XI 1 (X server) 

• CD-ROM or DVD drive (for installation) 

• E-mail (required), internet access (recommended) for product activation 

MEX-Files 

MEX-files are dynamically linked subroutines that MATLAB can 
automatically load and execute. They provide a mechanism by which you can 
cali your own C and Fortran subroutines from MATLAB as if they were built-in 
functions. 



For More Information "External Interfaces" in the MATLAB 
documentation provides information on how to write MEX-files. "External 
Interfaces Reference" in the MATLAB documentation describes the collection 
of these functions. Both of these are available from Help. 
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If you plan to build your own MEX-files, you will need a supported compiler. 
For the most up-to-date information about compilers, see the support área at 
the MathWorks Web site (www.mathworks.com). 



Installing and Actívatíng the Student Versión 

The following sections describe the steps you must follow to install and actívate 
MATLAB & Simulink Student Versión on a Macintosh computer. 



Note If you want to install MATLAB in a particular directory, you must have 
the appropriate permissions. For example, to install MATLAB in the 
Applications directory, you must have administrator status. To créate 
symbolic links in a particular directory, you must have the appropriate 
permissions. For information on setting permissions (privileges), see 
Macintosh Help (Command+? from the desktop). 



1 Insert the MATLAB & Simulink Student Versión CD in the CD-ROM drive. 
When the CD's icón appears on the desktop, double-click the icón to display 
the CD's contents. 

2 Double-click the InstallForMacOSX icón to begin the installation. 
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3 To install the MathWorks producís, you must enter a username and 
password for an administrator. Enter the username and password in the 
Authenticate dialog box and then click OK. 



Authenticate 



MathWorks ¡nstaller requires that you type your 

i." <\ 



password. 



Ñame: logín ñame 



Password: , 



(▼) Details 

© 



( Cancel f OK ) 



Note The user account you set up when you configured your Mac is an 
administrator account. If you don't know the username and password for an 
administrator, check with the person who administers your Macintosh. 
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4 When the Important Information window appears, review its contents. 
Click Continué to proceed with the installation. 



^O^ 



MathWorks Install 



Important Information 



Note: MATLAB requires Xll for Mac 05 X, which is uvdildble as an optional choice 
on the Mac 05 X installation CD. We recommend thot you install Xll for Mac OS X 
before installing MATLAB. 

It may be necessary to contact the MathWorks web site to activóte this 
insta l lation. 

MathWorks products are protected by U.S. patents (see www.mothworks.com/patents) 
and copyright laws. By entering into the Software License Agreement that 
follows within this installation process you will also agree to additional 
restrictions on your use of these programs. Any unauthorized use, reproduction, 
or distribution muy result in serious civil ond criminol pendlties. 



( Print... ) ( Save... ) 



( Quit ) ( Continué ) 



Note You will not be able to run MATLAB & Simulink Student Versión until 
you install Xll for Mac OS X. 



5 The Software License Agreement is displayed. If you agree to its terms, you 
can continué the installation. 

6 The Student Use Policy is displayed. If you qualify for the student license 
and agree to its terms, click Yes. 
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7 The default installation location is MATLAB_SV71 in the Applications folder 
on your system disk. To accept the default, click Continué. To change the 
location, click Choose Folder and then navigate to the desired location. 



Note Your installation folder ñame cannot contain spaces, the @ character, 
the % character, or the $ character. You cannot install into a folder named 
private, but you can have a folder named private on the path. To créate the 
installation folder in this location on your system, you must have 
administrator privileges. For information on setting privileges, see Macintosh 
Help (Command+? from the desktop). 



^O^ 



MathWorks Install 



Install Location 




Instüllatíon will Dccur ¡n "MATLAB_SV71" in the folder "Applications" 
on the disk "Mac Hard Drive", 



To choose a different folder, click Choose Folder. 



( Choose Folder... ) 



( Co Back ) { Continué ") 
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8 Select the producís you want to install from the list of producís available for 
installation and then click Continué. 



^O^ 



MathWorks Install 



Select producís for installation 



g MATLAB 
B Simulink 
S Symbollc Math Toolbox 



( DeselectAII ) (' Select All ) 



( Co Back j ( Continué } 



Note The installation process installs the online documentation for each 
product you install. This does not include documentation in PDF format, 
which is only available at the MathWorks Web site. 
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9 When the installation successfully completes, the installer displays the 
following Installation is complete window. This gives you options for 
accessingMATLAB. Checkyourdesired options foraccessingMATLAB, and 
click Continué to begin the activation process. 






oeo 



MathWorks Install 



Installation is complete 

Do you want to: 

@ Créate alias of MATLAB on the desktop 

Add MATLAB icón to the Dock 



The MATLAB installer is finished. This installation needs to be 
activated prior to use. Click Continué to begin the activation 
process. 



( Continué ) 



10 The activation process displays the Activation Overview, which describes 
the three steps in the process. During activation, you will: 

- Enter your serial number and e-mail address. 

- Provide registration information by creating a MathWorks account. 

- Provide proof of student status, unless previously provided. 

11 As the activation process proceeds, read the screens and enter the 
corresponding information. At the completion of the activation process, you 
will be able to use your Student Versión of MATLAB & Simulink. In certain 
cases, your software will be temporarily activated for 30 days until your 
proof of student status is verified. In these cases, you will be reminded that 
your activation is temporary and that you need to complete the activation 
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process. Once your proof of student status is verified, your activation is 
complete. 



Note If you encounter a problem during the activation process, check 
www.mathworks . com / academia /student_version/ activation. html for more 
information. 



12 To start MATLAB, double-click the MATLAB icón that the installer creates 
on your desktop. 



For More Information The Installation Guide for Mac OS X documentation 
provides additional installation information. This manual is available from 
Help. 



Installing Additional Toolboxes 

To purchase additional toolboxes, visit the MathWorks Store at 
www.mathworks .com/store. Once you purchase a toolbox, the product and its 
online documentation are downloaded to your computer. 

When you download a toolbox, you receive a file for the toolbox. Double-clicking 
the downloaded file's icón creates a folder that contains the installation 
program for the toolbox. To install the toolbox and its documentation, run the 
installation program by double-clicking its icón. After you successfully install 
the toolbox, all of its functionality and documentation will be available to you 
when you start MATLAB. 

Accessing the Online Documentation (Help) 

To access the online documentation (Help), select Full Product Family Help 
from the Help menú in the MATLAB Command Window. You can also type 
helpbrowser at the MATLAB prompt. The Help browser appears. 
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Mac OS X Documentatíon 

In general, the documentation for MathWorks producís is not specific for 
individual platforms unless the product is available only on a particular 
platform. For the Macintosh, when you access a product's documentation either 
in print or online through the Help browser, make sure you refer to the UNIX 
platform if there is different documentation for different platforms. 
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Installing on Linux 



This section describes the system requirements necessary to install MATLAB 
& Simulink Student Versión on a Linux computer. It also provides step-by-step 
instructions for installing and activating the software. 



System Requirements 



Note For more information on system requirements, visit 
www.mathworks . com/ academia /student_version/ requirements. html. 



MATLAB and Simulink 

• Linux Kernel 2.4.x; glibc (glibc6) 2.3.2 and 2.2.5 or Linux 2.6.x; (glibc6) 2.3.2 

• Intel Pentium III processor or later, Celeron processor or Intel Xeon 
processor family; AMD Athlon / Duron / Opteron ; or compatible 
processor 

• 256 MB RAM (512 MB recommended) 

• 500 MB disk space 

• 16-bit graphics or higher adaptor and display (24-bit recommended) 

• CD-ROM or DVD drive (for installation) 

• E-mail (required), internet access (recommended) for product activation 

MEX-Files 

MEX-files are dynamically linked subroutines that MATLAB can 
automatically load and execute. They provide a mechanism by which you can 
cali your own C and Fortran subroutines from MATLAB as if they were built-in 
functions. 



For More Information "External Interfaces" in the MATLAB 
documentation provides information on how to write MEX-files. "External 
Interfaces Reference" in the MATLAB documentation describes the collection 
of these functions. Both of these are available from Help. 
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If you plan to build your own MEX-files, you will need a supported compiler. 
For the most up-to-date information about compilers, see the support área at 
the Math Works Web site (www.mathworks . com) . 

Installing and Activating the Student Versión 

The following instructions describe how to install and actívate MATLAB & 
Simulink Student Versión on a Linux computer. 



Note On most systems, you will need root privileges to perform certain steps 
in the installation procedure. 



1 Insert the MATLAB & Simulink Student Versión CD in the CD-ROM drive. 

If your CD-ROM drive is not accessible to your operating system, you will 
need to mount the CD-ROM drive on your system. Créate a directory to be 
the mount point for it. 

mkdir /cdrom 

Mount a CD-ROM drive using the command 

$ mount /cdrom 

If your system requires that you have root privileges to mount a CD-ROM 
drive, this command should work on most systems. 

# mount -t iso9660 /dev/cdrom /cdrom 

To enable nonroot users to mount a CD-ROM drive, include the exec option 
in the entry for CD-ROM drives in your /etc/f stab file, as in the following 
example. 

/dev/cdrom /cdrom iso9660 noauto, ro,user,exec 

Note, however, that this option is often omitted from the /etc/f stab file for 
security reasons. 
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2 Créate an installation directory and move to it, using the cd command. For 
example, if you are going to install into the location 
/usr/local/matlab71_sv, use the commands 

cd /usr/local 
mkdir matlab71_sv 
cd matlab71_sv 

Subsequent instructions in this section refer to this directory as $MATLAB. 

3 Start the MathWorks Installer by running the install script. 
/cdrom/install_unix . sh 

The MathWorks Installer displays the Welcome dialog box. Click OK to 
proceed with the installation. 



rv" 



— 




MATLAB Installation 



Welcome to Uie MATLAB R1 4 Student Versión installer. 

(c) Copyright 1984-2005 by The MathWorks, Inc. 

(508) B47-7000 

The MathWorks proclucts are protectecl by one or more 

U.S. patente. See www.mathworks.com/patents. 



3 3 



OK 



Cancel 



Help 



Note If you need additional help on any step during this installation process, 
click the Help button at the bottom of the dialog box. 



4 The Software License Agreement is displayed. If you agree to its terms, you 
can continué the installation. 
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5 Verify the ñame of the directory in which you want to install MATLAB in the 
MATLAB Root Directory dialog box. You can edit the pathname in this 
dialog box. If the MATLAB root directory is correct, click OK to proceed with 
the installation. 



F 



ir-i-: c vi.- >y 



í'IaTLaE; F'/joL Dir 



MATLAB raot directory location: 



/iisr/lacaMmatlali?! _s\{ 



OK 



Cancel 






Help 



Ll 
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6 Select the producís you want to install in the Installation Options dialog 
box. 



Platforms: 



ítems to install: 



r Linux (x86) MATLAB 

Simulink 

J Mac 



Symbolic Matti Toolbox 



FU 



Remove > 



ítems not to install: 



FT. 



J 



t Add 



- 



Total Installation Size: 535300 kbytes 



OK 



J 



Cancel 



Help 



The producís you are licensed to install are listed in the ítems to install list 
box. If you do not want to install a product, select it in the list and click 
Remove. The installer moves the product to the ítems not to install list. 

To install the complete MATLAB & Simulink Student Versión, keep all the 
producís listed in the ítems to install list. 

After you select the producís you want ío insíall, click OK ío proceed wiíh 
íhe insíallaíion. 
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7 Specify in the Installation Data dialog box the directory in which you want 
to install symbolic links to the matlab and mex scripts. Choose a directory 
that is common to all users' paths, such as /usr/local/bin. You must be 
logged in as root to do this. If you choose not to set up these links, you can 
still run MATLAB; however, you must specify the full path to the MATLAB 
start-up script. Click OK to proceed with the installation. 











- 








_l Créate symbolic links te MATLAB and mex scripts in: 






'fsVÍTñSC'SMñ 


















OK Cancel 


Help 













8 Start the installation by clicking OK in the Begin Installation dialog box. 
During the installation, the installer displays information about the status 
of the installation. 

9 When the installation successfully completes, the activation process begins 
by displaying the Activation Overview, which describes the three steps in 
the process. During activation, you will: 

- Enter your serial number and e-mail address. 

- Provide registration information by creating a MathWorks account. 

- Provide proof of student status, unless previously provided. 

10 As the activation process proceeds, read the screens and enter the 
corresponding information. At the completion of the activation process, you 
will be able to use your Student Versión of MATLAB & Simulink. In certain 
cases, your software will be temporarily activated for 30 days until your 
proof of student status is verified. In these cases, you will be reminded that 
your activation is temporary and that you need to complete the activation 
process. Once your proof of student status is verified, your activation is 
complete. 
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Note If you encounter a problem during the activation process, check 
www.mathworks . com / academia /student_version/ activation. html for more 
information. 



11 To start MATLAB, enter the matlab command. If you did not set up symbolic 
links in a directory on your path, you must provide the full pathname to the 
matlab command 

$MATLAB/bin/matlab 

where $MATLAB represents your MATLAB installation directory. 

Installíng Additional Toolboxes 

To purchase additional toolboxes, visit the MathWorks Store at 
www.mathworks .com/store. Once you purchase a toolbox, the product and its 
online documentation are downloaded to your computer. When you download a 
toolbox on Linux, you receive a tar file (a standard, compressed archive 
format). 

To install the toolbox and documentation, you must 

1 Place the tar file in your installation directory ($MATLAB) and extract the 
files from the archive. Use the following syntax. 

tar -xf filename 

2 Start the MathWorks Installer. 

install 

After you successfully install the toolbox and documentation, all of its 
functionality will be available to you when you start MATLAB. 

Accessíng the Online Documentation (Help) 

To access the online documentation (Help), select Full Product Family Help 

from the Help menú in the MATLAB desktop. You can also type helpbrowser 
at the MATLAB prompt. 
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This chapter introduces MATLAB and Simulink, the documentation set, how to start and stop 
MATLAB, and the MATLAB desktop. 

About MATLAB and Simulink (p. 3-2) Overview of the producís 

MATLAB Documentation (p. 3-5) Overview of the documentation set 

Starting and Quitting MATLAB Steps to run and exit MATLAB 

(p. 3-7) 
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About MATLAB and Símulínk 
What Is MATLAB? 



MATLAB is a high-performance language for technical computing. It 
integrates computation, visualization, and programming in an easy-to-use 
environment where problems and solutions are exprés sed in familiar 
mathematical notation. Typical uses include 

• Math and computation 

• Algorithm development 

• Data acquisition 

• Modeling, simulation, and prototyping 

• Data analysis, exploration, and visualization 

• Scientific and engineering graphics 

• Application development, including graphical user interface building 

MATLAB is an interactive system whose basic data element is an array that 
does not require dimensioning. This allows you to solve many technical 
computing problems, especially those with matrix and vector formulations, in 
a fraction of the time it would take to write a program in a scalar noninteractive 
language such as C or Fortran. 

The ñame MATLAB stands for matrix laboratory. MATLAB was originally 
written to provide easy access to matrix software developed by the LINPACK 
and EISPACK projects. Today, MATLAB engines incorpórate the LAPACK 
and BLAS libraries, embedding the state of the art in software for matrix 
computation. 

MATLAB has evolved over a period of years with input from many users. In 
university environments, it is the standard instructional tool for introductory 
and advanced courses in mathematics, engineering, and science. In industry, 
MATLAB is the tool of choice for high-productivity research, development, and 
analysis. 

Toolboxes 

MATLAB features a family of add-on application-specific solutions called 
toolboxes. Very important to most users of MATLAB, toolboxes allow you to 
learn and apply specialized technology. Toolboxes are comprehensive 
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collections of MATLAB functions (M-files) that extend the MATLAB 
environment to solve particular classes of problems. Áreas in which toolboxes 
are available include signal processing, control systems, neural networks, 
fuzzy logic, wavelets, simulation, and many others. 

The MATLAB System 

The MATLAB system consists of five main parts: 

Development Environment. This is the set of tools and facilities that help you use 
MATLAB functions and files. Many of these tools are graphical user interfaces. 
It includes the MATLAB desktop and Command Window, a command history, 
an editor and debugger, and browsers for viewing help, the workspace, files, 
and the search path. 

The MATLAB Mathematical Function Library. This is a vast collection of computational 
algorithms ranging from elementary functions, like sum, sine, cosine, and 
complex arithmetic, to more sophisticated functions like matrix inverse, matrix 
eigenvalues, Bessel functions, and fast Fourier transforms. 

The MATLAB Language. This is a high-level matrix/array language with control 
flow statements, functions, data structures, input/output, and object-oriented 
programming features. It allows both "programming in the small" to rapidly 
créate quick and dirty throw-away programs, and "programming in the large" 
to créate large and complex application programs. 

Graphics. MATLAB has extensive facilities for displaying vectors and matrices 
as graphs, as well as annotating and printing these graphs. It includes 
high-level functions for two-dimensional and three-dimensional data 
visualization, image processing, animation, and presentation graphics. It also 
includes low-level functions that allow you to fully customize the appearance of 
graphics as well as to build complete graphical user interfaces on your 
MATLAB applications. 

The MATLAB Application Program Interface (API). This is a library that allows you to 
write C and Fortran programs that interact with MATLAB. It includes 
facilities for calling routines from MATLAB (dynamic linking), calling 
MATLAB as a computational engine, and for reading and writing MAT-files. 
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What Is Símulínk? 

Simulink is an interactive environment for modeling, simulating, and 
analyzing dynamic, multidomain systems. It lets you build a block diagram, 
simúlate the system's behavior, evalúate its performance, and refine the 
design. Simulink integrates seamlessly with MATLAB, providing you with 
immediate access to an extensive range of analysis and design tools. These 
benefits make Simulink the tool of choice for control system design, DSP 
design, Communications system design, and other simulation applications. 

Blocksets are collections of application-specific blocks that support múltiple 
design áreas, including electrical power-system modeling, digital signal 
processing, fixed-point algorithm development, and more. These blocks can be 
incorporated directly into your Simulink models. 

Real-Time Workshop® is a program that generates optimized, portable, and 
customizable ANSÍ C code from Simulink models. Generated code can run on 
PC hardware, DSPs, microcontrollers on bare-board environments, and with 
commercial or proprietary real-time operating systems. 

What Is Stateflow? 

Stateflow is an interactive design tool for modeling and simulating complex 
reactive systems. Tightly integrated with Simulink and MATLAB, Stateflow 
provides Simulink users with an elegant solution for designing embedded 
systems by giving them an efficient way to incorpórate complex control and 
supervisory logic within their Simulink models. 

With Stateflow, you can quickly develop graphical models of event-driven 
systems using finite state machine theory, statechart formalisms, and flow 
diagram notation. Together, Stateflow and Simulink serve as an executable 
specification and virtual prototype of your system design. 



Note Your MATLAB & Simulink Student Versión includes a demo versión of 
Stateflow. 
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MATLAB Documentation 



MATLAB provides extensive documentation, in both printed and online 
format, to help you learn about and use all of its features. If you are a new user, 
start with the MATLAB specific sections in this book. It covers all the primary 
MATLAB features at a high level, including many examples. 

The MATLAB online help provides task-oriented and reference information 
about MATLAB features. MATLAB documentation is also available in printed 
form and in PDF format. 

MATLAB Online Help 

To view the online documentation, select MATLAB Help from the Help menú 
in MATLAB. The MATLAB documentation is organized into these main topics: 

• Desktop Tools and Development Environment — Startup and shutdown, the 
desktop, and other tools that help you use MATLAB 

• Mathematics — Mathematical operations and data analysis 

• Programming — The MATLAB language and how to develop MATLAB 
applications 

• Graphics — Tools and techniques for plotting, graph annotation, printing, 
and programming with Handle Graphics® 

• 3-D Visualization — Visualizing surface and volume data, transparency, and 
viewing and lighting techniques 

• Creating Graphical User Interfaces — GUI-building tools and how to write 
callback functions 

• External Interfaces — MEX-files, the MATLAB engine, and interfacing to 
Java, COM, and the serial port 

MATLAB also includes reference documentation for all MATLAB functions: 

• Functions — Categorical List — Lists all MATLAB functions grouped into 
categories 

• Handle Graphics Property Browser — Provides easy access to descriptions of 
graphics object properties 

• External Interfaces Reference — Covers those functions used by the 
MATLAB external interfaces, providing information on syntax in the calling 
language, description, arguments, return valúes, and examples 
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The MATLAB online documentation also includes 

• Examples — An index of examples included in the documentation 

• Reléase Notes — New features and known problems in the current reléase 

• Printable Documentation — PDF versions of the documentation suitable for 
printing 

For more information about using the Help browser, see Chapter 8, "Desktop 
Tools and Development Environment." 



Note References to UNIX in the documentation include both Linux and 
Mac OS X. 
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Startíng and Quitting MATLAB 



Starting MATLAB 

On Windows platforms, start MATLAB by double-clicking the MATLAB 
shortcut icón & on your Windows desktop. 

On Mac OS X platforms, start MATLAB by double-clicking the MATLAB icón 
on your desktop. 

On Linux platforms, start MATLAB by typing matlab at the operating system 
prompt. 

You can customize MATLAB startup. For example, you can change the 
directory in which MATLAB starts or automatically execute MATLAB 
statements in a script file named startup. m. 



For More Information See "Starting MATLAB" in the Desktop Tools and 
Development Environment documentation. 



Quitting MATLAB 

To end your MATLAB session, select File -> Exit MATLAB in the desktop, or 
type quit in the Command Window. You can run a script file named f inish . m 
each time MATLAB quits that, for example, executes functions to save the 
workspace, or displays a quit confirmation dialog box. 



For More Information See "Quitting MATLAB" in the Desktop Tools and 
Development Environment documentation. 



MATLAB Desktop 

When you start MATLAB, the MATLAB desktop appears, containing tools 
(graphical user interfaces) for managing files, variables, and applications 
associated with MATLAB. 
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The following illustration shows the default desktop. You can customize the 
arrangement of tools and documents to suit your needs. For more information 
about the desktop tools, see Chapter 8, "Desktop Tools and Development 
Environment." 



Enter MATLAB functions at the Command Window prompt. 



# MATLAB 



File Edit Debug Desktop Window Help 



Deí * * B « o- % tf 



3L2£J 



? X 



a tf s\ s a - 



All Files 



File Type 



La 



CDseaJernp 
¡^bucky.m 
Jcaution.mdl 
±J 



Folder Ap_±] 

M-file No— 1 

Model No ■ 

i ► I 



d:\myrnfiles 



Command Window 



"U ,1 a 



Current Directory I Workspace | 



? X 

lfl-% — 4/16/04 10:00 AM — % ^J 
-ele 

-forraat long e 
-cd d: /mymf iles/sea te 



To 
>> 



get started, select MATLAB Help or Demos 



The Command Hislory maintains a record of 
the MATLAB functions you ran. 
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This chapter introduces you to MATLAB by teaching you how to handle matrices 
Matrices and Magic Squares (p. 4-2) 



Expressions (p. 4-10) 

Working with Matrices (p. 4-14) 



More About Matrices and Arrays 
(p. 4-18) 

Controlling Command Window Input 
and Output (p. 4-28) 



Enter matrices, perform matrix operations, and access 
matrix elements. 

Work with variables, numbers, operators, functions, and 
expressions. 

Genérate matrices, load matrices, créate matrices from 
M-files and concatenation, and delete matrix rows and 
columns. 

Use matrices for linear algebra, work with arrays, 
multivariate data, scalar expansión, and logical 
subscripting, and use the f ind function. 

Change output format, suppress output, enter long lines, 
and edit at the command line. 



1 Matrices and Arrays 



Matrices and Magic Squares 



In MATLAB, a matrix is a rectangular array of numbers. Special meaning is 
sometimes attached to 1-by-l matrices, which are scalars, and to matrices with 
only one row or column, which are vectors. MATLAB has other ways of storing 
both numeric and nonnumeric data, but in the beginning, it is usually best to 
think of everything as a matrix. The operations in MATLAB are designed to be 
as natural as possible. Where other programming languages work with 
numbers one at a time, MATLAB allows you to work with entire matrices 
quickly and easily. A good example matrix, used throughout this book, appears 
in the Renaissance engraving Melencolia I by the Germán artist and amateur 
mathematician Albrecht Dürer. 
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This image is filled with mathematical symbolism, and if you look carefully, 
you will see a matrix in the upper right córner. This matrix is known as a magic 
square and was believed by many in Dürer's time to have genuinely magical 
properties. It does turn out to have some fascinating characteristics worth 
exploring. 




Enteríng Matrices 

The best way for you to get started with MATLAB is to learn how to handle 
matrices. Start MATLAB and follow along with each example. 

You can enter matrices into MATLAB in several different ways: 

• Enter an explicit list of elements. 

• Load matrices from external data files. 

• Genérate matrices using built-in functions. 

• Créate matrices with your own functions in M-files. 

Start by entering Dürer's matrix as a list of its elements. You only have to 
follow a few basic conventions: 

• Sepárate the elements of a row with blanks or commas. 

• Use a semicolon, ; , to indicate the end of each row. 

• Surround the entire list of elements with square brackets, [ ] . 
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To enter Dürer's matrix, simply type in the Command Window 
A = [16 3 2 13; 5 10 11 8; 9 6 7 12; 4 15 14 1] 

MATLAB displays the matrix you just entered. 
A = 



16 


3 


2 


13 


5 


10 


11 


8 


9 


6 


7 


12 


4 


15 


14 


1 



This matrix matches the numbers in the engraving. Once you have entered the 
matrix, it is automatically remembered in the MATLAB workspace. You can 
refer to it simply as A. Now that you have A in the workspace, take a look at 
what makes it so interesting. Why is it magic? 

sum, transpose, and díag 

You are probably already aware that the special properties of a magic square 
have to do with the various ways of summing its elements. If you take the sum 
along any row or column, or along either of the two main diagonals, you will 
always get the same number. Let us verify that using MATLAB. The first 
statement to try is 

sum(A) 

MATLAB replies with 

ans = 

34 34 34 34 

When you do not specify an output variable, MATLAB uses the variable ans, 
short for answer, to store the results of a calculation. You have computed a row 
vector containing the sums of the columns of A. Sure enough, each of the 
columns has the same sum, the magic sum, 34. 

How about the row sums? MATLAB has a preference for working with the 
columns of a matrix, so the easiest way to get the row sums is to transpose the 
matrix, compute the column sums of the transpose, and then transpose the 
result. 

MATLAB has two transpose operators. The apostrophe operator (e.g., A' ) 
performs a complex conjúgate transposition. It flips a matrix about its main 
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diagonal, and also changes the sign of the imaginary component of any complex 
elements of the matrix. The apostrophe-dot operator (e.g., A' .), transposes 
without affecting the sign of complex elements. For matrices containing all real 
elements, the two operators return the same result. 



So 



A 1 








produces 








ans = 








16 


5 


9 


4 


3 


10 


6 


15 


2 


11 


7 


14 


13 


8 


12 


1 



and 

sum(A' ) ' 

produces a column vector containing the row sums 

ans = 
34 
34 
34 
34 

The sum of the elements on the main diagonal is obtained with the sum and the 
diag functions. 

diag(A) 

produces 

ans = 
16 

10 
7 
1 

and 

sum(diag(A) ) 

produces 
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ans = 
34 

The other diagonal, the so-called antidiagonal, is not so important 
mathematically, so MATLAB does not have a ready-made function for it. But a 
function originally intended for use in graphics, f liplr, flips a matrix from left 
to right. 

sum(diag(fliplr(A))) 

ans = 
34 

You have verified that the matrix in Dürer's engraving is indeed a magic 
square and, in the process, have sampled a few MATLAB matrix operations. 
The following sections continué to use this matrix to illustrate additional 
MATLAB capabilities. 

Subscripts 

The element in row i and column j of A is denoted by A ( i , j ) . For example, 
A(4 , 2 ) is the number in the fourth row and second column. For our magic 
square, A (4, 2) is 15. So to compute the sum of the elements in the fourth 
column of A, type 

A(1,4) + A(2,4) + A(3,4) + A(4,4) 

This produces 

ans = 
34 

but is not the most elegant way of summing a single column. 

It is also possible to refer to the elements of a matrix with a single subscript, 
A(k). This is the usual way of referencing row and column vectors. But it can 
also apply to a fully two-dimensional matrix, in which case the array is 
regarded as one long column vector formed from the columns of the original 
matrix. So, for our magic square, A(8) is another way of referring to the valué 
15 stored in A (4, 2). 

If you try to use the valué of an element outside of the matrix, it is an error. 

t = A(4,5) 
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Index exceeds matrix dimensions. 

On the other hand, if you store a valué in an element outside of the matrix, the 
size increases to accommodate the newcomer. 



= A; 










(4,5) = 


= 17 








16 


3 


2 


13 





5 


10 


11 


8 





9 


6 


7 


12 





4 


15 


14 


1 


17 



The Colon Operator 

The colon, :, is one of the most important MATLAB operators. It occurs in 
several different forms. The expression 

1 :10 

is a row vector containing the integers from 1 to 10, 

1 23456789 10 

To obtain nonunit spacing, specify an increment. For example, 
100:-7:50 

is 

100 93 86 79 72 65 58 51 

and 

0:pi/4:pi 

is 

0.7854 1.5708 2.3562 3.1416 

Subscript expressions involving colons refer to portions of a matrix. 
A(1:k,j) 

is the first k elements of the jth column of A. So 
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sum(A(1 :4,4)) 

computes the sum of the fourth column. But there is a better way. The colon by 
itself refers to all the elements in a row or column of a matrix and the keyword 
end refers to the last row or column. So 

sum(A( : ,end) ) 

computes the sum of the elements in the last column of A. 

ans = 
34 

Why is the magic sum for a 4-by-4 square equal to 34? If the integers from 1 to 
16 are sorted into four groups with equal sums, that sum must be 

sum(1 :16)/4 

which, of course, is 

ans = 
34 

The magic Functíon 

MATLAB actually has a built-in function that creates magic squares of almost 
any size. Not surprisingly, this function is named magic. 

B = magic(4) 



B = 












16 


2 


3 


13 




5 


11 


10 


8 




9 


7 


6 


12 




4 


14 


15 


1 



This matrix is almost the same as the one in the Dürer engraving and has all 
the same "magic" properties; the only difference is that the two middle columns 
are exchanged. 

To make this B into Dürer's A, swap the two middle columns. 

A = B( : , [1 3 2 4]) 
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This says, for each of the rows of matrix B, reorder the elements in the order 1, 
3, 2, 4. It produces 



A = 












16 


3 


2 


13 




5 


10 


11 


8 




9 


6 


7 


12 




4 


15 


14 


1 



Why would Dürer go to the trouble of rearranging the columns when he could 
have used MATLAB ordering? No doubt he wanted to include the date of the 
engraving, 1514, at the bottom of his magic square. 
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Like most other programming languages, MATLAB provides mathematical 
expressions, but unlike most programming languages, these expressions 
involve entire matrices. The building blocks of expressions are 

• "Variables" on page 4-10 

• "Numbers" on page 4-11 

• "Operators" on page 4-11 

• "Functions" on page 4-12 

See also "Examples of Expressions" on page 4-13. 

Variables 

MATLAB does not require any type declarations or dimensión statements. 
When MATLAB encounters a new variable ñame, it automatically creates the 
variable and allocates the appropriate amount of storage. If the variable 
already exists, MATLAB changes its contents and, if necessary, allocates new 
storage. For example, 

num_students = 25 

creates a 1-by-l matrix named num_students and stores the valué 25 in its 
single element. To view the matrix assigned to any variable, simply enter the 
variable ñame. 

Variable ñames consist of a letter, followed by any number of letters, digits, or 
underscores. MATLAB is case sensitive; it distinguishes between uppercase 
and lowercase letters. A and a are not the same variable. 

Although variable ñames can be of any length, MATLAB uses only the first N 
characters of the ñame, (where N is the number returned by the function 
namelengthmax), and ignores the rest. Henee, it is important to make each 
variable ñame unique in the first N characters to enable MATLAB to 
distinguish variables. 

N = namelengthmax 
N = 

63 
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The genvarname function can be useful in creating variable ñames that are 
both valid and unique. 

Numbers 

MATLAB uses conventional decimal notation, with an optional decimal point 
and leading plus or minus sign, for numbers. Scientific notation uses the letter 
e to specify a power-of-ten scale factor. Imaginary numbers use either i or j as 
a suffix. Some examples of legal numbers are 

3 -99 0.0001 

9.6397238 1.60210e-20 6.02252e23 

1i -3.14159] 3e5i 

All numbers are stored internally using the long format specified by the IEEE 

floating-point standard. Floating-point numbers have a finite precisión of 

roughly 16 significant decimal digits and a finite range of roughly 10" 308 to 
-j^q+308 

Operators 

Expressions use familiar arithmetic operators and precedence rules. 

+ Addition 

Subtraction 

* Multiplication 

/ División 

\ Left división (described in "Matrices and Linear 

Algebra" in the MATLAB documentation) 

Power 

Complex conjúgate transpose 

( ) Specify evaluation order 
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Functions 

MATLAB provides a large number of standard elementary mathematical 
functions, including abs, sqrt, exp, and sin. Taking the square root or 
logarithm of a negative number is not an error; the appropriate complex result 
is produced automatically. MATLAB also provides many more advanced 
mathematical functions, including Bessel and gamma functions. Most of these 
functions accept complex arguments. For a list of the elementary mathematical 
functions, type 

help elfun 

For a list of more advanced mathematical and matrix functions, type 

help specfun 
help elmat 

Some of the functions, like sqrt and sin, are built in. Built-in functions are 
part of the MATLAB core so they are very efficient, but the computational 
details are not readily accessible. Other functions, like gamma and sinh, are 
implemented in M-files. 

There are some differences between built-in functions and other functions. For 
example, for built-in functions, you cannot see the code. For other functions, 
you can see the code and even modify it if you want. 

Several special functions provide valúes of useful constants. 

pi 3.14159265... 

i Imaginary unit, J-l 

j Same as i 

—52 
eps Floating-point relative precisión, £ = 2 

-1022 
realmin Smallest floating-point number, 2 

1023 
realmax Largest floating-point number, (2 - £)2 

Inf Infinity 

NaN Not-a-number 



4-12 



Expressions 



Infinity is generated by dividing a nonzero valué by zero, or by evaluating well 
defined mathematical expressions that overflow, i.e., exceed realmax. 
Not-a-number is generated by trying to evalúate expressions like 0/0 or 
I nf - 1 nf that do not have well defined mathematical valúes. 

The function ñames are not reserved. It is possible to overwrite any of them 
with a new variable, such as 

eps = 1 .e-6 

and then use that valué in subsequent calculations. The original function can 
be restored with 

clear eps 

Examples of Expressions 

You have already seen several examples of MATLAB expressions. Here are a 
few more examples, and the resulting valúes. 

rho = (1+sqrt(5) )/2 
rho = 

1 .6180 

a = abs(3+4i) 
a = 

5 

z = sqrt (besselk(4/3, rho-i) ) 
z = 

0.3730+ 0.3214Í 

huge = exp (log( realmax ) ) 
huge = 

1 .7977e+308 

toobig = pi*huge 
toobig = 
Inf 
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Workíng with Matrices 

This section introduces you to other ways of creating matrices: 

• "Generating Matrices" on page 4-14 

• "The load Function" on page 4-15 

• "M-Files" on page 4-15 

• "Concatenation" on page 4-16 

• "Deleting Rows and Columns" on page 4-17 

Generating Matrices 

MATLAB provides four functions that genérate basic matrices. 

zeros All zeros 

ones All ones 

rand Uniformly distributed random elements 

randn Normally distributed random elements 



re 


are some examples. 


Z 


= zeros(2,4) 




Z 


= 




















F 


= 5*ones(3,3) 




F 


= 






5 5 


5 




5 5 


5 




5 5 


5 


N 


= f ix(10*rand(- 


1,10)) 


N 


= 






9 2 


6 4 


R 


= randn(4,4) 
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R 


= 














0.6353 


0.0860 


-0 


.3210 


-1 


.2316 




-0.6014 


-2.0046 


1 


.2366 


1 


.0556 




0.5512 


-0.4931 


-0 


.6313 


-0 


.1132 




-1 .0998 


0.4620 


-2 


.3252 





.3792 



The load Functíon 

The load function reads binary files containing matrices genera ted by earlier 
MATLAB sessions, or reads text files containing numeric data. The text file 
should be organized as a rectangular table of numbers, separated by blanks, 
with one row per line, and an equal number of elements in each row. For 
example, outside of MATLAB, créate a text file containing these four lines. 



16.0 


3.0 


2.0 


13.0 


5.0 


10.0 


11 .0 


8.0 


9.0 


6.0 


7.0 


12.0 


4.0 


15.0 


14.0 


1 .0 



Store the file under the ñame magik . dat. Then the statement 
load magik.dat 

reads the file and creates a variable, magik, containing our example matrix. 

An easy way to read data into MATLAB in many text or binary formats is to 
use Import Wizard. 

M-Fíles 

You can créate your own matrices using M-/ü/es, which are text files containing 
MATLAB code. Use the MATLAB Editor or another text editor to créate a file 
containing the same statements you would type at the MATLAB command 
line. Save the file under a ñame that ends in . m. 

For example, créate a file containing these five lines. 



A = [ 








16.0 


3.0 


2.0 


13.0 


5.0 


10.0 


11 .0 


8.0 


9.0 


6.0 


7.0 


12.0 


4.0 


15.0 


14.0 


1.0 ]; 
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Store trie file under the ñame magik . m. Then the statement 

magik 

reads the file and creates a variable, A, containing our example matrix. 



Concatenatíon 

Concatenation is the process of joining small matrices to make bigger ones. In 
fact, you made your first matrix by concatenating its individual elements. The 
pair of square brackets, [ ] , is the concatenation operator. For an example, start 
with the 4-by-4 magic square, A, and form 

B = [A A+32; A+48 A+16] 

The result is an 8-by-8 matrix, obtained by joining the four submatrices. 
B = 



16 


3 


2 


13 


48 


35 


34 


45 


5 


10 


11 


8 


37 


42 


43 


40 


9 


6 


7 


12 


41 


38 


39 


44 


4 


15 


14 


1 


36 


47 


46 


33 


64 


51 


50 


61 


32 


19 


18 


29 


53 


58 


59 


56 


21 


26 


27 


24 


57 


54 


55 


60 


25 


22 


23 


28 


52 


63 


62 


49 


20 


31 


30 


17 



This matrix is halfway to being another magic square. Its elements are a 
rearrangement of the integers 1 :64. Its column sums are the correct valué for 
an 8-by-8 magic square. 

sum(B) 

ans = 

260 260 260 260 260 260 260 260 

But its row sums, sum(B ' ) ' , are not all the same. Further manipulation is 
necessary to make this a valid 8-by-8 magic square. 
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Deletíng Rows and Columns 

You can delete rows and columns from a matrix using just a pair of square 
brackets. Start with 



X = A; 






Then, to delete the second column of X, use 


X(:,2) = 


[] 




This changes X to 




X = 






16 


2 


13 


5 


11 


8 


9 


7 


12 


4 


14 


1 



If you delete a single element from a matrix, the result is not a matrix anymore. 
So, expressions like 

X(1,2) = [] 

result in an error. However, using a single subscript deletes a single element, 
or sequence of elements, and reshapes the remaining elements into a row 
vector. So 

X(2:2:10) = [] 

results in 

X = 

16 9 2 7 13 12 1 



4-17 



1 Matrices and Ar 



rays 



More About Matrices and Arrays 



This section shows you more about working with matrices and arrays, focusing 
on 

• "Linear Algebra" on page 4-18 

• "Arrays" on page 4-21 

• "Multivariate Data" on page 4-24 

• "Scalar Expansión" on page 4-25 

• "Logical Subscripting" on page 4-26 

• "The find Function" on page 4-27 

Linear Algebra 

Informally, the terms matrix and array are often used interchangeably. More 
precisely, a matrix is a two-dimensional numeric array that represents a linear 
transformation . The mathematical operations defined on matrices are the 
subject oí linear algebra. 

Dürer's magic square 

A = [16 3 2 13 

5 10 11 8 

9 6 7 12 

4 15 14 1 ] 

provides several examples that give a taste of MATLAB matrix operations. You 
have already seen the matrix transpose, A'. Adding a matrix to its transpose 
produces a symmetric matrix. 

A + A' 



ns = 








32 


8 


11 


17 


8 


20 


17 


23 


11 


17 


14 


26 


17 


23 


26 


2 



The multiplication symbol, *, denotes the matrix multiplication involving inner 
producís between rows and columns. Multiplying the transpose of a matrix by 
the original matrix also produces a symmetric matrix. 
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A'*A 

ans = 

378 212 206 360 

212 370 368 206 

206 368 370 212 

360 206 212 378 



The determinant of this particular matrix happens to be zero, indicating that 
the matrix is singular. 

d = det(A) 





The reduced row echelon form of A is not the identity. 
R = rref(A) 

R = 

10 1 

10-3 

13 



Since the matrix is singular, it does not have an inverse. If you try to compute 
the inverse with 

X = inv(A) 

you will get a warning message 

Warning: Matrix is cióse to singular or badly scaled. 

Results may be inaccurate. RC0ND = 9.796086e-018. 

Roundoff error has prevented the matrix inversión algorithm from detecting 
exact singularity. But the valué of rcond, which stands for reciprocal condition 
estímate, is on the order of eps, the floating-point relative precisión, so the 
computed inverse is unlikely to be of much use. 
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The eigenvalues of the magic square are interesting. 
e = eig(A) 



34.0000 

8.0000 

0.0000 
-8.0000 

One of the eigenvalues is zero, which is another consequence of singularity. 
The largest eigenvalue is 34, the magic sum. That is because the vector of all 
ones is an eigenvector. 

v = ones(4,1 ) 



1 
1 
1 
1 

A*v 

ans = 
34 
34 
34 
34 

When a magic square is scaled by its magic sum, 
P = A/34 

the result is a doubly stochastic matrix whose row and column sums are all 1. 

P = 






.4706 





.0882 





.0588 





.3824 





.1471 





.2941 





.3235 





.2353 





.2647 





.1765 





.2059 





.3529 





.1176 





.4412 





.4118 





.0294 
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Such matrices represent the transition probabilities in a Markov process. 
Repeated powers of the matrix represent repeated steps of the process. For our 
example, the fifth power 

P~5 
is 






,2507 





.2495 





.2494 





.2504 





.2497 





.2501 





.2502 





.2500 





.2500 





.2498 





.2499 





.2503 





.2496 





.2506 





.2505 





.2493 



This shows that as k approaches infinity, all the elements in the k th power, 
p , approach 1/4 . 

Finally, the coefficients in the characteristic polynomial 

poly(A) 

are 

1 -34 -64 2176 

This indicates that the characteristic polynomial 

det(A-?J) 
is 

?i 4 -34A 3 -64A 2 + 2176?i 

The constant term is zero, because the matrix is singular, and the coefficient of 
the cubic term is -34, because the matrix is magic! 

Arrays 

When they are taken away from the world of linear algebra, matrices become 
two-dimensional numeric arrays. Arithmetic operations on arrays are done 
element by element. This means that addition and subtraction are the same for 
arrays and matrices, but that multiplicative operations are different. MATLAB 
uses a dot, or decimal point, as part of the notation for multiplicative array 
operations. 
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The list of operators includes 



Addition 

Subtraction 

Element-by-element multiplication 

Element-by-element división 

Element-by-element left división 

Element-by-element power 

Unconjugated array transpose 



If the Dürer magic square is multiplied by itself with array multiplication 
A.*A 

the result is an array containing the squares of the integers from 1 to 16, in an 
unusual order. 

ans = 



256 


9 


4 


169 


25 


100 


121 


64 


81 


36 


49 


144 


16 


225 


196 


1 



Building Tables 

Array operations are useful for building tables. Suppose n is the column vector 

n = (0:9)'; 
Then 

pows = [n n."2 2."n] 
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builds a table of squares and powers of 2. 



pows 









1 


1 


1 


2 


2 


4 


4 


3 


9 


8 


4 


16 


16 


5 


25 


32 


6 


36 


64 


7 


49 


128 


8 


64 


256 


9 


81 


512 



The elementary math functions opérate on arrays element by element. So 

format short g 
x = (1 :0.1 :2) ' ; 
logs = [x log10(x) ] 

builds a table of logarithms. 
logs = 



1 


.0 





1 


.1 


0.04139 


1 


.2 


0.07918 


1 


.3 


0.11394 


1 


.4 


0.14613 


1 


.5 


0.17609 


1 


.6 


0.20412 


1 


.7 


0.23045 


1 


.8 


0.25527 


1 


.9 


0.27875 


2 


.0 


0.30103 
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Mu Itíva ríate Data 

MATLAB uses column-oriented analysis for multivariate statistical data. Each 
column in a data set represents a variable and each row an observation. The 
( i , ] ) th element is the ith observation of the j th variable. 

As an example, consider a data set with three variables: 

• Heart rate 

• Weight 

• Hours of exercise per week 

For five observations, the resulting array might look like 

D = [ 72 134 3.2 

81 201 3.5 
69 156 7.1 

82 148 2.4 
75 170 1.2] 

The first row contains the heart rate, weight, and exercise hours for patient 1, 
the second row contains the data for patient 2, and so on. Now you can apply 
many MATLAB data analysis functions to this data set. For example, to obtain 
the mean and standard deviation of each column, use 

mu = mean(D), sigma = std(D) 

mu = 

75.8 161.8 3.48 

sigma = 

5.6303 25.499 2.2107 

For a list of the data analysis functions available in MATLAB, type 
help datafun 

If you have access to the Statistics Toolbox, type 
help stats 
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Scalar Expansión 

Matrices and scalars can be combined in several different ways. For example, 
a scalar is subtracted from a matrix by subtracting it from each element. The 
average valué of the elements in our magic square is 8.5, so 

B = A - 8.5 

forms a matrix whose column sums are zero. 

B = 

7.5 -5.5 -6.5 4.5 

-3.5 1.5 2.5 -0.5 

0.5 -2.5 -1 .5 3.5 

-4.5 6.5 5.5 -7.5 

sum(B) 

ans = 



With scalar expansión, MATLAB assigns a specified scalar to all Índices in a 
range. For example, 



B(1 :2,2:3) = 


= 






•oes out a portion of B. 






B = 
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Logícal Subscríptíng 

The logical vectors created from logical and relational operations can be used 
to reference subarrays. Suppose X is an ordinary matrix and L is a matrix of the 
same size that is the result of some logical operation. Then X(L) specifies the 
elements of X where the elements of L are nonzero. 

This kind of subscripting can be done in one step by specifying the logical 
operation as the subscripting expression. Suppose you have the following set of 
data. 

x = [2.1 1.7 1.6 1.5 NaN 1.9 1.8 1.5 5.1 1.8 1.4 2.2 1.6 1.8]; 

The NaN is a marker for a missing observation, such as a failure to respond to 
an item on a questionnaire. To remove the missing data with logical indexing, 
use isf inite(x), whichis true for allfinite numerical valúes and false for NaN 
and Inf. 

x = x(isf inite(x) ) 
x = 

2.1 1.7 1.6 1.5 1.9 1.8 1.5 5.1 1.8 1.4 2.2 1.6 1.8 

Now there is one observation, 5.1, which seems to be very different from the 
others. It is an outlier. The following statement removes outliers, in this case 
those elements more than three standard deviations from the mean. 

x = x(abs(x-mean(x) ) <= 3*std(x)) 
x = 

2.1 1.7 1.6 1.5 1.9 1.8 1.5 1.8 1.4 2.2 1.6 1.8 

For another example, highlight the location of the prime numbers in Dürer's 
magic square by using logical indexing and scalar expansión to set the 
nonprimes to 0. (See "The magic Function" on page 4-8.) 



A(-isprime 


(A)) 


= 




A = 
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The fínd Functíon 

The f ind function determines the Índices of array elements that meet a given 
logical condition. In its simplest form, f ind returns a column vector of Índices. 
Transpose that vector to obtain a row vector of índices. For example, start 
again with Dürer's magic square. (See "The magic Function" on page 4-8.) 

k = f ind(isprime(A) ) ' 

picks out the locations, using one-dimensional indexing, of the primes in the 
magic square. 

k = 

2 5 9 10 11 13 

Display those primes, as a row vector in the order determined by k, with 
A(k) 

ans = 

5 3 2 11 7 13 

When you use k as a left-hand-side index in an assignment statement, the 
matrix structure is preserved. 

A(k) = NaN 

A = 

16 NaN NaN NaN 

NaN 10 NaN 8 

9 6 NaN 12 

4 15 14 1 
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Controlling Command Window Input and Output 

So far, you have been using the MATLAB command line, typing functions and 
expressions, and seeing the results printed in the Command Window. This 
section describes 

• "The format Function" on page 4-28 

• "Suppressing Output" on page 4-30 

• "Entering Long Statements" on page 4-30 

• "Command Line Editing" on page 4-30 

The format Function 

The format function controls the numeric format of the valúes displayed by 
MATLAB. The function affects only how numbers are displayed, not how 
MATLAB computes or saves them. Here are the different formats, together 
with the resulting output produced from a vector x with components of 
different magnitudes. 

Note To ensure proper spacing, use a fixed-width font, such as Courier. 

x = [4/3 1 .2345e-6] 
format short 

1.3333 0.0000 
format short e 

1 .3333e+000 1 .2345e-006 
format short g 

1.3333 1.2345e-006 
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format long 

1 .33333333333333 0.00000123450000 
format long e 

1 .333333333333333e+000 1 . 234500000000000e-006 
format long g 

1.33333333333333 1.2345e-006 

format bank 

1.33 0.00 
format rat 

4/3 1/810045 
format hex 

3ff 5555555555555 3eb4b6231 abf d271 

If the largest element of a matrix is larger than 10 3 or smaller than 10" 3 , 
MATLAB applies a common scale factor for the short and long formats. 

In addition to the format functions shown above 

format compact 

suppresses many of the blank lines that appear in the output. This lets you 
view more information on a screen or window. If you want more control over 
the output format, use the sprintf and f printf functions. 
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Suppressíng Output 

If you simply type a statement and press Return or Enter, MATLAB 
automatically displays the results on screen. However, if you end the line with 
a semicolon, MATLAB performs the computation but does not display any 
output. This is particularly useful when you genérate large matrices. For 
example, 

A = magic(100) ; 

Enteríng Long Statements 

If a statement does not fit on one line, use an ellipsis (three periods), . . ., 
followed by Return or Enter to indicate that the statement continúes on the 
next line. For example, 

s = 1 -1/2 + 1/3 -1/4 + 1/5 - 1/6 + 1/7 ... 
- 1/8 + 1/9 - 1/10 + 1/11 - 1/12; 

Blank spaces around the =, +, and - signs are optional, but they improve 
readability. 

Command Une Edítíng 

Various arrow and control keys on your keyboard allow you to recall, edit, and 
reuse statements you have typed earlier. For example, suppose you mistakenly 
enter 

rho = (1 + sqt(5))/2 

You have misspelled sqrt. MATLAB responds with 
Undefined function or variable ' sqt ' . 

Instead of retyping the entire line, simply press the T key. The statement you 
typed is redisplayed. Use the <- key to move the cursor over and insert the 
missing r. Repeated use of the T key recalls earlier lines. Typing a few 
characters and then the T key finds a previous line that begins with those 
characters. You can also copy previously executed statements from the 
Command History. For more information, see "Command History" on page 8-5. 

Following is the list of arrow and control keys you can use in the Command 
Window. If the preference you select for Command line key bindings is 
Emacs (MATLAB standard), you can also use the Ctrl+key combinations 
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shown. See also general keyboard shortcuts for desktop tools in the Desktop 
Tools and Development Environment documentation. 



Key 



Control Key for Emacs 
(MATLAB standard) 
Preference 



Operation 



t 



i 



Ctrl+ -> 

Ctrl+ <- 

Home 

End 

Ctrl+Home 

Ctrl+End 

Esc 

Delete 

Backspace 

Shift+Home 
Shift+End 



Ctrl+P 

Ctrl+N 

Ctrl+B 
Ctrl+F 



Ctrl+A 
Ctrl+E 



Ctrl+U 
Ctrl+D 
Ctrl+H 
Ctrl+K 



Recall previous line. Works only at command line. 

Recall next line. Works only at command line if 
you previously used the up arrow or Ctrl+P. 

Move back one character. 

Move /brward one character. 

Move right one word. 

Move Zeft one word. 



Move to beginning of command line. 
Move to end of command line. 
Move to top of Command Window. 
Move to end of Command Window. 
Clear command line. 



Delete character at cursor in command line. 

Delete character before cursor in command line. 

Cut contents (&Í11) to end of command line. 

Highlight to beginning of command line. 

Highlight to end of last line. Can start at any line 
in the Command Window. 



4-31 



1 Matrices and Arrays 



4-32 




Graphics 



This chapter provides an introduction to plotting data in MATLAB. 

Overview of MATLAB Plotting (p. 5-2) Créate plots, include múltiple data sets, specify property 

valúes, and save figures. 



Editing Plots (p. 5-16) 

Mesh and Surface Plots (p. 5-52) 
Images (p. 5-58) 
Printing Graphics (p. 5-60) 
Handle Graphics (p. 5-62) 
Animations (p. 5-71) 



Edit plots interactively and using functions, and use the 
property editor. 

Visualize functions of two variables. 

Work with images. 

Print and export figures. 

Work with graphics objects and set object properties. 

Créate moving graphics. 



3 Graphics 



Overvíew of MATLAB Plottíng 



MATLAB provides a wide variety of techniques to display data graphically. 
Interactive tools enable you to manipúlate graphs to achieve results that reveal 
the most information about your data. You can also annotate and print graphs 
for presentations, or export graphs to standard graphics formats for 
presentation in web browsers or other media. 



For More Information "Graphics" and "3-D Visualization" in the MATLAB 
documentation provide in-depth coverage of MATLAB graphics and 
visualization tools. Access these topics from the Help browser. 



The Plottíng Process 

The process of visualizing data typically involves a series of operations. This 
section provides a "big picture" view of the plotting process and contains links 
to sections that have examples and specific details about performing each 
operation. 

Creating a Graph 

The type of graph you choose to créate depends on the nature of your data and 
what you want to reveal about the data. MATLAB predefines many graph 
types, such as line, bar, histogram, and pie graphs. There are also 3-D graphs, 
such as surfaces, slice planes, and streamlines. 

There are two basic ways to créate graphs in MATLAB: 

• Use plotting tools to créate graphs interactively. 

See "Examples — Using MATLAB Plotting Tools" on page 5-20. 

• Use the command interface to enter commands in the Command Window or 
créate plotting programs. 

See "Basic Plotting Functions" on page 5-38. 

You might find it useful to combine both approaches. For example, you might 
issue a plotting command to créate a graph and then modify the graph using 
one of the interactive tools. 
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Exploring Data 

Once you créate a graph, you can extract specific information about the data, 
such as the numeric valué of a peak in a plot, the average valué of a series of 
data, or you can perform data fitting. 



For More Information See "Data Exploration Tools" and "Basic Fitting 
Options" in the MATLAB documentation. 



Editing the Graph Components 

Graphs are composed of objects, which have properties you can change. These 
properties affect the way the various graph components look and behave. 

For example, the axes used to define the coordínate system of the graph has 
properties that define the limits of each axis, the scale, color, etc. The line used 
to créate a line graph has properties such as color, type of marker used at each 
data point (if any), line style, etc. 

Note that the data used to créate a line graph are properties of the line. You 
can, therefore, change the data without actually creating a new graph. 

See "Editing Plots" on page 5-16. 

Annotating Graphs 

Annotations are the text, arrows, callouts, and other labels added to graphs to 
help viewers see what is important about the data. You typically add 
annotations to graphs when you want to show them to other people or when you 
want to save them for later reference. 



For More Information See "Annotating Graphs" in the MATLAB 
documentation or select Annotating Graphs from the figure Help menú. 
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Printing and Exporting Graphs 

You can print your graph on any printer connected to your computer. The print 
previewer enables you to view how your graph will look when printed. It 
enables you to add headers, footers, a date, and so on. The page setup dialog 
lets you control the size, layout, and other characteristics of the graph (select 
Page Setup from the figure File menú). 

Exporting a graph means creating a copy of it in a standard graphics file 
format, such as TIF, JPEG, or EPS. You can then import the file into a word 
processor, include it in an HTML document, or edit it in a drawing package 
select Export Setup from the figure File menú). 



For More Information See the print command reference page and 
"Printing and Exporting" in the MATLAB documentation or select Printing 
and Exporting from the figure Help menú. 



Saving Graphs to Reload into MATLAB 

There are two ways to save graphs that enable you to save the work you have 
invested in their preparation: 

• Save the graph as a FIG-file (select Save from the figure File menú). 

• Genérate MATLAB code that can recréate the graph (select Genérate 
M-File from the figure File menú). 

FIG-Files. FIG-files are a binary format that saves a figure in its current state. 
This means that all graphics objects and property settings are stored in the file 
when you créate it. You can reload the file into a different MATLAB session, 
even if you are running MATLAB on a different type of computer. When you 
load a FIG-file, MATLAB creates a new figure in the same state as the one you 
saved. 

Note that the states of any figure tools (i.e., any items on the toolbars) are not 
saved in a FIG-file; only the contents of the graph are saved. 
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Genemted Code. You can use the MATLAB M-code generator to créate code that 
recreates the graph. Unlike a FIG-file, the generated code does not contain any 
data. You must pass the data to the generated function when you run the code. 

Note that studying the generating code for a graph is a good way to learn how 
to program with MATLAB. 



For More Information See the print command reference page and "Saving 
Your Work" in the MATLAB documentation. 



Graph Components 

MATLAB displays graphs in a special window known as a figure. To créate a 
graph, you need to define a coordinate system. Therefore every graph is placed 
within axes, which are contained by the figure. 

The actual visual representation of the data is achieved with graphics objects 
like lines and surfaces. These objects are drawn within the coordinate system 
defined by the axes, which MATLAB automatically creates specifically to 
accommodate the range of the data. The actual data is stored as properties of 
the graphics objects. 

See "Handle Graphics" on page 5-62 for more information about graphics object 
properties. 

The following picture shows the basic components of a typical graph. 
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Figure Tools 

The figure is equipped with sets of tools that opérate on graphs. The figure 
Tools menú provides access to many graph tools. 



-J Figure 1 
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For More Information See "MATLAB Plotting Tools" in the MATLAB 
documentation or select Plotting Tools from the figure Help menú. 
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Access to Tools 

You can access the figure toolbars and the plotting tools from the View menú, 
as shown in the following picture. 
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Figure Toolbars 

Figure toolbars provide easy access to many graph modification features. There 
are three toolbars. When you place the cursor over a particular tool, a text box 
pops up with the tool ñame. The following picture shows the three toolbars 
displayed with the cursor over the Data Cursor tool. 



- J Figure 1 




File Edit View Insert Tools Desktop Window Help 


■íl 


D^yal k|<*^o®|^iJjD0|no 


'< K $ k ¿I # Í \i feata?ursor1 O | % © 


■ |\ W^ t DO|-w 4 



For More Information See "Figure Toolbars" in the MATLAB 
documentation. 
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Plotting Tools 

Plotting tools are attached to figures and créate an environment for creating 
graphs. These tools enable you to do the following: 

• Select from a wide variety of graph types 

• Set the properties of graphics objects 

• Annotate graphs with text, arrows, etc. 

• Créate and arrange subplots in the figure 

• Drag and drop data into graphs 

Display the plotting tools from the View menú or by clicking in the figure 
toolbar, as shown in the following picture. 

Enable plotting tools from the View 
menú or toolbar. 
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There are three components to the plotting tools: 

• Figure Palette — Specify and arrange subplots, access workspace variables 
for plotting or editing, and add annotations. 

• Plot Browser — Select objects in the graphics hierarchy, control visibility, 
and add data to axes. 

• Property Editor — Change key properties of the selected object. Click 
Inspector for access to all object properties. 
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The following picture shows a figure with the plotting tools enabled. 
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Plotting Tools and MATLAB Commands 

You can enable the plotting tools on any graph, even if you created it using 
MATLAB commands. For example, suppose you créate the following graph. 



t = 0:pi/20:2*pi; 
y = exp(sin(t)) ; 
plotyy(t,y,t,y, ' plot ' 
xlabel( 'X Axis ' ) 
ylabel( 'Plot Y Axis' ) 
title( 'Two Y Axes' ) 



stem ' ) 




This graph contains two y -axes, one for each plot type (lineseries and stem 
graphs). The plotting tools make it easy to select any of the objects that the 
graph contains and set their properties. 
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For example, adding a label for they-axis that corresponds to the stem plot is 
easily accomplished by selecting that axes in the Plot Browser and setting the 
YLabel property in the Property Editor. 
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Arrangíng Graphs Within a Figure 

You can place a number of axes within a figure by selecting the layout you want 
from the Figure Palette. For example, the following picture shows how to 
specify four 2-D axes in the figure. 



Click to add one axes to bottom 
of current layout. 

Click and drag right lo specify axes layout. 




O 



DD 

EB 



□DDD 



\ Line 
\ Arrow 
^ Double Arrow 
X Text Arrow 
T Text Box 
I I Rectangle 
O Ellipse 



nnnn 



0.6 
0.4 



1 = 

0.2 
.. 1 n _ 




Select the axes you want to target for plotting. You can also use the subplot 
command to créate múltiple axes. 
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Selecting Plot Types 

You can use the Plot Catalog to select from a variety of techniques for plotting 
data. To access the Plot Catalog, 

1 Select the variables you want to plot in the Figure Palette. 

2 Right-click to display the context menú. 

3 Select More Plots to display the catalog. 
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MATLAB displays the Plot Catalog with the selected variables ready to plot, 
once you select a plot type. 



Specify variables to plot. 



See a description of etich plot type. 
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Description 




2-D Line Graph 

2-D line graph using linear 
axes 

Vectors créate a single line; 
matrices créate one line per 
column. 



Plotted Variables 

* Single variable -- plot a 
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Select a calegory of graphs and then choose a specific lype. 
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MATLAB automatically formats graphs by setting the scale of the axes, adding 
tick marks along axes, and using colors and line styles to distinguish the data 
plotted in the graph. However, if you are creating graphics for presentation, 
you can change the default formatting or add descriptive labels, titles, legends 
and other annotations to help explain your data. 

Plot Edítíng Mode 

Plot editing mode enables you to perform point-and-click editing of the graphics 
objects in your graph. 

Enabling Plot Edit Mode 

To enable plot edit mode, click the arrowhead in the figure toolbar: 

Enable plot edit mode. 
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You can also select Edit Plot from the figure Tools menú. 

Setting Object Properties 

Once you have enabled plot edit mode, you can select objects by clicking on 
them in the graph. Selection handles appear and indicate that the object is 
selected. Select múltiple objects using Shlft+cllck. 

Right-click with the pointer over the selected object to display the object's 
context menú: 
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The context menú provides quick access to the most commonly used operations 
and properties. 



Using the Property Editor 

In plot edit mode, double-clicking on an object in a graph starts the Property 
Editor with that object's major properties displayed. The Property Editor 
provides access to the most used object properties. It is updated to display the 
properties of whatever object you select. 
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Click to display Property Inspector. 
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Accessing All Properties — Property Inspector 

The Property Inspector is a tool that enables you to access all object properties. 
If you do not find the property you want to set in the Property Editor, click the 
Inspector button to display the Property Inspector. You can also use the 
inspect command to start the Property Inspector. 

The following picture shows the Property Inspector displaying the properties of 
a graph's axes. It lists each property and provides a text field or other 
appropriate device (such as a color picker) from which you can set the valué of 
the property. 

As you select different objects, the Property Inspector is updated to display the 
properties of the current object. 
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Usíng Functíons to Edít Graphs 

If you prefer to work from the MATLAB command line, or if you are creating 
an M-file, you can use MATLAB commands to edit the graphs you créate. You 
can use the set and get commands to change the properties of the objects in a 
graph. For more information about using graphics commands, see "Handle 
Graphics" on page 5-62. 
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Examples — Using MATLAB Plotting Tools 



Suppose you want to graph the function y = x 3 over the x domain -1 to 1. The 
first step is to genérate the data to plot. 

It is simple to evalúate a function because MATLAB can distribute arithmetic 
operations over all elements of a multivalued variable. 

For example, the following statement creates a variable x that contains valúes 
ranging from -1 to 1 in increments of 0.1 (you could also use the linspace 
function to genérate data for x). The second statement raises each valué in x to 
the third power and stores these valúes in y. 

x = -1:.1:1; % Define the range of x 

y = x."3; % Raise each element in x to the third power 

Now that you have generated some data, you can plot it using the MATLAB 
plotting tools. To start the plotting tools, type 

plottools 
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MATLAB displays a figure with plotting tools attached. 
Variables ¡n workspace Figure plotting área 
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Plotting Two Variables 

A simple line graph is a suitable way to display x as the independent variable 
and y as the dependent variable. To do this, select both variables (click to 
select, then Shift-click to select again), then right-click to display the context 
menú. 
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Select plot(x,y) from the menú. MATLAB creates the line graph in the figure 
área. The black squares indicate that the line is selected and you can edit its 
properties with the Property Editor. 
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Changing the Appearance 

Next change the line properties so that the graph displays only the data point. 
Use the Property Editor to set following properties: 

• Line to no line 

• Marker to o 

• Marker size to 4.0 

• Marker fill color to red 
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--Set Marker f 

color to red. 
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Adding More Data to the Graph 

You can add more data to the graph by defining more variables or by specifying 
an expression that MATLAB uses to genérate data for the plot. This second 
approach makes it easy to explore variations of the data already plotted. 

To add data to the graph, select the axes in the Plot Browser and click the Add 
Data button. When you are using the plotting tools, MATLAB always adds 
data to the existing graph, instead of replacing the graph, as it would if you 
issued repeated plotting commands. That is, the plotting tools are in a hold on 
state. 




Plot Browser 




. Add Data to Axes 



1. Select axes. 



*l 



Plot Type: 
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plot 
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an 


espress 


on 





OK I Cancel 



-3. Enter expression 



2. Click Add Data. 



Add Data.. 



The picture above shows how to configure the Add Data to Axes dialog to 
créate a line plot of y = x 4 , which is added to the existing plot of y = x 3 . 
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Add Data... 



Changing the Type of Graph 

The plotting tools enable you to easily view your data with a variety of plot 
types. The following picture shows the same data as above converted to stem 
plots. To change the plot type, 

1 Select the plotted data in the Plot Browser. 

2 Select Stem in the Plot Type menú. 
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Selecl boih sets of data. 



Select Stem as the Plot 
Type. 
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Modifying the Graph Data Source 

You can link graph data to variables in your workspace. When you change the 
valúes contained in the variables, you can then update the graph to use the new 
data without having to créate a new graph. (See also the ref resh function.) 

Suppose you have the following data: 

x = linspace( -pi,pi,25) ; % 25 points between -71 and n 
y = sin(x) ; 

Using the plotting tools, créate a graph of y = sin(x). 
plottools 




ATLAB sets the data source lo the 
selected variables. 
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New Valúes for the Data Source 

The data that defines the graph is linked to the x and y variables in the base 
workspace. If you assign new valúes to these variables and click the Refresh 
Data button, MATLAB updates the graph to use the new data. 

x = linspace( -2*pi,2*pi,25) ; % 25 points between -2pi and 2pi 
y = sin(x); % Recalculate y based on the new x valúes 
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Click Refresh Data to 

update the plot. 
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Preparíng Graphs for Presentaron 



Suppose you plot the following data and want to créate a graph that presents 
certain information about the data. 

x = -10: .005:40; 

y = [1 .5*cos(x)+4*exp( - .01 *x) . *cos(x)+exp( .07*x) . *sin(3*x) ] ; 

plot(x.y) 

This picture shows the graph created by the code above. 




-10 -5 5 10 15 20 25 30 35 40 



Now suppose you want to save copies of the graph by 

• Printing the graph on a local printer so you have a copy for your notebook 

• Exporting the graph to an Encapsulated PostScript (EPS) file to incorpórate 
into a word processor document 
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Modífy the Graph to Enhance the Presentation 

To obtain a better view, zoom in on the graph using horizontal zoom. 

After enabling zoom mode from the figure toolbar, right-click to display the 
context menú. Select Horizontal Zoom (2-D Plots Only) from the Zoom 
Options. 



~¿ Figure 1 
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Zoom Out Alt-Click 
Reset to Original View 
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Y 



Left-click to zoom in on a región of the graph and use the panning tool to 
position the points of interest where you want them on the graph. 

Label some key points using data tips. 
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"} Figure 1 
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Delete 



Export Cursor Data to Workspace, , 




Finally, add text annotations, axis labels, and a title. You can add the title and 
axis labels using the following commands: 

title ('y = 1.5cos(x) + 4e"{-0.01x}cos(x) + e"{0.07x}sin(3x) ' ) 
xlabel( 'X Axis ' ) 
ylabel( 'Y Axis ' ) 

Note that the text string passed to the title command uses T E X syntax to 
produce the exponents. See the text String property for more information on 
using T E X syntax to produce mathematical symbols. 

The graph is now ready to print and export. 
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Printíng the Graph 

Before printing the graph, use the print previewer to see how the graph will be 
laid out on the page. Display the graph in the print previewer by selecting 
Print Preview from the figure File menú. 

• Click the Header button to add some descriptive text to the top of the page. 
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• Next click the Page Setup button. The Page Setup dialog enables you to set 
a number of properties that control how the page is printed. 

Note that MATLAB recalculates the valúes of the axes tick marks because 
the printed graph is larger than the one displayed on the computer screen. 
To forcé MATLAB to use the same tick marks and limits, select Keep screen 
limits and ticks from the Axes and Figure tab in the Page Setup dialog. 

• Click OK to accept the setting and dismiss the dialog. 



Page Setup - Figure 1 
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• Click Print on the print previewer to send the graph to your default printer. 

The Page Setup dialog provides many other options for controlling how printed 
graphs look. Click Help for more information. 
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Exporting the Graph 



Exporting a graph is the process of creating a standard graphics file format of 
the graph (such as EPS or TIFF), which you can then import into other 
applications like word processors, drawing packages, etc. 

This example exports the graph as an EPS file with the following 
requirements: 

• The size of the picture when imported into the word processor document 
should be four inches wide and three inches high. 

• All the text in the figure should have a size of 8 points. 

Specifying the Size of the Graph 

To set the size, select Export Setup from the figure File menú. 
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Specifying the Font Size 

To set the font size of all the text in the graph, select Fonts in the Export 
Setup dialog's Properties selector. Then select Use fixed font size and enter 
8 in the text box. 
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Selecting the File Format 

Once you finish setting options for the exported graph, click the Export button. 
MATLAB displays a Save As dialog that enables you to specify a ñame for the 
file as well as select the type of file format you want to use. 

For this example, select EPS, as shown in the following picture. 




File ñame: 



Save asíype: 



MATLAB Figure p.fig] 
Adobe lllustrator file (".ai 
Bitimap file f".brni 




Select EPS from the 

Enhanced nnetafile (".emf) % drop-down menú. 

JPEG image f".jpg] 

MATLAB Figure [*.f¡g] 

Paintbrush 24-bit file [*.pcx] 

Portable Bitrmap file ("pbm] 

Portable Document Format (".pdf) 

You can import the saved file into any application that supports EPS files. The 
Save as type drop-down menú lists other options for file types. 

You can also use the print command to print figures on your local printer or to 
export graphs to standard file types. 



For More Information See the print command reference page and 
"Printing and Exporting" in the MATLAB documentation or select Printing 
and Exporting from the figure Help menú. 
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Basic Plottíng Functíons 



This section describes important graphics functions and provides examples of 
some typical applications. The plotting tools, described in previous sections, 
make use of MATLAB plotting functions and use these functions to genérate 
code for graphs. 

"Creating a Plot" on page 5-38 

"Múltiple Data Sets in One Graph" on page 5-40 

"Specifying Line Styles and Colors" on page 5-41 

"Plotting Lines and Markers" on page 5-41 

"Imaginary and Complex Data" on page 5-43 

"Adding Plots to an Existing Graph" on page 5-44 

"Figure Windows" on page 5-46 

"Múltiple Plots in One Figure" on page 5-46 

"Controlling the Axes" on page 5-48 

"Axis Labels and Titles" on page 5-49 

"Saving Figures" on page 5-51 



Creating a Plot 

The plot function has different forms, depending on the input arguments. If y 
is a vector, plot (y) produces a piecewise linear graph of the elements of y 
versus the index of the elements of y. If you specify two vectors as arguments, 
plot (x, y) produces a graph of y versus x. 

For example, these statements use the colon operator to créate a vector of x 
valúes ranging from to 2n, compute the sine of these valúes, and plot the 
result. 

x = 0:pi/100:2*pi; 
y = sin(x) ; 
plot(x,y) 
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Now label the axes and add a title. The characters \pi créate the symbol n. See 
the text string property for more symbols. 

xlabel( x = 0:2\pi' ) 

ylabel( ' Sine of x' ) 

title('Plot of the Sine Function ' , ' FontSize ' , 12) 





Plot of the Sine Function 




1 






i V 


1 


0.8 


~ / \ 


- 


0.6 


- / \ 


- 


0.4 


- / \ 


- 


0.2 


-/ \ 


- 


X 


/ \ 










o 






o> 


\ 


/ 


c 


\ 




co 


\ 




-0.2 


\ 


/ 


-0.4 


\ 


/ 


-0.6 


\ / 


- 


-0.8 

_1 


- , , , x y 


- 



3 4 
x = 0:2n 



5-39 



3 Graphics 



Múltiple Data Sets ¡n One Graph 

Múltiple x-y pair arguments créate múltiple graphs with a single cali to plot. 
MATLAB automatically cycles through a predefined (but user settable) list of 
colors to allow discrimination among sets of data. See the axes ColorOrder and 
LineStyleOrder properties. 

For example, these statements plot three related functions of x, with each 
curve in a sepárate distinguishing color. 

x = 0:pi/100:2*pi; 

y = sin(x) ; 

y2 = sin(x- .25) ; 

y3 = sin(x- .5) ; 

plot(x,y,x,y2,x,y3) 

The legend command provides an easy way to identify the individual plots. 
legend( 'sin(x) ' , ' sin(x- .25) ' , ' sin(x- .5) ' ) 
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For More Information See "Defining the Color of Lines for Plotting" in 
"Axes Properties" in the MATLAB documentation. 



Specífyíng Une Styles and Colors 

It is possible to specify color, line styles, and markers (such as plus signs or 
circles) when you plot your data using the plot command. 

plot (x,y, ' color_style_marker' ) 

color_style_marker is a string containing from one to four characters 
(enclosed in single quotation marks) constructed from a color, a line style, and 
a marker type: 

• Color strings are ' c ' , ' m ' , 'y', ' r ' , ' g ' , ' b ' , ' w ' , and ' k ' . These correspond 
to cyan, magenta, yellow, red, green, blue, white, and black. 

• Line style strings are ' - ' for solid, ' - - ' for dashed, ' : ' for dotted, ' - . ' for 
dash-dot. Omit the line style for no line. 

• The marker types are ' + ', 'o', '*', and ' x ' , and the filled marker types are 
' s ' for square, ' d ' for diamond, ' " ' for up triangle, ' v ' for down triangle, 

' > ' for right triangle, ' < ' for left triangle, ' p ' for pentagram, ' h ' for 
hexagram, and none for no marker. 

You can also edit color, line style, and markers interactively. See "Editing 
Plots" on page 5-16 for more information. 

Plotting Lines and Markers 

If you specify a marker type but not a line style, MATLAB draws only the 
marker. For example, 

plot(x,y, 'ks') 

plots black squares at each data point, but does not connect the markers with 
a line. 
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The statement 

plot(x,y, ' r:+") 

plots a red dotted line and places plus sign markers at each data point. 

Placing Markers at Every Tenth Data Point 

You might want to use fewer data points to plot the markers than you use to 
plot the lines. This example plots the data twice using a different number of 
points for the dotted line and marker plots. 

x1 = 0:pi/100:2*pi; 

x2 = 0:pi/10:2*pi; 

plot (x1 , sin(x1 ) , ' r : ' ,x2,sin(x2) , ' r+' ) 
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Imagínary and Complex Data 

When the arguments to plot are complex, the imaginary part is ignored except 
when you pass plot a single complex argument. For this special case, the 
command is a shortcut for a graph of the real part versus the imaginary part. 
Therefore, 

plot(Z) 

where Z is a complex vector or matrix, is equivalent to 
plot (real(Z) ,imag(Z) ) 

For example, 

t = 0:pi/10:2*pi; 
plot(exp(i*t) , ' -o' ) 
axis equal 

draws a 20-sided polygon with little circles at the vértices. The command 
axis equal makes the individual tick-mark increments on the x- and j-axes 
the same length, which makes this plot more circular in appearance. 
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Addíng Plots to an Exístíng Graph 

The hold command enables you to add plots to an existing graph. When you 
type 

hold on 

MATLAB does not replace the existing graph when you issue another plotting 
command; it adds the new data to the current graph, rescaling the axes if 
necessary. 

For example, these statements first créate a contour plot of the peaks function, 
then superimpose a pseudocolor plot of the same function. 

[x,y,z] = peaks; 
pcolor (x,y,z) 
shading interp 
hold on 

contour(x,y,z,20, ' k ' ) 
hold off 

The hold on command causes the pcolor plot to be combined with the contour 
plot in one figure. 
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For More Information See "Creating Specialized Plots" in the MATLAB 
documentation for information on a variety of graph types. 
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Figure Windows 

Graphing functions automatically open a new figure window if there are no 
figure windows already on the screen. If a figure window exists, MATLAB uses 
that window for graphics output. If there are múltiple figure windows open, 
MATLAB targets the one that is designated the "current figure" (the last figure 
used or clicked in). 

To make an existing figure window the current figure, you can click the mouse 
while the pointer is in that window or you can type 

f igure(n) 

where n is the number in the figure title bar. The results of subsequent 
graphics commands are displayed in this window. 

To open a new figure window and make it the current figure, type 

figure 

Clearing the Figure for a New Plot 

When a figure already exists, most plotting commands clear the axes and use 
this figure to créate the new plot. However, these commands do not reset figure 
properties, such as the background color or the colormap. If you have set any 
figure properties in the previous plot, you might want to use the clf command 
with the reset option, 

clf reset 
before creating your new plot to restore the figure's properties to their defaults. 



For More Information See "Figure Properties" and "Graphics Window 
the Figure" in the MATLAB documentation for information on figures. 



Múltiple Plots in One Figure 

The subplot command enables you to display múltiple plots in the same 
window or print them on the same piece of paper. Typing 

subplot (m,n,p) 
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partitions the figure window into an m-by-n matrix of small subplots and selects 
the pth subplot for the current plot. The plots are numbered along first the top 
row of the figure window, then the second row, and so on. For example, these 
statements plot data in four different subregions of the figure window. 

t = 0:pi/10:2*pi; 
[X,Y,Z] = cylinder(4*cos(t) ) ; 
subplot (2, 2,1 ) ; mesh(X) 
subplot(2,2,2) ; mesh(Y) 
subplot(2,2,3) ; mesh(Z) 
subplot(2,2,4) ; mesh(X,Y,Z) 



*) Figure No. 1 
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Controllíng the Axes 

The axis command provides a number of options for setting the scaling, 
orientation, and aspect ratio of graphs. You can also set these options 
interactively. See "Editing Plots" on page 5-16 for more information. 

Setting Axis Limits 

By default, MATLAB finds the máxima and minima of the data and chooses the 
axis limits to span this range. The axis command enables you to specify your 
own limits: 

axis([xmin xmax ymin ymax]) 

or for three-dimensional graphs, 

axis([xmin xmax ymin ymax zmin zmax]) 

Use the command 
axis auto 

to reenable MATLAB automatic limit selection. 

Setting Axis Aspect Ratio 

axis also enables you to specify a number of predefined modes. For example, 

axis square 

makes the x-axis and j-axis the same length. 
axis equal 

makes the individual tick mark increments on the x-axes and j-axes the same 
length. This means 

plot(exp(i*[0:pi/10:2*pi])) 

followed by either axis square or axis equal turns the oval into a proper 
circle. 

axis auto normal 
returns the axis scaling to its default automatic mode. 
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Setting Axis Visibility 

You can use the axis command to make the axis visible or invisible. 

axis on 

makes the axes visible. This is the default. 
axis off 

makes the axes invisible. 

Setting Grid Lines 

The grid command toggles grid lines on and off. The statement 

grid on 

turns the grid lines on, and 

grid off 

turns them back off again. 



For More Information See the axis and axes reference pages and "Axes 
Properties" in the MATLAB documentation. 



Axis La be I s and Tifies 

The xlabel, ylabel, and zlabel commands add x-, y-, and z-axis labels. The 
title command adds a title at the top of the figure and the text function 
inserís text anywhere in the figure. 
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You can produce mathematical symbols using LaTeX notation in the text 
string, as the following example illustrates. 

t = -pi:pi/100:pi; 

y = sin(t) ; 

plot(t,y) 

axis( [ -pi pi -1 1 ] ) 

xlabel('-\pi \leq {\itt} \leq \pi') 

ylabel( ' sin(t) ' ) 

title( 'Graph of the sine function') 

text (1 , -1 /3, ' {\itNote the odd symmetry . } ' ) 



Graph of the sine function 




You can also set these options interactively. See "Editing Plots" on page 5-16 
for more information. 

Note that the location of the text string is defined in axes units (i.e., the same 
units as the data). See the annotation function for a way to place text in 
normalized figure units. 
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Savíng Figures 

Save a figure by selecting Save from the File menú to display a file save dialog. 
MATLAB saves the data it needs to recréate the figure and its contents (i.e., 
the entire graph) in a file with a .f ig extensión. 

To save a figure using a standard graphics format, such as TIFF, for use with 
other applications, select Export Setup from the File menú. You can also save 
from the command line — use the saveas command, including any options to 
save the figure in a different format. 

See "Exporting the Graph" on page 5-35 for an example. 

Saving Workspace Data 

You can save the variables in your workspace using the Save Workspace As 
item in the figure File menú. You can reload saved data using the Import Data 
item in the figure File menú. MATLAB supports a variety of data file formats, 
including MATLAB data files, which have a .mat extensión. 

Generating M-Code to Recréate a Figure 

You can genérate MATLAB code that recreates a figure and the graph it 
contains by selecting the Genérate M-File item from the figure File menú. 
This option is particularly useful if you have developed a graph using plotting 
tools and want to créate a similar graph using the same or different data. 

Saving Figures That Are Compatible with Previous Versión of MATLAB 

Créate backward-compatible FIG-files by following these two steps. 

• Ensure that any plotting functions used to créate the contents of the figure 
are called with the ' v6 ' argument, where applicable. 

• Use the ' - v6 ' option with the hgsave command. 



For More Information See "Plot Objects and Backward Compatibility" in 
the MATLAB documentation more information. 
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Mesh and Surface Plots 



MATLAB defines a surface by the z-coordinates of points above a grid in the x-y 
plañe, using straight lines to connect adjacent points. The mesh and surf 
plotting functions display surfaces in three dimensions. mesh produces 
wireframe surfaces that color only the lines connecting the defining points. 
surf displays both the connecting lines and the faces of the surface in color. 

The figure colormap and figure properties determine how MATLAB colors the 
surface. 

Vísualízíng Functions of Two Variables 

To display a function of two variables, z = f(x,y), 

• Genérate X and Y matrices consisting of repeated rows and columns, 
respectively, over the domain of the function. 

• Use X and Y to evalúate and graph the function. 

The meshgrid function transforms the domain specified by a single vector or 
two vectors x and y into matrices X and Y for use in evaluating functions of two 
variables. The rows of X are copies of the vector x and the columns of Y are 
copies of the vector y. 

Example — Graphing the sinc Function 

This example evaluates and graphs the two-dimensional sinc function, sin(r)/r, 
between thex and j directions. R is the distance from the origin, which is at the 
center of the matrix. Adding eps (a MATLAB command that returns a small 
floating-point number) avoids the indeterminate 0/0 at the origin. 

[X,Y] = meshgrid(-8: .5:8) ; 
R = sqrt(X.~2 + Y. "2) + eps; 
Z = sin(R) ./R; 
mesh(X,Y,Z, 'EdgeColor' , 'black' ) 
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By default, MATLAB colors the mesh using the current colormap. However, 
this example uses a single-colored mesh by specifying the EdgeColor surface 
property. See the surface reference page for a list of all surface properties. 

You can créate a mesh with see-through faces by disabling hidden line removal. 

hidden off 

See the hidden reference page for more information on this option. 

Example — Colored Surface Plots 

A surface plot is similar to a mesh plot except that MATLAB colors the 
rectangular faces of the surface. The color of each faces is determined by the 
valúes of Z and the colormap (a colormap is an ordered list of colors). These 
statements graph the sinc function as a surface plot, specify a colormap, and 
add a color bar to show the mapping of data to color. 

surf (X,Y,Z) 

colormap hsv 
colorbar 
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See the colormap reference page for information on colormaps. 



For More Information See "Creating 3-D Graphs" in the MATLAB 
documentation for more information on surface plots. 



Transparent Surfaces 

You can make the faces of a surface transparent to a varying degree. 
Transparency (referred to as the alpha valué) can be specified for the whole 
object or can be based on an alphamap, which behaves in a way analogous to 
colormaps. For example, 

surf (X,Y,Z) 
colormap hsv 
alpha( .4) 
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produces a surface with a face alpha valué of 0.4. Alpha valúes range from 
(completely transparent) to 1 (not transparent). 




For More Information See "Transparency" in the MATLAB documentation 
for more information on using transparency. 



Surface Plots with Lighting 

Lighting is the technique of illuminating an object with a directional light 
source. In certain cases, this technique can make subtle differences in surface 
shape easier to see. Lighting can also be used to add realism to 
three-dimensional graphs. 

This example uses the same surface as the previous examples, but colors it red 
and removes the mesh lines. A light object is then added to the left of the 
"camera" (the camera is the location in space from where you are viewing the 
surface). 



surf (X,Y,Z, ' FaceColor ' , ' red ' , ' EdgeColor ' 
camlight left; lighting phong 



none 
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Manipulating the Surface 

The figure toolbar and the camera toolbar provide ways to explore 3-D graphics 
interactively. Display the camera toolbar by selecting Camera Toolbar from 
the figure View menú. 

The following picture shows both toolbars with the Rotate 3D tool selected. 



Figure 1 
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These tools enables you to move the camera around the surface object, zoom, 
add lighting, and perform other viewing opera tions without issuing commands. 
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The following picture shows the surface viewed by orbiting the camera toward 
the bottom using the Rotate 3D tool. A scene light has been added to illuminate 
the underside of the surface, which is not lit by the light added in the previous 
section. 




For More Information See "Lighting as a Visualization Tool" and "View 
Control with the Camera Toolbar" in the MATLAB documentation for 
information on these techniques. 
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Two-dimensional arrays can be displayed as images, where the array elements 
determine brightness or color of the images. For example, the statements 

load durer 

whos 

Ñame Size Bytes Class 

X 648x509 2638656 double array 

caption 2x28 112 char array 

map 128x3 3072 double array 

load the file durer . mat, adding three variables to the workspace. The matrix X 
is a 648-by-509 matrix and map is a 128-by-3 matrix that is the colormap for this 
image. 

MAT-files, such as durer .mat, are binary files that can be created on one 
platform and later read by MATLAB on a different platform. 

The elements of X are integers between 1 and 128, which serve as Índices into 
the colormap, map. Then 

image(X) 
colormap(map) 
axis image 

reproduces Dürer's etching shown at the beginning of this book. A 
high-resolution sean of the magic square in the upper right córner is available 
in another file. Type 

load detail 

and then use the up arrow key on your keyboard to reexecute the image, 
colormap, and axis commands. The statement 

colormap(hot) 

adds some twentieth century colorization to the sixteenth century etching. The 
function hot generates a colormap containing shades of reds, oranges, and 
yellows. Typically a given image matrix has a specific colormap associated with 
it. See the colormap reference page for a list of other predefined colormaps. 
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Readíng and Wrítíng Images 

You can read standard image files (TIFF, JPEG, BMP, etc.) into MATLAB 
using the imread function. The type of data returned by imread depends on the 
type of image you are reading. 

You can write MATLAB data to a variety of standard image formats using the 
imwrite function. See the reference pages for these functions for more 
information and examples. 



For More Information See "Displaying Bit-Mapped Images" in the 
MATLAB documentation for information on the image processing capabilities 
of MATLAB. 
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Printing Graphics 



You can print a MATLAB figure directly on a printer connected to your 
computer or you can export the figure to one of the standard graphics file 
formats supported by MATLAB. There are two ways to print and export 
figures: 

• Use the Print or Export Setup options under the File menú 

• Use the print command to print or export the figure 

See "Preparing Graphs for Presentation" on page 5-29 for an example. 

Printing from the Menú 

There are four menú options under the File menú that pertain to printing: 

• The Page Setup option displays a dialog box that enables you to adjust 
characteristics of the figure on the printed page. 

• The Print Setup option displays a dialog box that sets printing defaults, but 
does not actually print the figure. 

• The Print Preview option enables you to view the figure the way it will look 
on the printed page. 

• The Print option displays a dialog box that lets you select standard printing 
options and print the figure. 

Generally, use Print Preview to determine whether the printed output is what 
you want. If not, use the Page Setup dialog box to change the output settings. 
Select the Page Setup dialog box Help button to display information on how to 
set up the page. 

Exporting Figure to Graphics Files 

The Export Setup option in the File menú enables you to set a variety of figure 
characteristics, such as size and font type, as well as apply predefined 
templates to achieve standard-looking graphics files. After setup, you can 
export the figure to a number of standard graphics file formats. 
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Using the Print Command 

The print command provides more flexibility in the type of output sent to the 
printer and allows you to control printing from M-files. The result can be sent 
directly to your default printer or stored in a specified file. A wide variety of 
output formats, including TIFF, JPEG, and PostScript, is available. 

For example, this statement saves the contents of the current figure window as 
color Encapsulated Level 2 PostScript in the file called magicsquare . eps. It 
also includes a TIFF preview, which enables most word processors to display 
the picture. 

print -depsc2 -tiff magicsquare.eps 

To save the same figure as a TIFF file with a resolution of 200 dpi, use the 
command 

print -dtiff -r200 magicsquare. tiff 

If you type print on the command line, 
print 

MATLAB prints the current figure on your default printer. 



For More Information See the print reference page and "Printing and 
Exporting" in the MATLAB documentation for more information on printing. 
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Handle Graphics 



Handle Graphics refers to a system of graphics objects that MATLAB uses to 
implement graphing and visualization functions. Each object created has a 
fixed set of properties. You can use these properties to control the behavior and 
appearance of your graph. 

When you cali a plotting function, MATLAB creates the graph using various 
graphics objects, such as a figure window, axes, lines, text, and so on. MATLAB 
enables you to query the valué of each property and set the valúes of most 
properties. 

For example, the following statement creates a figure with a white background 
color and without displaying the figure toolbar. 

f igure( 'Color' , 'white ' , 'Toolbar' , ' none ' ) 

Using the Handle 

Whenever MATLAB creates a graphics object, it assigns an identifier (called a 
handle) to the object. You can use this handle to access the object's properties 
with the set and get functions. For example, the following statements créate 
a graph and return a handle to a lineseries object in h. 

x = 1:10; 
y = x.~3; 
h = plot(x,y) ; 

You can use the handle h to set the properties of the lineseries object. For 
example, you can set its Color property: 

set (h, 'Color' , ' red ' ) 

You can also specify properties when you cali the plotting function. 
h = plot(x,y, 'Color ' , ' red' ) ; 

When you query the lineseries properties, 

get(h, 'LineWidth' ) 

MATLAB returns the answer: 

ans = 

0.5000 
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Use the handle to see what properties a particular object contains. 
qet(h) 



Graphics Objects 



Graphics objects are the basic elements used to display graphs and user 
interface components. These objects are organized into a hierarchy, as shown 
by the following diagram. 



UI Objects 



Root 



Figure 



Axes 



Hidden 
Annotation Axes 



Core Objects Plot Objects Group Objects Annotation Objects 



Key Graphics Objects 

When you cali a function to créate a graph, MATLAB creates a hierarchy of 
graphics objects. For example, calling the plot function creates the following 
graphics objects: 

• Lineseries plot objects — Represent the data passed to the plot function 

• Axes — Provide a frame of reference and scaling for the plotted lineseries 

• Text — Label the axes tick marks and are used for titles and annotations 

• Figures — Are the windows that contain axes toolbars, menus, etc. 

Different types of graphs use different objects to represent data: however, all 
data objects are contained in axes and all objects (except root) are contained in 
figures. 

The root is an abstract object that primarily stores information about your 
computer or MATLAB state. You cannot créate an instance of the root object. 
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For More Information See "Handle Graphics Objects" in the MATLAB 
documentation for information on graphics objects. 



User interface objects are used to créate graphical user interfaces (GUIs). 
These objects include components like push buttons, editable text boxes, and 
list boxes. 



For More Information See "Creating Graphical User Interfaces" in the 
MATLAB documentation for more information on user interface objects. 



Creating Objects 

Plotting functions (like plot and surf ) cali the appropriate low-level function 
to draw their respective graph. For information about an object's properties, 
see the Handle Graphics Property Browser in the MATLAB online 
documentation. 

Commands for Working with Objects 

This table lists commands commonly used when working with objects. 

Function Purpose 

allchild Find all children of specified objects 

ancestor Find ancestor of graphics object 

copyobj Copy graphics object 

de le te Delete an object 



f indall Find all graphics objects (including hidden handles) 

f indobj Find the handles of objects having specified property 

valúes 



gca Return the handle of the current axes 
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Function Purpose 



gcf Return the handle of the current figure 

gco Return the handle of the current object 

get Query the valúes of an object's properties 

ishandle True if valué is valid object handle 



set Set the valúes of an object's properties 

Settíng Object Properties 

All object properties have default valúes. However, you might find it useful to 
change the settings of some properties to customize your graph. There are two 
ways to set object properties: 

• Specify valúes for properties when you créate the object. 

• Set the property valué on an object that already exists. 

Setting Properties from Plotting Commands 

You can specify object property valué pairs as arguments to many plotting 
functions, such as plot, mesh, and surf . 

For example, plotting commands that créate lineseries or surfaceplot objects 
enable you to specify property name/property valué pairs as arguments. The 
command 

surf (x,y,z, ' FaceColor ' , ' interp' , . . . 
' FaceLighting ' , ' gouraud ' ) 

plots the data in the variables x, y, and z using a surfaceplot object with 
interpolated face color and employing the Gouraud face light technique. You 
can set any of the object's properties this way. 

Setting Properties of Existing Objects 

To modify the property valúes of existing objects, you can use the set command 
or the Property Editor. This section describes how to use the set command. See 
"Using the Property Editor" on page 5-17 for more information. 
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Most plotting functions return the handles of the objects that they créate so you 
can modify the objects using the set command. For example, these statements 
plot a five-by-five matrix (creating five lineseries, one per column) and then set 
the Marker property to a square and the MarkerFaceColor property to green. 

h = plot (magic(5) ) ; 

set (h, 'Marker' , ' s ' , 'MarkerFaceColor' , 'g ' ) 

In this case, h is a vector containing five handles, one for each of the five 
lineseries in the graph. The set statement sets the Marker and 
MarkerFaceColor properties of all lineseries to the same valúes. 

Setting Múltiple Property Valúes 

If you want to set the properties of each lineseries to a different valué, you can 
use cell arrays to store all the data and pass it to the set command. For 
example, créate a plot and save the lineseries handles. 

h = plot (magic(5) ) ; 

Suppose you want to add different markers to each lineseries and color the 
marker's face color the same color as the lineseries. You need to define two cell 
arrays — one containing the property ñames and the other containing the 
desired valúes of the properties. 

The prop_name cell array contains two elements. 

prop_name(1) = {'Marker'}; 
prop_name(2) = {'MarkerFaceColor'}; 

The prop_values cell array contains 10 valúes: five valúes for the Marker 
property and five valúes for the MarkerFaceColor property. Notice that 
prop_values is a two-dimensional cell array. The first dimensión indicates 
which handle in h the valúes apply to and the second dimensión indicates 
which property the valué is assigned to. 

prop_values(1 , 1 ) = {'s'}; 

prop_values(1 ,2) = {get (h(1 ), 'Color ')} ; 

prop_values(2, 1 ) = { ' d ' } ; 

prop_values(2,2) = {get (h(2) ,' Color' )} ; 

prop_values(3, 1 ) = {'o'}; 

prop_values(3,2) = {get (h(3) ,' Color' )} ; 

prop_values(4, 1 ) = {'p'}; 

prop_values(4,2) = {get (h(4) ,' Color' )} ; 
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prop_values(5,1 ) 
prop_values(5,2) 



{'h'}j 

{get(h(5), 'Color')}; 



The MarkerFaceColor is always assigned the valué of the corresponding line's 
color (obtained by getting the lineseries Color property with the get 
command). 

After defining the cell arrays, cali set to specify the new property valúes. 

set (h,prop_name,prop_values) 
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Specifying the Axes or Figure 

MATLAB always creates an axes or figure if one does not exist when you issue 
a plotting command. However, when you are creating a graphics M-file, it is 
good practice to explicitly créate and specify the parent axes and figure, 
particularly if others will use your program. Specifying the parent prevenís the 
following problems: 

• Your M-file overwrites the graph in the current figure. Note that a figure 
becomes the current figure whenever a user clicks on it. 

• The current figure might be in an unexpected state and not behave as your 
program expects. 

The following examples shows a simple M-file that plots a function and the 
mean of the function over the specified range. 

function myfunc(x) 

% x = -10:. 005:40; Here's a valué you can use for x 

y = [1.5*cos(x) + 6*exp(-.1*x) + exp( .07*x) . *sin(3*x) ] ; 

ym = mean(y) ; 

hfig = figure( ' Ñame ',' Function and Mean',... 

'Pointer' , 'fullcrosshair' ) ; 
hax = axes( ' Parent ', hfig) ; 
plot (hax,x,y) 
hold on 

plot (hax, [min(x) max(x)],[ym ym] , 'Color ',' red ' ) 
hold off 

ylab = get (hax, ' YTick ' ) ; 
set (hax, ' YTick ' ,sort ( [ylab ym])) 

title ('y = 1.5cos(x) + 6e"{-0.1x} + e"{0 .07x}sin(3x) ' ) 
xlabel('X Axis'); ylabel('Y Axis') 
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y = 1 .5cos(x) + 6e °' 1x + e° 07x sin(3x) 




Finding the Handles of Existíng Objects 

The f indobj function enables you to obtain the handles of graphics objects by 
searching for objects with particular property valúes. With f indobj you can 
specify the valúes of any combination of properties, which makes it easy to pick 
one object out of many. f indobj also recognizes regular expressions (regexp). 

For example, you might want to find the blue line with square marker having 
blue face color. You can also specify which figures or axes to search, if there are 
more than one. The following sections provide examples illustrating how to use 
f indobj. 

Finding All Objects of a Certain Type 

Because all objects have a Type property that identifies the type of object, you 
can find the handles of all occurrences of a particular type of object. For 
example, 

h = f indobj ( 'Type ',' patch ') ; 
finds the handles of all patch objects. 
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Finding Objects with a Particular Property 

You can specify múltiple properties to narrow the search. For example, 

h = f indobj ( 'Type ',' line ', 'Color ',' r' ,' LineStyle ',':') ; 
finds the handles of all red dotted lines. 

Limiting the Scope of the Search 

You can specify the starting point in the object hierarchy by passing the handle 
of the starting figure or axes as the first argument. For example, 

h = f indobj (gca, 'Type ' , ' text ' , ' String ' , ' \pi/2 ' ) ; 
finds the string k/2 only within the current axes. 

Using findobj as an Argument 

Because findobj returns the handles it finds, you can use it in place of the 
handle argument. For example, 

set(findobj ( 'Type' , 'line' , 'Color' , ' red '), 'LineStyle',':') 
finds all red lines and sets their line style to dotted. 
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Animations 



MATLAB provides three ways of genera ting moving, animated graphics: 

• "Erase Mode Method" on page 5-71 — Continually erase and then redraw the 
objects on the screen, making incremental changes with each redraw. 

• "Creating Movies" on page 5-73 — Save a number of different pictures and 
then play them back as a movie. 

• Using AVI files. See avif ile for more information and examples. 

Erase Mode Method 

Using the EraseMode property is appropriate for long sequences of simple plots 
where the change from frame to frame is minimal. Here is an example showing 
simulated Brownian motion. Specify a number of points, such as 

n = 20 

and a temperature or velocity, such as 
s = .02 

The best valúes for these two parameters depend upon the speed of your 
particular computer. Genérate n random points with (x,y) coordinates between 
-V2 and +V2. 

x = rand(n,1 ) -0.5; 
y = rand(n,1 ) -0.5; 

Plot the points in a square with sides at -1 and +1. Save the handle for the 
vector of points and set its EraseMode to xor. This tells the MATLAB graphics 
system not to redraw the entire plot when the coordinates of one point are 
changed, but to restore the background color in the vicinity of the point using 
an exclusive or operation. 

h = plot(x,y, ' . '); 

axis([-1 1 -1 1]) 

axis square 

grid off 

set(h, 'EraseMode' , 'xor' , ' MarkerSize ' ,18) 
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Now begin the animation. Here is an infinite while loop, which you can 
eventually exit by typing Ctrl+C. Each time through the loop, add a small 
amount of normally distributed random noise to the coordinates of the points. 
Then, instead of creating an entirely new plot, simply change the XData and 
YData properties of the original plot. 

while 1 

drawnow 

x = x + s*randn(n, 1 ) ; 

y = y + s*randn(n, 1 ) ; 

set(h, 'XData' ,x, 'YData' ,y) 
end 

See how long it takes for one of the points to get outside the square and how 
long before all the points are outside the square. 
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Creatíng Movíes 

If you increase the number of points in the Brownian motion example to n = 
300 and s = . 02, the motion is no longer very fluid; it takes too much time to 
draw each time step. It becomes more effective to save a predetermined number 
of frames as bitmaps and to play them back as a movie. 

First, decide on the number of frames, 

nframes = 50; 

Next, set up the first plot as before, except using the default EraseMode 
(normal). 

x = rand(n,1 ) -0.5; 
y = rand(n,1 ) -0.5; 
h = plot(x,y, ' . ' 
set(h, 'MarkerSize' ,18) ; 
axis([-1 1 -1 1]) 
axis square 
grid off 

Genérate the movie and use getf rame to capture each frame. 

f or k = 1 ¡nframes 

x = x + s*randn(n, 1 ) ; 

y = y + s*randn(n, 1 ) ; 

set(h, 'XData' ,x, 'YData' ,y) 

M(k) = getframe; 
end 

Finally, play the movie 30 times, 
movie (M, 30) 
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This chapter introduces programming constructs, multidimensional arrays and other data 
structures, and MATLAB functions. 



Flow Control (p. 6-2) 

Other Data Structures (p. 6-9) 

Scripts and Functions (p. 6-19) 



Use flow control constructs including if , switch and 
case, for, while, continué, and break. 

Work with multidimensional arrays, cell arrays, 
character and text data, and structures. 

Write scripts and functions, use global variables, pass 
string arguments to functions, use eval to evalúate text 
expressions, vectorize code, preallocate arrays, reference 
functions using handles, and use functions that opérate 
on functions. 
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Flow Control 

MATLAB has several flow control constructs: 

• "if' on page 6-2 

• "switch and case" on page 6-4 

• "for" on page 6-5 

• "while" on page 6-6 

• "continué" on page 6-6 

• "break" on page 6-7 

• "try - catch" on page 6-8 

• "return" on page 6-8 

¡f 

The if statement evaluates a logical expression and executes a group of 
statements when the expression is true. The optional elseif and else 
keywords provide for the execution of altérnate groups of statements. An end 
keyword, which matches the if , terminates the last group of statements. The 
groups of statements are delineated by the four keywords — no braces or 
brackets are involved. 

The MATLAB algorithm for generating a magic square of order n involves 
three different cases: when n is odd, when n is even but not divisible by 4, or 
when n is divisible by 4. This is described by 

if rem(n,2) -= 

M = odd_magic(n) 
elseif rem(n,4) -= 

M = single_even_magic(n) 
else 

M = double_even_magic(n) 
end 

In this example, the three cases are mutually exclusive, but if they weren't, the 
first true condition would be executed. 
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It is important to understand how relational opera tors and if statements work 
with matrices. When you want to check for equality between two variables, you 
might use 

if A == B, ... 

This is valid MATLAB code, and does what you expect when A and B are 
scalars. But when A and B are matrices, A == B does not test ¿fthey are equal, 
it tests where they are equal; the result is another matrix of O's and l's showing 
element-by-element equality. (In fact, if A and B are not the same size, then A 
== B is an error.) 

B(1,1) = 0; 



A = magic 


(4); 


B = 


= A; 


A == B 








ans = 











1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



The proper way to check for equality between two variables is to use the 
isequal function, 

if isequal(A,B) , ... 

isequal returns a scalar logical valué of 1 (representing true) or (f alse), 
instead of a matrix, as the expression to be evaluated by the if function. Using 
the A and B matrices from above, you get 

isequal(A, B) 
ans = 



Here is another example to emphasize this point. If A and B are scalars, the 
following program will never reach the "unexpected situation". But for most 
pairs of matrices, including our magic squares with interchanged columns, 
none of the matrix conditions A > B, A < B, or A == B is true for all elements 
and so the else clause is executed. 
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if A > B 

' greater ' 
elseif A < B 

' less ' 
elseif A == B 

' equal ' 
else 

error ( ' Unexpected situation ' ) 
end 

Several functions are helpful for reducing the results of matrix comparisons to 
scalar conditions for use with if , including 

isequal 
isempty 
all 
any 

switch and case 

The switch statement executes groups of statements based on the valué of a 
variable or exprés sion. The keywords case and otherwise delinéate the 
groups. Only the first matching case is executed. There must always be an end 
to match the switch. 

The logic of the magic squares algorithm can also be described by 

switch (rem(n,4)==0) + (rem(n,2)==0) 
case 

M = odd_magic(n) 
case 1 

M = single_even_magic(n) 
case 2 

M = double_even_magic(n) 
otherwise 

error('This is impossible ' ) 
end 
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Note Unlike the C language switch statement, MATLAB switch does not 
fall through. If the first case statement is true, the other case statements do 
not execute. So, break statements are not required. 



for 

The for loop repeats a group of statements a fixed, predetermined number of 
times. A matching end delineates the statements. 

for n = 3:32 

r(n) = rank(magic(n) ) ; 
end 
r 

The semicolon terminating the inner statement suppresses repeated printing, 
and the r after the loop displays the final result. 



6-5 



6 Prog 



rammmg 



It is a good idea to indent the loops for readability, especially when they are 
nested. 

for i = 1 :m 

for ] = 1 :n 

H(i,j) = 1/(i+j); 

end 
end 

while 

The while loop repeats a group of statements an indefinite number of times 
under control of a logical condition. A matching end delineates the statements. 

Here is a complete program, illustrating while, if , else, and end, that uses 
interval bisection to find a zero of a polynomial. 

a = 0; f a = -Inf ; 
b = 3; fb = Inf; 
while b-a > eps*b 

x = (a+b)/2; 

fx = x"3-2*x-5; 

if sign(fx) == sign(fa) 
a = x; fa = fx; 

else 

b = x; fb = fx; 

end 
end 
x 

The result is a root of the polynomial x 3 - 2x - 5, namely 

x = 

2.09455148154233 

The cautions involving matrix comparisons that are discussed in the section on 
the if statement also apply to the while statement. 

continué 

The continué statement passes control to the next iteration of the for loop or 
while loop in which it appears, skipping any remaining statements in the body 
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of the loop. In nested loops, continué passes control to the next iteration of the 
f or loop or while loop enclosing it. 

The example below shows a continué loop that counts the lines of code in the 
file magic . m, skipping all blank lines and comments. A continué statement is 
used to advance to the next line in magic .m without incrementing the count 
whenever a blank line or comment line is encountered. 

fid = fopen( 'magic. m' ,' r' ) ; 
count = 0; 
while -feof(fid) 

line = fgetl(f id) ; 

if isempty (line) | strncmp(line, '%' ,1 ) 
continué 

end 

count = count + 1 ; 
end 
disp(sprintf ( '%d lines ', count) ) ; 

break 

The break statement lets you exit early from a for loop or while loop. In nested 
loops, break exits from the innermost loop only. 

Here is an improvement on the example from the previous section. Why is this 
use of break a good idea? 

a = 0; f a = -Inf ; 

b = 3; fb = Inf; 

while b-a > eps*b 

x = (a+b)/2; 

fx = x"3-2*x-5; 

if fx == 

break 
elseif sign(fx) == sign(fa) 

a = x; fa = fx; 
else 

b = x; fb = fx; 
end 
end 
x 
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try - catch 

The general form of a try-catch statement sequence is 

try 

statement 

statement 
catch 

statement 

statement 
end 

In this sequence the statements between try and catch are executed until an 
error occurs. The statements between catch and end are then executed. Use 
lasterr to see the cause of the error. If an error occurs between catch and end, 
MATLAB terminates execution unless another try-catch sequence has been 
established. 

return 

return terminates the current sequence of commands and returns control to 
the invoking function or to the keyboard. return is also used to terminate 
keyboard mode. A called function normally transfers control to the function 
that invoked it when it reaches the end of the function. You can insert a return 
statement within the called function to forcé an early termination and to 
transfer control to the invoking function. 
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This section introduces you to some other data structures in MATLAB, 
including 

• "Multidimensional Arrays" on page 6-9 

• "Cell Arrays" on page 6-11 

• "Characters and Text" on page 6-13 

• "Structures" on page 6-16 

Multidimensional Arrays 

Multidimensional arrays in MATLAB are arrays with more than two 
subscripts. One way of creating a multidimensional array is by calling zeros, 
ones, rand, or randn with more than two arguments. For example, 

R = randn(3,4,5) ; 

creates a 3-by-4-by-5 array with a total of 3x4x5 = 60 normally distributed 
random elements. 

A three-dimensional array might represent three-dimensional physical data, 
say the temperature in a room, sampled on a rectangular grid. Or it might 
represent a sequence of matrices, A*, or samples of a time-dependent matrix, 
A(t). In these latter cases, the (i, j)th element of the Mh matrix, or the í k th 
matrix, is denoted by A ( i , j , k ) . 

MATLAB and Dürer's versions of the magic square of order 4 differ by an 
interchange of two columns. Many different magic squares can be generated by 
interchanging columns. The statement 

p = perms(1 :4) ; 

generates the 4! = 24 permutations of 1 :4. The kth permutation is the row 
vector p ( k , : ) . Then 

A = magic(4) ; 

M = zeros(4,4,24) ; 

for k = 1 :24 

M(:,:,k) = A(:,p(k,:)); 
end 
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stores trie sequence of 24 magic squares in a three-dimensional array, M. The 
size of M is 

size(M) 



ans 



24 



13 


16 


2 


- ' 3 


8 


5 


U' ' 


10 


12 


9 ' ' 
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1. ' 
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14 


15 


' 
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16 
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. - 13 ' 




16 
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13 
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8 
12 


16 


2 


3 


13 


10 


5 
9 
4 


11 1 

7 
14 1 



6 
5 


8 

12 

1 


6 
15 


1 


.-'' 


,-' 












- 





Note The order of the matrices shown in this illustration might differ from 
your results. The perms function always returns all permutations of the input 
vector, but the order of the permutations might be different for different 
MATLAB versions. 



The statement 

sum(M,d) 

computes sums by varying the dth subscript. So 
sum(M, 1 ) 
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is a l-by-4-by-24 array containing 24 copies of the row vector 
34 34 34 34 



34 34 34 34 

and 

sum(M,2 



is a 4-by-l-by-24 array containing 24 copies of the column vector 

-5/1 



34 
34 
34 
34 



Finally, 

S = sum(M,3) 

adds the 24 matrices in the sequence. The result has size 4-by-4-by-l, so it looks 
like a 4-by-4 array. 



S = 

204 204 204 204 

204 204 204 204 

204 204 204 204 

204 204 204 204 



Cell Arrays 

Cell arrays in MATLAB are multidimensional arrays whose elements are 
copies of other arrays. A cell array of empty matrices can be created with the 
cell function. But, more often, cell arrays are created by enclosing a 
miscellaneous collection of things in curly braces, { }. The curly braces are also 
used with subscripts to access the contents of various cells. For example, 

C = {A sum(A) prod(prod(A) ) } 

produces a l-by-3 cell array. The three cells contain the magic square, the row 
vector of column sums, and the product of all its elements. When C is displayed, 
you see 

C = 

[4x4 double] [1x4 double] [20922789888000] 
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This is because the first two cells are too large to print in this limited space, but 
the third cell contains only a single number, 16!, so there is room to print it. 

Here are two important points to remember. First, to retrieve the contents of 
one of the cells, use subscripts in curly braces. For example, C{1 } retrieves the 
magic square and C{3} is 16. Second, cell arrays contain copies of other arrays, 
not pointers to those arrays. If you subsequently change A, nothing happens to 
C. 

You can use three-dimensional arrays to store a sequence of matrices of the 
same size. Cell arrays can be used to store a sequence of matrices oidifferent 
sizes. For example, 

M = cell(8,1) ; 
for n = 1 :8 

M{n} = magic(n) ; 
end 

M 

produces a sequence of magic squares of different order. 
M = 





1 


[ 2x2 


double 


[ 3x3 


double 


[ 4x4 


double 


[ 5x5 


double 


[ 6x6 


double 


[ 7x7 


double 


[ 8x8 


double 
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You can retrieve the 4-by-4 magic square matrix with 

M{4} 

Characters and Text 

Enter text into MATLAB using single quotes. For example, 
s = 'Helio' 

The result is not the same kind of numeric matrix or array you have been 
dealing with up to now. It is a l-by-5 character array. 
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Internally, the characters are stored as numbers, but not in floating-point 
format. The statement 



a = double(s) 

converts the character array to a numeric matrix containing floating-point 
representations of the ASCII codes for each character. The result is 



72 101 108 108 111 

The statement 
s = char(a) 

reverses the conversión. 

Converting numbers to characters makes it possible to investígate the various 
fonts available on your computer. The printable characters in the basic ASCII 
character set are represented by the integers 32 : 1 27. (The integers less than 
32 represent nonprintable control characters.) These integers are arranged in 
an appropriate 6-by-16 array with 

F = reshape(32:127,16,6) ' ; 

The printable characters in the extended ASCII character set are represented 
by F+128. When these integers are interpreted as characters, the result 
depends on the font currently being used. Type the statements 

char(F) 
char(F+128) 

and then vary the font being used for the Command Window. Select 
Preferences from the File menú to change the font. If you include tabs in lines 
of code, use a fixed-width font, such as Monospaced, to align the tab positions 
on different lines. 

Concatenation with square brackets joins text variables together into larger 
strings. The statement 

h = [s, ' world ' ] 

joins the strings horizontally and produces 

h = 

Helio world 
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The statement 

v = [s; 'world ' ] 

joins the strings vertically and produces 

v = 

Helio 
world 

Note that a blank has to be inserted before the ' w ' in h and that both words in 
v have to have the same length. The resulting arrays are both character arrays; 
h is 1-by-ll and v is 2-by-5. 

To manipúlate a body of text containing lines of different lengths, you have two 
choices — a padded character array or a cell array of strings. When creating a 
character array, you must make each row of the array the same length. (Pad 
the ends of the shorter rows with spaces.) The char function does this padding 
for you. For example, 

S = char( ' A' , ' rolling ' , ' stone ' , ' gathers ' , 'momentum. ' ) 

produces a 5-by-9 character array. 

S = 

A 

rolling 

stone 

gathers 

momentum. 

Alternatively, you can store the text in a cell array. For example, 

C = { 'A' ;' rolling ';' stone '; 'gathers ';' momentum. ' } 

creates a 5-by-l cell array that requires no padding because each row of the 
array can have a different length. 

C = 

'A' 

' rolling ' 
' stone ' 
'gathers ' 
'momentum. ' 
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You can convert a padded character array to a cell array of strings with 
C = cellstr(S) 

and reverse the process with 
S = char(C) 

Structures 

Structures are multidimensional MATLAB arrays with elements accessed by 
textual field designators. For example, 

S.name = ' Ed Plum ' ; 
S.score = 83; 
S. grade = 'B+' 

creates a scalar structure with three fields. 

S = 

ñame: ' Ed Plum' 

score: 83 

grade: 'B+' 

Like everything else in MATLAB, structures are arrays, so you can insert 
additional elements. In this case, each element of the array is a structure with 
several fields. The fields can be added one at a time, 

S(2).name = 'Toni Miller'; 
S(2) .score = 91 ; 
S(2) .grade = 'A- ' ; 

or an entire element can be added with a single statement. 

S(3) = struct ( ' ñame ' , ' Jerry García',... 
' score ' ,70, 'grade ' , ' C ) 

Now the structure is large enough that only a summary is printed. 

S = 

1x3 struct array with fields: 

ñame 

score 

grade 



6-16 



Other Data Structur 



There are several ways to reassemble the various fields into other MATLAB 
arrays. They are all based on the notation of a comma-separated list. If you type 

S.score 

it is the same as typing 

S(1).score, S(2).score, S(3).score 

This is a comma-separated list. Without any other punctuation, it is not very 
useful. It assigns the three scores, one at a time, to the default variable ans and 
dutifully prints out the result of each assignment. But when you endose the 
expression in square brackets, 

[S. score] 

it is the same as 

[S(1). score, S(2). score, S(3) . score] 

which produces a numeric row vector containing all the scores. 

ans = 

83 91 70 

Similarly, typing 
S.name 

just assigns the ñames, one at a time, to ans. But enclosing the expression in 
curly braces, 

{S.name} 
creates a l-by-3 cell array containing the three ñames. 



ans = 






'Ed Plum 


'Toni Miller' 


'Jerry Garcia 


And 






char (S.name) 
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calis the char function with three arguments to créate a character array from 
the ñame fields, 

ans = 
Ed Plum 
Toni Miller 
Jerry García 

Dynamic Field Ñames 

The most common way to access the data in a structure is by specifying the 
ñame of the field that you want to reference. Another means of accessing 
structure data is to use dynamic field ñames. These ñames express the field as 
a variable expression that MATLAB evaluates at run-time. The 
dot-parentheses syntax shown here makes expression a dynamic field ñame. 

structName. (expression) 

Index into this field using the standard MATLAB indexing syntax. For 
example, to evalúate expression into a field ñame and obtain the valúes of that 
field at columns 1 through 25 of row 7, use 

structName . (expression) (7,1 :25) 

Dynamic Field Ñames Example. The avgscore function shown below computes an 
average test score, retrieving information from the testscores structure using 
dynamic field ñames: 

function avg = avgscore(testscores, student, first, last) 
for k = first:last 

scores(k) = testscores. (student) .week(k) ; 
end 
avg = sum(scores) / (last - first +1); 

You can run this function using different valúes for the dynamic field student. 

avgscore(testscores, 'Ann Lañe', 1, 20) 
ans = 

83.5000 

avgscore(testscores, 'William King ' , 1, 20) 
ans = 

92.1000 



6-18 



Scripts and Functions 



Scripts and Functions 

Topics covered in this section are 



• "Scripts" on page 6-20 

• "Functions" on page 6-21 

• "Global Variables" on page 6-25 

• "Passing String Arguments to Functions" on page 6-26 

• "The eval Function" on page 6-27 

• "Function Handles" on page 6-28 

• "Function Functions" on page 6-28 

• "Vectorization" on page 6-31 

• "Preallocation" on page 6-31 

MATLAB is a powerful programming language as well as an interactive 
computational environment. Files that contain code in the MATLAB language 
are called M-files. You créate M-files using a text editor, then use them as you 
would any other MATLAB function or command. 

There are two kinds of M-files: 

• Scripts, which do not accept input arguments or return output arguments. 
They opérate on data in the workspace. 

• Functions, which can accept input arguments and return output arguments. 
Internal variables are local to the function. 

If you're a new MATLAB programmer, just créate the M-files that you want to 
try out in the current directory. As you develop more of your own M-files, you 
will want to organize them into other directories and personal toolboxes that 
you can add to your MATLAB search path. 

If you duplicate function ñames, MATLAB executes the one that occurs first in 
the search path. 

To view the contents of an M-file, for example, myf unction . m, use 

type myfunction 
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Scripts 

When you invoke a script, MATLAB simply executes the commands found in 
the file. Scripts can opérate on existing data in the workspace, or they can 
créate new data on which to opérate. Although scripts do not return output 
arguments, any variables that they créate remain in the workspace, to be used 
in subsequent computations. In addition, scripts can produce graphical output 
using functions like plot. 

For example, créate a file called magicrank . m that contains these MATLAB 
commands. 

% Investígate the rank of magic squares 
r = zeros(1 ,32) ; 
for n = 3:32 

r(n) = rank(magic(n) ) ; 
end 
r 
bar(r) 

Typing the statement 
magicrank 

causes MATLAB to execute the commands, compute the rank of the first 30 
magic squares, and plot a bar graph of the result. After execution of the file is 
complete, the variables n and r remain in the workspace. 
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Functions 

Functions are M-files that can accept input arguments and return output 
arguments. The ñames of the M-file and of the function should be the same. 
Functions opérate on variables within their own workspace, sepárate from the 
workspace you access at the MATLAB command prompt. 

A good example is provided by rank. The M-file rank . m is available in the 
directory 

toolbox/matlab/matfun 

You can see the file with 
type rank 
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Here is the file. 

function r = rank(A,tol) 

% RANK Matrix rank. 

% RANK(A) provides an estimate of the number of linearly 

% independent rows or columns of a matrix A. 

% RANK(A,tol) is the number of singular valúes of A 

% that are larger than tol. 

% RANK(A) uses the default tol = max(size(A) ) * norm(A) * eps. 

s = svd(A) ; 
if nargin==1 

tol = max(size(A) ' ) * max(s) * eps; 
end 
r = sum(s > tol) ; 

The first line of a function M-file starts with the keyword function. It gives the 
function ñame and order of arguments. In this case, there are up to two input 
arguments and one output argument. 

The next several lines, up to the first blank or executable line, are comment 
lines that provide the help text. These lines are printed when you type 

help rank 

The first line of the help text is the Hl line, which MATLAB displays when you 
use the lookf or command or request help on a directory. 

The rest of the file is the executable MATLAB code defining the function. The 
variable s introduced in the body of the function, as well as the variables on the 
first line, r, A and tol, are all local to the function; they are sepárate from any 
variables in the MATLAB workspace. 

This example illustrates one aspect of MATLAB functions that is not ordinarily 
found in other programming languages — a variable number of arguments. 
The rank function can be used in several different ways. 

rank(A) 

r = rank(A) 

r = rank(A,1 .e-6) 
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Many M-files work this way. If no output argument is supplied, the result is 
stored in ans. If the second input argument is not supplied, the function 
computes a default valué. Within the body of the function, two quantities 
named nargin and nargout are available which tell you the number of input 
and output arguments involved in each particular use of the function. The rank 
function uses nargin, but does not need to use nargout. 

Types of Functions 

MATLAB offers several different types of functions to use in your 
programming. 

Anonymous Functions 

An anonymous function is a simple form of MATLAB function that does not 
require an M-file. It consists of a single MATLAB expression and any number 
of input and output arguments. You can define an anonymous function right at 
the MATLAB command line, or within an M-file function or script. This gives 
you a quick means of creating simple functions without having to créate M-files 
each time. 

The syntax for creating an anonymous function from an expression is 

f = @(arglist)expression 

The statement below creates an anonymous function that finds the square of a 
number. When you cali this function, MATLAB assigns the valué you pass in 
to variable x, and then uses x in the equation x . "2. 

sqr = @(x) x. "2; 

To execute the sqr function defined above, type 

a = sqr(5) 
a = 
25 

Primary and Subfunctions 

All functions that are not anonymous must be defined within an M-file. Each 
M-file has a required primary function that appears first in the file, and any 
number oí subfunctions that follow the primary. Primary functions have a 
wider scope than subfunctions. That is, primary functions can be invoked from 
outside of their M-file (from the MATLAB command line or from functions in 
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other M-files) while subfunctions cannot. Subfunctions are visible only to the 
primary function and other subfunctions within their own M-file. 

The rank function shown in the section on "Functions" on page 6-21 is an 
example of a primary function. 

Prívate Functions 

A prívate function is a type of primary M-file function. Its unique characteristic 
is that it is visible only to a limited group of other functions. This type of 
function can be useful if you want to limit access to a function, or when you 
choose not to expose the implementation of a function. 

Prívate functions reside in subdirectories with the special ñame prívate . They 
are visible only to functions in the parent directory. For example, assume the 
directory newmath is on the MATLAB search path. A subdirectory of newmath 
called prívate can contain functions that only the functions in newmath can 
cali. 

Because prívate functions are invisible outside the parent directory, they can 
use the same ñames as functions in other directories. This is useful if you want 
to créate your own versión of a particular function while retaining the original 
in another directory. Because MATLAB looks for prívate functions before 
standard M-file functions, it will find a prívate function named test .m before 
a nonprivate M-file named test . m. 

Nested Functions 

You can define functions within the body of any MATLAB M-file function. 
These are said to be nested within the outer function. A nested function 
contains any or all of the components of any other M-file function. In this 
example, function B is nested in function A. 

function x = A(p1 , p2) 

B(p2) 

function y = B(p3) 

end 
end 
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Like other functions, a nested function has its own workspace where variables 
used by the function are stored. But it also has access to the workspaces of all 
functions in which it is nested. So, for example, a variable that has a valué 
assigned to it by the primary function can be read or overwritten by a function 
nested at any level within the primary. Similarly, a variable that is assigned in 
a nested function can be read or overwritten by any of the functions containing 
that function. 

Function Overloading 

Overloaded functions act the same way as overloaded functions in most 
computer languages. Overloaded functions are useful when you need to créate 
a function that responds to different types of inputs accordingly. For instance, 
you might want one of your functions to accept both double-precision and 
integer input, but to handle each type somewhat differently. You can make this 
difference invisible to the user by creating two sepárate functions having the 
same ñame, and designating one to handle double types and one to handle 
integers. When you cali the function, MATLAB chooses which M-file to 
dispatch to based on the type of the input arguments. 

Global Variables 

If you want more than one function to share a single copy of a variable, simply 
declare the variable as global in all the functions. Do the same thing at the 
command line if you want the base workspace to access the variable. The global 
declaration must occur before the variable is actually used in a function. 
Although it is not required, using capital letters for the ñames of global 
variables helps distinguish them from other variables. For example, créate an 
M-file called f alling .m. 

function h = falling(t) 

global GRAVITY 

h = 1/2*GRAVITY*t.~2; 

Then interactively enter the statements 

global GRAVITY 

GRAVITY = 32; 

y = falling((0:.1:5)'); 
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The two global statements make the valué assigned to GRAVITY at the 
command prompt available inside the function. You can then modify GRAVITY 
interactively and obtain new solutions without editing any files. 

Passíng Stríng Arguments to Functíons 

You can write MATLAB functions that accept string arguments without the 
parentheses and quotes. That is, MATLAB interprets 

foo abe 

as 

f oo( 'a','b','c') 

However, when you use the unquoted form, MATLAB cannot return output 
arguments. For example, 

legend apples oranges 

creates a legend on a plot using the strings apples and oranges as labels. Ifyou 
want the legend command to return its output arguments, then you must use 
the quoted form. 

[legh,objh] = legend( ' apples ',' oranges ') ; 

In addition, you cannot use the unquoted form if any of the arguments is not a 
string. 

Constructing String Arguments in Code 

The quoted form enables you to construct string arguments within the code. 
The following example processes múltiple data files, Augustl .dat, 
August2 .dat, and so on. It uses the function int2str, which converts an 
integer to a character, to build the filename. 

for d = 1 :31 

s = ['August' int2str(d) '.dat']; 
load(s) 

% Code to process the contents of the d-th file 
end 
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A Cautionary Note 

While the unquoted syntax is convenient, it can be used incorrectly without 
causing MATLAB to genérate an error. For example, given a matrix A, 

A = 

-6 -1 

6 2 -16 

-5 20 -10 

The eig command returns the eigenvalues of A. 

eig(A) 

ans = 
-3.0710 

-2. 4645+17. 6008Í 
-2. 4645-17. 6008Í 

The following statement is not allowed because A is not a string; however, 
MATLAB does not genérate an error. 

eig A 

ans = 

65 

MATLAB actually takes the eigenvalue of the ASCII numeric equivalent of the 
letter A (which is the number 65). 

The eval Function 

The eval function works with text variables to implement a powerful text 
macro facility. The expression or statement 

eval(s) 

uses the MATLAB interpreter to evalúate the expression or execute the 
statement contained in the text string s. 

The example of the previous section could also be done with the following code, 
although this would be somewhat less efficient because it involves the full 
interpreter, not just a function cali. 



6-27 



6 Prog 



rammmg 



for d = 1 :31 

s = ['load August' int2str(d) '.dat']; 

eval(s) 

% Process the contents of the d-th file 
end 

Functíon Handles 

You can créate a handle to any MATLAB function and then use that handle as 
a means of referencing the function. A function handle is typically passed in an 
argument list to other functions, which can then execute, or evalúate, the 
function using the handle. 

Construct a function handle in MATLAB using the at sign, @, before the 
function ñame. The following example creates a function handle for the sin 
function and assigns it to the variable f handle. 

fhandle = ©sin; 

You can cali a function by means of its handle in the same way that you would 
cali the function using its ñame. The syntax is 

f handle(arg1 , arg2, ...); 

The function plot_f handle, shown below, receives a function handle and data, 
generates y-axis data using the function handle, and plots it. 

function x = plot_fhandle(f handle, data) 
plot(data, f handle(data) ) 

When you cali plot_f handle with a handle to the sin function and the 
argument shown below, the resulting evaluation produces a sine wave plot. 

plot_f handle(@sin, -pi:0.01:pi) 

Function Functions 

A class of functions called "function functions" works with nonlinear functions 
of a scalar variable. That is, one function works on another function. The 
function functions include 

• Zero finding 

• Optimization 
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• Quadrature 

• Ordinary differential equations 

MATLAB represents the nonlinear function by a function M-file. For example, 
here is a simplified versión of the function humps from the matlab/demos 
directory. 

function y = humps(x) 

y = 1 ./((x-.3) ."2 + .01) + 1 ./((x-.9) .'2 + .04) - 6; 

Evalúate this function at a set of points in the interval < x < 1 with 

x = 0: .002:1 ; 
y = humps(x) ; 

Then plot the function with 
plot(x,y) 
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The graph shows that the function has a local mínimum near x = 0.6. The 
function f minsearch finds the minimizer, the valué of x where the function 
takes on this minimum. The first argument to f minsearch is a function handle 
to the function being minimized and the second argument is a rough guess at 
the location of the minimum. 

p = fminsearch(@humps, .5) 

P = 

0.6370 

To evalúate the function at the minimizer, 
humps(p) 

ans = 

11 .2528 

Numerical analysts use the terms quadrature and integration to distinguish 
between numerical approximation of definite integráis and numerical 
integration of ordinary differential equations. MATLAB quadrature routines 
are quad and quadl. The statement 

Q = quadl (@humps, 0, 1 ) 

computes the área under the curve in the graph and produces 

Q = 

29.8583 

Finally, the graph shows that the function is never zero on this interval. So, if 
you search for a zero with 

z = f zero(@humps, .5) 

you will find one outside the interval 

z = 

-0.1316 
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Vectorízatíon 

One way to make your MATLAB programs run faster is to vectorize the 
algorithms you use in constructing the programs. Where other programming 
languages might use f or loops or DO loops, MATLAB can use vector or matrix 
operations. A simple example involves creating a table of logarithms. 

x = .01 ; 

for k = 1 : 1001 

y(k) = log10(x); 
x = x + .01 ; 
end 

A vectorized versión of the same code is 

x = .01 : .01 :10; 
y = log10(x) ; 

For more complicated code, vectorization options are not always so obvious. 



For More Information See "Improving Performance and Memory Usage" in 
the MATLAB Programming documentation for other techniques that you can 
use. 



Preallocatíon 

If you can't vectorize a piece of code, you can make your for loops go faster by 
preallocating any vectors or arrays in which output results are stored. For 
example, this code uses the function zeros to preallocate the vector created in 
the for loop. This makes the for loop execute significantly faster. 

r = zeros(32,1 ) ; 
fon n = 1 :32 

r(n) = rank(magic(n) ) ; 
end 

Without the preallocation in the previous example, the MATLAB interpreter 
enlarges the n vector by one element each time through the loop. Vector 
preallocation eliminates this step and results in faster execution. 
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This chapter introduces GUIDE, the MATLAB graphical user interface development environment, 
which is a tool that you can use to créate graphical user interfaces. 

What Is GUIDE? (p. 7-2) Introduces GUIDE, the MATLAB graphical user 

interface design environment 

Laying Out a GUI (p. 7-3) Briefly describes the GUIDE Layout Editor 

Programming a GUI (p. 7-6) Introduces callbacks to define behavior of the GUI 

components 



m Creating Graphical User Interfaces 



What \s GUIDE? 
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GUIDE, the MATLAB graphical user interface development environment, 
provides a set of tools for creating graphical user interfaces (GUIs). These tools 
greatly simplify the process of designing and building GUIs. You can use the 
GUIDE tools to 



• 



• 



Lay out the GUI. 

Using the GUIDE Layout Editor, you can lay out a GUI easily by clicking and 
dragging GUI components — such as panels, buttons, text fields, sliders, 
menus, and so on — into the layout área. GUIDE stores the GUI layout in a 
FIG-file. 

Program the GUI. 

GUIDE automatically generates an M-file that controls how the GUI 
opera tes. The M-file initializes the GUI and contains a framework for the 
most commonly used callbacks for each component — the commands that 
execute when a user clicks a GUI component. Using the M-file editor, you can 
add code to the callbacks to perform the functions you want. 



Laying Out a GUI 



Layíng Out a GUI 



Starting GUIDE 

Start GUIDE by typing guide at the MATLAB command prompt. This displays 
the GUIDE Quick Start dialog, as shown in the following figure. 



Créate New GUI | Open Exlstlncj GUI | 



GUIDE templates 



•#kBlankGUI(DBfault) 
•^kGUIwlth Ulcontrols 
■^ GUI wlth tees and Menu 
^k Modal Question Díalos 



BLANK 



V Isave orí startup asi 



D:\WorMGUIDEUjntltled.flg 



OK 



Cancel 



Help 



From the GUIDE Quick Start dialog, you can 

• Créate a new GUI from one of the GUIDE templates — prebuilt GUIs that 
you can modify for your own purposes. 

• Open an existing GUI. 
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The Layout Editor 

When you open a GUI in GUIDE, it is displayed in the Layout Editor, which is 
the control panel for all of the GUIDE tools. The following figure shows the 
Layout Editor with a blank GUI témplate. 
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You can lay out your GUI by dragging components, such as panels, push 
buttons, pop-up menus, or axes, from the component palette, at the left side of 
the Layout Editor, into the layout área. For example, if you drag a push button 
into the layout área, it appears as in the following figure. 
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You can also use the Layout Editor to créate menus and set basic properties of 
the GUI components. 

To get started using the Layout Editor and setting property valúes, see 
"Creating a GUI" in the MATLAB documentation. 
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Programmíng a GUI 



After laying out the GUI and setting component properties, the next step is to 
program the GUI. You program the GUI by coding one or more callbacks for 
each of its components. Callbacks are functions that execute in response to 
some action by the user. A typical action is clicking a push button. 

A GUI's callbacks are found in an M-file that GUIDE generates automatically. 
GUIDE adds templates for the most commonly used callbacks to this M-file, 
but you may want to add others. Use the M-file Editor to edit this file. 

The following figure shows the Callback témplate for a push button. 
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h Executes 011 button press in pushbuttonl. 

function pushbuttonl_Callback (hObject, eventdata, handles) 

% hObject handle to pushbuttonl (see GCBO) 

% eventdata reserved - to toe defined in a future versión of HílTLAB 

% handles stuucture with handles and user: data [see GUIDATA) 
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To learn more about programming a GUI, see "Creating a GUI" in the 
MATLAB documentation. 
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This chapter introduces the MATLAB development environment and the tools you can use to manage 
your work in MATLAB. 



Desktop Overview (p. 8-2) 

Command Window and Command 
History (p. 8-4) 

Help Browser (p. 8-6) 



Access tools, arrange the desktop, and set preferences. 
Run functions and enter variables. 

Find and view documentation and demos. 



Current Directory Browser and Search Manage and use M-files with MATLAB. 
Path (p. 8-9) 



Workspace Browser and Array Editor 
(p. 8-11) 

Editor/Debugger (p. 8-13) 

M-Lint Code Check and Profiler 
Reports (p. 8-15) 

Other Development Environment 
Features (p. 8-18) 



Work with variables in MATLAB. 

Créate and debug M-files (MATLAB programs). 
Improve and tune your M-files. 

Interface with source control systems, and publish M-file 
results. 
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Desktop Overview 



Use desktop tools to manage your work in MATLAB. You can also use 
MATLAB functions to perform the equivalent of most of the features found in 
the desktop tools. 

The following illustration shows the default configuration of the MATLAB 
desktop. You can modify the setup to meet your needs. 



Menus change, depending on the tool you 
are currently using. 
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To ¡jet started, select MATLAB Help or Demos 
» 



T3 




For More Information For an overview of the desktop tools, watch the video 
tutorials, accessible by typing demo matlab desktop. For complete details, see 
the MATLAB Desktop Tools and Development Environment documentation. 
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Arrangíng the Desktop 

These are some common ways to customize the desktop: 

• Show or hide desktop tools via the Desktop menú. 

• Resize any tool by dragging one of its edges. 

• Move a tool outside of the desktop by clicking the undock button ¿I in the 
tool's title bar. 

• Reposition a tool within the desktop by dragging its title bar to the new 
location. Tools can occupy the same position, as shown for the Current 
Directory and Workspace browser in the preceding illustration, in which 
case, you access a tool via its tab. 

• Change fonts and other options by using File -> Preferences. 



Start Button 

The MATLAB Start button provides easy access to tools, demos, shortcuts, and 
documentation. Click the Start button to see the options. 
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For More Information See "Desktop" in the MATLAB Desktop Tools and 
Development Environment documentation. 
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Command Window and Command Hístory 



Command Window 

Use the Command Window to enter variables and to run functions and M-file 
scripts. 
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Press the up arrow key 1" to recall a statement you previously typed. Edit the 
statement as needed and then press Enter to run it. For more information 
about entering statements in the Command Window, see "Controlling 
Command Window Input and Output" on page 4-28. 



For More Information See "Running Functions" in the MATLAB Desktop 
Tools and Development Environment documentation for complete details. 
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Command History 

Statements you enter in the Command Window are logged in the Command 
History. From the Command History, you can view previously run statements, 
as well as copy and execute selected statements. You can also créate an M-file 
from selected statements. 



Timestamp 
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To save the input and output from a MATLAB session to a file, use the diary 
function. 



For More Information See "Command History" in the MATLAB Desktop 
Tools and Development Environment documentation, and the reference page 
for the diary function. 
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Help Browser 



Use the Help browser to search for and view documentation and demos for all 
your MathWorks producís. The Help browser is an HTML viewer integrated 
into the MATLAB desktop. 

To open the Help browser, click the help button f" in the desktop toolbar. 

The Help browser consists of two panes, the Help Navigator, which you use to 
find information, and the display pane, where you view the information. These 
are the key features: 

• Contents tab — View the titles and tables of contents of the documentation. 

• Index tab — Find specific index entries (selected keywords) in the 
documentation. 

• Search tab — Look for specific words in the documentation. 

• Demos tab — View and run demonstrations for your MathWorks producís. 

While viewing the documentation, you can 

• Browse to other pages — Use the arrows at the tops and bottoms of the pages 
to move through the document, or use the back and forward buttons in the 
toolbar to go to previously viewed pages. 

• Bookmark pages — Use the Favorites menú. 

• Print a page — Click the print button in the toolbar. 

• Find a term in the page — Click the find icón (#4) in the toolbar. 

• Copy or evalúate a selection — Select text, such as code from an example, 
then right-click and use a context menú item to copy the selection or evalúate 
(run) it. 
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Tnbs ¡n tlie Help Navigator pane provide difieren! ways to find information. 
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Other Forms of Help 

In addition to the Help browser, you can use help functions. To get help for a 
specific function, use the doc function. For example, doc f ormat displays 
documentation for the f ormat function in the Help browser. 

To see a briefer form of the documentation for a function, type help followed by 
the function ñame. The resulting help text appears in the Command Window. 
It shows function ñames in all capital letters to distinguish them from the 
surrounding text. When you use the function ñames, type them in lowercase or 
they will not run. Some functions actually consist of both uppercase and 
lowercase letters, and the help text clearly indicates that. For those functions, 
match the case used in the help function. 

Other means for getting help include contacting Technical Support 

(http : / /www . mathworks . com/ support) and participating in the newsgroup for 

MATLAB users, comp . sof t - sys . matlab. 



For More Information See "Help for Using MATLAB" in the MATLAB 
Desktop Tools and Development Environment documentation, and the 
reference pages for the doc and help functions. 
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Current Directory Browser and Search Path 

MATLAB file operations use the current directory and the search path as 
reference points. Any file you want to run must either be in the current 
directory or on the search path. 

Current Directory 

A quick way to view or change the current directory is by using the current 
directory field in the desktop toolbar, shown here. 



| D:Vnymfiles 



"Z J a 



To search for, view, open, and make changes to MATLAB related directories 
and files, use the MATLAB Current Directory browser. Alternatively, you can 
use the functions dir, cd, and delete. Use the Visual Directory and Directory 
Reports to help you manage M-files. 
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For More Information See "File Management Operations" in the MATLAB 
Desktop Tools and Development Environment documentation, and the 
reference pages for the dir, cd, and delete functions. 



Sea re h Path 

MATLAB uses a search path to find M-files and other MATLAB related files, 
which are organized in directories on your file system. Any file you want to run 
in MATLAB must reside in the current directory or in a directory that is on the 
search path. When you créate M-files and related files for MATLAB, add the 
directories in which they are located to the MATLAB search path. By default, 
the files supplied with MATLAB and other MathWorks producís are included 
in the search path. 

To see which directories are on the search path or to change the search path, 
select File -> Set Path and use the resulting Set Path dialog box. 
Alternatively, you can use the path function to view the search path, addpath 
to add directories to the path, and rmpath to remove directories from the path. 



For More Information See "Search Path" in the MATLAB Desktop Tools 
and Development Environment documentation, and the reference pages for 
the path, addpath, and rmpath functions. 
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Workspace Browser and Array Editor 



The MATLAB workspace consists of the set of variables (named arrays) built 
up during a MATLAB session and stored in memory. You add variables to the 
workspace by using functions, running M-files, and loading saved workspaces. 



Workspace Browser 



To view the workspace and information about each variable, use the 
Workspace browser, or use the functions who and whos. 
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To delete variables from the workspace, select the variables and select Edit -> 
Delete. Alternatively, use the clear function. 

The workspace is not maintained after you end the MATLAB session. To save 
the workspace to a file that can be read during a later MATLAB session, select 
File -> Save, or use the save function. This saves the workspace to a binary file 
called a MAT-file, which has a .mat extensión. You can use options to save to 
different formats. To read in a MAT-file, select File -> Import Data, or use the 
load function. 



For More Information See "MATLAB Workspace" in the MATLAB Desktop 
Tools and Development Environment documentation, and the reference pages 
for the who, clear, save, and load functions. 
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Array Editor 

Double-click a variable in the Workspace browser, or use openvar 
variablename, to see it in the Array Editor. Use the Array Editor to view and 
edit a visual representation of variables in the workspace. 



View and change valúes 
of array elements. 



!=£! Array Editor - ans 



Arrange the display of 

array documenls. 



File Edit View Graphics Debug Desktor: Wíndow Help "* " x 




Use documenl bar to view other variables you bave open in the Array Editor. 



For More Information See "Viewing and Editing Workspace Variables with 
the Array Editor" in the MATLAB Desktop Tools and Development 
Environment documentation, and the reference page for the openvar function. 
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Edítor/Debugger 



Use the Editor/Debugger to créate and debug M-files, which are programs you 
write to run MATLAB functions. The Editor/Debugger provides a graphical 
user interface for text editing, as well as for M-file debugging. To créate or edit 
an M-file use File -> New or File -> Open, or use the edit function. 



Set breakpoints 

whereyouwant 
execution to 
pause so you 
can examii 
variables. 



Holdthe 
cursor over a, 
variable and 
its current 
valué appears 
(known as a 
data tip). 



Commenl selecled lines and specify 
indenting style using tbe Text menú. 



Find and replace text. 



Arrange the display of 
documents in the Editor. 



fl Editor - d:\mymfiles\collatz.m 



File Edit Text Cell Tools Debug Desktop Wlndow Help 



D té 



juen" 



/.ítem 



1 1 coll... J_\ \Ó jj 



-|qi| x| 

X 



function seguence=collatz (n) 

% Collatz pnonlem. Genérate a sequence of integens resolving to 1 

% For any positive integer, n: 

% Divide n by 2 i£ n 13 even 

% Hultiply n by 3 and add 1 if n is odd 



Repeat f^ — tuna — raenl t 

. n: lxl double 




t is lí 



sequence = ñT; 
next vaioe = n; 
Sile next_value > 1 

if rem(next_value,2) ==0 

next_value = next_value/2; 
else 

next_value = 3*next_value+l; 
end 

seqúense = [sequence, next_value]; 
and 

_lJ 



| collatz. m x collatzplot.m 



Use the document bar to access other 
documents open in the Editor/Debugger. 



collstz 



Ln 1£ 



Col 4 



OVR 



8-13 



O Desktop Tools and Development Environment 



You can use any text editor to créate M-files, such as Emacs. Use preferences 
(accessible from the desktop File menú) to specify that editor as the default. If 
you use another editor, you can still use the MATLAB Editor/Debugger for 
debugging, or you can use debugging functions, such as dbstop, which sets a 
breakpoint. 



If you just need to view the contents of an M-file, you can display the contents 
in the Command Window using the type function. 



For More Information See "Editing and Debugging M-Files" in the 
MATLAB Desktop Tools and Development Environment documentation, and 
the function reference pages for edit, type, and debug. 
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M-Lint Code Check and Profiler Reports 



MATLAB provides tools to help you manage and improve your M-files, 
including the M-Lint Code Check and Profiler Reports. 

M-Lint Code Check Report 

The M-Lint Code Check Report displays potential errors and problems, as well 
as opportunities for improvement in your M-files. The term "lint" is used by 
similar tools in other programming languages such as C. 

Access the M-Lint Code Check Report and other directory reports from the 
Current Directory browser. You run a report for all files in the current 
directory. Alternatively, you can use the mlint function to get results for a 
single M-file. 



Directory reports 



nJfc Current Directory - d:\mymfiles 



File Edit Vlew Debug Desktop Wlndow Help 



d:\mymfiles 



All Files L. 



□ html 
CU mydemos 
Ü seajemp 
ibucky.rn 
icaution.mdl 
i collatz. m 



UJfi#i© 



File Type Last Modified 



Folder 
Folder 

Folder 
M-file 
Model 
M-file 



Mar 3 ,2004.. 
Feb3,2004.. 
Apr5,2004.. 
Nov27,1997. 
Nov13,1997. 
Mar 5 ,2004.. 




TODO/FIXME Report 
© Help Report 
Q Contents Report 
\% Dependency Report 
'Sq File Comparison Report 
O Coverage Report 



£ 



//. 



In MATLAB, the M-Lint Code Check Report displays a message for each line 
of an M-file it determines might be improved. For example, a common M-Lint 
message is that a variable is defined but never used in the M-file. 
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Desktop Tools and Development Environment 



The report 
displays a line 
number and 
message for 
each potential 
problem or 
improvement 
opportunily. 



C lie k a line 
number to operu 
ihe M-fíle ¡n tbe 
Editor at that 
line. 



Í V M-Lint Code Checker Report 



File Edit View Go Debug Desktop Window Help 



HBE 



M Lint Cade Checker Report 



Contents 
No messages 




bucky 

No messages 





lencrthofline 
13 messages 



22: The valué assigned here to variable 

'nothandle ' is never used 

23: HUMEL(x) is usually faster than 

PR0D[5IZE(x)) 

24: Array 'notline' is constructed using 

subscripting. Consider pneallocating fon 

speed 



For More Information See "Tuning and Managing M-Files" in the 
MATLAB Desktop Tools and Development Environment documentation, and 
the reference page for the mlint function. 
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Profiler 

MATLAB includes the Profiler to help you improve the performance of your 
M-files. Run a MATLAB statement or an M-file in the Profiler and it produces 
a report of where the time is being spent. Access the Profiler from the Desktop 
menú, or use the prof ile function. 



»Jt Profiler 



File Edit Debug Desktop Window Help 



*■ -► 



Start Profiling Run this code: collatz 



^« 




Profiletime: seo 



Profiler for Improving Performance 

One way to improve the performance of your M-files is using profiling tools. 
MATLAB pmvldes the M-file Profiler, a graphical user interface that ¡s based on 
the results returned by the prof ile function. Use the Profiler to help you 
determine where you can modify your code to make performance improvements. 

For details on how to use the Profiler, see the Profiler documentation. 



For More Information See "Tuning and Managing M-Files" in the 
MATLAB Desktop Tools and Development Environment documentation, and 
the reference page for the prof ile function. 
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Other Development Environment Features 

Additional development environment features include 

• Source Control — Access your source control system from within MATLAB. 

• Publishing Results — Use the Editor's cell features to publish M-files and 
results to popular output formats including HTML and Microsoft Word. You 
can also use MATLAB Notebook to access MATLAB functions from within 
Microsoft Word. 



For More Information See "Source Control" and "Publishing Results" in 
the MATLAB Desktop Tools and Development Environment documentation. 
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Introducing the Symbolic 
Math Toolbox 



This chapter introduces you to the Symbolic Math Toolbox, and describes how to créate and use 
symbolic objects. It covers the following topics: 



What Is the Symbolic Math Toolbox? 
(p. 9-2) 

Symbolic Objects (p. 9-3) 



Creating Symbolic Variables and 
Expressions (p. 9-5) 

The subs Command (p. 9-8) 

Symbolic and Numeric Conversions 
(p. 9-10) 

Creating Symbolic Math Functions 
(p. 9-15) 



Overview of the toolbox 

Describes symbolic objects and how they differ from 
standard MATLAB data types 

How to créate symbolic objects 



How to substitute numbers for variables 

How to convert between symbolic objects and numeric 
valúes 

How to créate functions that opérate on symbolic objects 
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What Is the Symbolic Math Toolbox? 



The Symbolic Math Toolbox incorporates symbolic computation into the 
numeric environment of MATLAB. The toolbox is a collection of more than 100 
MATLAB functions that provide access to the Maple kernel using a syntax and 
style that is a natural extensión of the MATLAB language. The Symbolic Math 
Toolbox supplements MATLAB numeric and graphical facilities with several 
other types of mathematical computation, which are summarized in following 
table. 



Facility 



Covers 



Calculus 



Linear Algebra 



Simplification 

Solution of 
Equations 



Variable-Precisión 
Arithmetic 

Transforms 



Differentiation, integration, limits, summation, and 
Taylor series 

Inverses, determinants, eigenvalues, singular valué 
decomposition, and canonical forms of symbolic 
matrices 



Methods of simplifying algebraic expressions 

Symbolic and numerical solutions to algebraic and 
differential equations 

Numerical evaluation of mathematical expressions 
to any specified accuracy 

Fourier, Laplace, z-transform, and corresponding 
inverse transforms 



The computational engine underlying the toolboxes is the kernel of Maple, a 
system developed primarily at the University of Waterloo, Canadá and, more 
recently, at the Eidgenossiche Technische Hochschule, Zürich, Switzerland. 
Maple is marketed and supported by Waterloo Maple, Inc. 
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Symbolíc Objects 



The Symbolic Math Toolbox defines a new MATLAB data type called a 
symbolic object. (See "Programming and Data Types" in the MATLAB 
documentation for an introduction to MATLAB classes and objects.) Internally, 
a symbolic object is a data structure that stores a string representation of the 
symbol. The Symbolic Math Toolbox uses symbolic objects to represent 
symbolic variables, expressions, and matrices. The actual computations 
involving symbolic objects are performed primarily by Maple, mathematical 
software developed by Waterloo Maple, Inc. 

The following example illustrates the difference between a standard MATLAB 
data type, such as double, and the corresponding symbolic object. The 
MATLAB command 

sqrt(2) 

returns a floating-point decimal number: 

ans = 

1 .4142 

On the other hand, if you convert 2 to a symbolic object using the sym command, 
and then take its square root by entering 

a = sqrt (sym(2) ) 

the result is 

a = 
2"(1/2) 

MATLAB gives the result 2" (1 12), which means 2 , using symbolic notation 
for the square root operation, without actually calculating a numerical valué. 
MATLAB records this symbolic expression in the string that represents 
2 " ( 1 / 2 ) . You can always obtain the numerical valué of a symbolic object with 
the double command: 

double(a) 
ans = 

1 .4142 

Notice that the result is indented, which tells you it has data type double. 
Symbolic results are not indented. 
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When you créate a fraction involving symbolic objects, MATLAB records the 
numerator and denominator. For example: 

sym(2) /sym(5) 

ans = 

2/5 

MATLAB performs arithmetic on symbolic objects differently than it does on 
standard data types. If you add two fractions that are of data type double, 
MATLAB gives the answer as a decimal fraction. For example: 

2/5 + 1/3 
ans = 
0.7333 

If you add the same fractions as symbolic objects, MATLAB finds their common 
denominator and combines them by the usual procedure for adding rational 
numbers: 

sym(2) /sym(5) + sym(1 ) /sym(3) 

ans = 

11/15 

The Symbolic Math Toolbox enables you to perform a variety of symbolic 
calculations that arise in mathematics and science. These are described in 
detail in Chapter 10, "Using the Symbolic Math Toolbox." 
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Creatíng Symbolic Variables and Expressions 

The sym command lets you construct symbolic variables and expressions. For 
example, the commands 

x = sym( 'x ' ) 

a = sym( ' alpha 1 ) 

créate a symbolic variable x that prints as x and a symbolic variable a that 
prints as alpha. 

Suppose you want to use a symbolic variable to represent the golden ratio 

1±V5 
y 2 

The command 

rho = sym('(1 + sqrt(5))/2') 

achieves this goal. Now you can perform various mathematical operations on 
rho. For example, 

f = rho"2 - rho - 1 

returns 

f = 

(1/2+1/2*5" (1/2)) "2 -3/2 -1/2*5" (1/2) 

You can simplify this answer by entering 
simplify(f) 

which returns 

ans = 


2 
Now suppose you want to study the quadratic function f = ax +bx + c . One 

approach is to enter the command 
f = sym( 'a*x"2 + b*x + c ) 
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2 
which assigns the symbolic exprés sion ax + bx + c to the variable f . However, 

in this case, the Symbolic Math Toolbox does not créate variables 

corresponding to the terms of the exprés sion, a , b , c , and x . To perform 

symbolic math operations (e.g., integration, differentiation, substitution, etc.) 

on f , you need to créate the variables explicitly. A better alternative is to enter 

the commands 

a = sym( ' a ' ) 
b = sym( ' b ' ) 
c = sym( ' c ' ) 
x = sym( ' x ' ) 

or simply 

syms a b c x 

Then enter 

f = sym( 'a*x"2 + b*x + c' ) 



In general, you can use sym or syms to créate symbolic variables. We 
recommend you use syms because it requires less typing. 



Note To créate a symbolic exprés sion that is a constant, you must use the 
sym command. For example, to créate the expression whose valué is 5, enter 
f = sym ( ' 5 ' ) . Note that the command f = 5 does not define f as a symbolic 
expression. 



If you set a variable equal to a symbolic expression, and then apply the syms 
command to the variable, MATLAB removes the previously defined expression 
from the variable. For example, 

syms a b 
f = a + b 

returns 

f = 
a+b 
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If you then enter 

syms f 
f 

MATLAB returns 

f = 
f 

You can use the syms command to clear variables of definitions that you 
assigned to them previously in your MATLAB session. However, syms does not 
clear the properties of the variables in the Maple workspace. See "Clearing 
Variables in the Maple Workspace" on page 9-12 for more information. 

The fíndsym Command 

To determine what symbolic variables are present in an expression, use the 
f indsym command. For example, given the symbolic expressions f and g 
defined by 

syms a b n t x z 

f = x"n; g = sin(a*t + b) ; 

you can find the symbolic variables in f by entering 

f indsym(f ) 
ans = 
n, x 

Similarly, you can find the symbolic variables in g by entering 

f indsym(g) 
ans = 
a, b, t 
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The subs Command 



You can substitute a numerical valué for a symbolic variable using the subs 
command. For example, to substitute the valué % = 2 in the symbolic 
exprés sion, 

f = 2*x~2 - 3*x + 1 

enter the command 
subs(f ,2) 

This returns f(2): 
ans = 



Note To substitute a matrix A into the symbolic expression f , use the 
command polyvalm(sym2poly(f ) , A) , which replaces all occurrences of xby 
A, and replaces the constant term of f with the constant times an identity 
matrix. 



When your expression contains more than one variable, you can specify the 
variable for which you want to make the substitution. For example, to 
substitute the valué x = 3 in the symbolic expression, 

syms x y 

f = x"2*y + 5*x*sqrt(y) 

enter the command 
subs(f, x, 3) 

This returns 

ans = 
9*y+15*y~(1/2) 
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On the other hand, to substitutey = 3, enter 

subs(f, y, 3) 
ans = 
3*x~2+5*x*3~(1/2) 

The Default Symbolíc Variable 

If you do not specify a variable to substitute for, MATLAB chooses a default 
variable according to the following rule. For one-letter variables, MATLAB 
chooses the letter closest to x in the alphabet. If there are two letters equally 
cióse to x, MATLAB chooses the one that comes later in the alphabet. In the 
preceding function, subs (f , 3) returns the same answer as 
subs(f, x, 3). 

You can use the f indsym command to determine the default variable. For 
example, 

syms s t 
g = s + t; 
f indsym(g,1 ) 

returns the default variable: 

ans = 
t 

See "Substitutions" on page 10-50 to learn more about substituting for 
variables. 
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Symbolic and Numeric Conversions 



Consider the ordinary MATLAB quantity 
t = 0.1 

The sym function has four options for returning a symbolic representation of 
the numeric valué stored in t. The ' f ' option 

sym(t, 'f ' ) 

returns a symbolic floating-point representation 
'1 .999999999999a' *2~( -4) 

The ' r ' option 
sym(t, 'r' ) 

returns the rational form 
1/10 

This is the default setting for sym. That is, calling sym without a second 
argument is the same as using sym with the ' r ' option: 

sym(t) 

ans = 
1/10 

The third option ' e ' returns the rational form of t plus the difference between 
the theoretical rational expression for t and its actual (machine) floating-point 
valué in terms of eps (the floating-point relative accuracy): 

sym(t, 'e' ) 

ans = 
1/10+eps/40 

The fourth option ' d ' returns the decimal expansión of t up to the number of 
significant digits specified by digits: 

sym(t, 'd' ) 
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ans = 
.10000000000000000555111512312578 

The default valué of digits is 32 (henee, sym(t , ' d ' ) returns a numberwith 32 
significant digits), but if you prefer a shorter representation, use the digits 
command as follows: 

digits(7) 

sym(t, 'd' ) 

ans = 
.1000000 

A particularly effective use of sym is to convert a matrix from numeric to 
symbolic form. The command 

A = hilb(3) 
generates the 3-by-3 Hilbert matrix: 



1.0000 0, 


.5000 





.3333 


0.5000 0. 


.3333 





.2500 


0.3333 0, 


.2500 





.2000 


applying sym to A 








A = sym(A) 









you can obtain the symbolic (infinitely precise) form of the 3-by-3 Hilbert 
matrix: 



[ 1, 1/2, 1/3] 
[ 1/2, 1/3, 1/4] 
[ 1/3, 1/4, 1/5] 



9-11 



y Introducing the Symbolic Math Toolbox 



Constructíng Real and Complex Variables 

The sym command allows you to specify the mathematical properties of 
symbolic variables by using the ' real ' option. That is, the statements 

x = sym( ' x ' , ' real ' ) ; y = sym( ' y ' , ' real ' ) ; 

or more efficiently 

syms x y real 
z = x + i*y 

créate symbolic variables x and y that have the added mathematical property 
of being real variables. Specifically this means that the expression 

f = x"2 + y"2 

is strictly nonnegative. Henee, z is a complex variable and can be manipulated 
as such. Thus, the commands 

conj(x), conj(z), expand(z*conj (z) ) 

return 

x, x-i*y, x"2+y"2 

respectively. The con j command is the complex conjúgate operator for the 
toolbox. If conj (x) == x returns 1, then x is a real variable. 

Clearing Variables in the Maple Workspace 

When you declare x to be real with the command 

syms x real 

a becomes a symbolic object in the MATLAB workspace and a positive real 
variable in the Maple kernel workspace. If you later want to remove the real 
property from x, enter 

syms a unreal 

Note that entering 
clear x 

only clears x in the MATLAB workspace. If you then enter syms x, MATLAB 
still treats x as a positive real number. 
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Creatíng Abstract Functíons 

If you want to créate an abstract (i.e., indeterminant) function f(x) , type 
f = sym( 'f (x) ' ) 

Then f acts like fix) and can be manipulated by the toolbox commands. For 
example, to construct the first difference ratio, type 

df = (subs(f , 'x' , 'x+h' ) - f )/'h' 
or 

syms x h 

df = (subs(f ,x,x+h)-f )/h 

which returns 

df = 

(f(x+h)-f (x))/h 

This application of sym is useful when computing Fourier, Laplace, and 
z-transforms. 

Using sym to Access Maple Functions 

Similarly, you can access Maple's factorial function k! using sym: 

kf ac = sym( ' k! ' ) 

To compute 6! or n!, type 

syms k n 

subs(kf ac, k,6) , subs(kf ac,k,n) 

ans = 
720 

ans = 

n! 

Example: Creating a Symbolic Matrix 

A circulant matrix has the property that each row is obtained from the previous 
one by cyclically permuting the entries one step forward. You can créate the 
circulant matrix A whose elements are a, b, and c, using the commands 
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syms abe 

A= [abe; bea; cab] 

which return 

A = 

[ a, b, c ] 

[ b, c, a ] 

[ c, a, b ] 

Since A is circulant, the sum over each row and column is the same. To check 
this for the first row and second column, enter the command 

sum(A(1 , : ) ) 

which returns 

ans = 
a+b+c 

The command 

sum(A(1,:)) == sum(A(:,2)) % This is a logical test. 

returns 

ans = 
1 

Now replace the (2,3) entry of A with beta and the variable b with alpha. The 
commands 

syms alpha beta; 

A(2,3) = beta; 

A = subs(A,b, alpha) 

return 

A = 

[ a, alpha, c] 

[ alpha, c, beta] 

[ c, a, alpha] 

From this example, you can see that using symbolic objeets is very similar to 
using regular MATLAB numeric objeets. 
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Creatíng Symbolic Math Functions 

There are two ways to créate functions: 



• Use symbolic expressions 

• Créate an M-file 

Usíng Symbolic Expressions 

The sequence of commands 

syms x y z 

r = sqrt(x"2 + y"2 + z"2) 

t = atan(y/x) 

f = sin(x*y) / (x*y) 

generates the symbolic expressions r, t, and f . You can use dif f , int, subs, 
and other Symbolic Math Toolbox functions to manipúlate such expressions. 

Creating an M-File 

M-files permit a more general use of functions. Suppose, for example, you want 
to créate the sinc function sin (x) /x. To do this, créate an M-file in the @sym 
directory: 

function z = sinc(x) 

%SINC The symbolic sinc function 

% sin(x)/x. This function 

% accepts a sym as the input argument. 

if isequal(x,sym(0) ) 

z = 1; 
else 

z = sin(x)/x; 
end 

You can extend such examples to functions of several variables. See 
"Programming and Data Types" in the MATLAB documentation for a more 
detailed discussion on object-oriented programming. 
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Using the Symbolic Math 
Toolbox 



This chapter explains how to use the Symbolic Math Toolbox to perform many common mathematical 
operations. It covers the following topics: 



Calculus (p. 10-2) 

Simplifications and Substitutions 
(p. 10-41) 

Variable-Precisión Arithmetic 
(p. 10-57) 

Linear Algebra (p. 10-62) 
Solving Equations (p. 10-86) 



Differentiation, integration, limits, summation, and 
Taylor series 

Methods of simplifying algebraic expressions 

Numerical evaluation of mathematical expressions to any 
specified accuracy 

Inverses, determinants, eigenvalues, singular valué 
decomposition, and canonical forms of symbolic matrices 

Symbolic and numerical solutions to algebraic and 
differential equations 
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Calculus 



The Symbolic Math Toolbox provides functions to do the basic operations of 
calculus. The following sections describe these functions: 

• "Differentiation" on page 10-2 

• "Limits" on page 10-8 

• "Integration" on page 10-11 

• "Symbolic Summation" on page 10-18 

• "Taylor Series" on page 10-18 

• "Calculus Example" on page 10-20 

• "Extended Calculus Example" on page 10-28 

Differentiation 

To illustrate how to take derivatives using the Symbolic Math Toolbox, first 
créate a symbolic expression: 

syms x 

f = sin(5*x) 

The command 

diff (f) 

differentiates f with respect to x: 

ans = 
5*cos(5*x) 

As another example, let 
g = exp(x) *cos(x) 

where exp(x) denotes e x , and differentiate g: 

diff(g) 

ans = 

exp(x) *cos(x)-exp(x)*sin(x) 
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To take the second derivative of g, enter 

diff(g,2) 

ans = 

-2*exp(x) *sin(x) 

You can get the same result by taking the derivative twice: 

diff (diff(g)) 
ans = 
-2*exp(x)*sin(x) 

In this example, MATLAB automatically simplifies the answer. However, in 
some cases, MATLAB might not simply an answer, in which case you can use 
the simplif y command. For an example of this, see "More Examples" on 
page 10-5. 

Note that to take the derivative of a constant, you must first define the 
constant as a symbolic expression. For example, entering 

c = sym( '5 ' ) ; 
diff (c) 

returns 

ans = 



If you just enter 
diff (5) 

MATLAB returns 

ans = 
[] 

because 5 is not a symbolic expression. 

Derivatives of Expressions with Several Variables 

To differentiate an expression that contains more than one symbolic variable, 
you must specify the variable that you want to differentiate with respect to. 
The diff command then calculates the partial derivative of the expression 
with respect to that variable. For example, given the symbolic expression 
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syms s t 

f = sin(s*t) 

the command 

diff (f ,t) 

calculates the partial derivative df/dt . The result is 

ans = 
cos(s*t) *s 

To differentiate f with respect to the variable s, enter 
diff (f ,s) 

which returns: 

ans = 
cos(s*t) *t 

If you do not specify a variable to differentiate with respect to, MATLAB 
chooses a default variable by the same rule described in "The subs Command" 
on page 9-8. For one-letter variables, the default variable is the letter closest to 
x in the alphabet. In the preceding example, diff ( f ) takes the derivative of f 
with respect to t because t is closer to x in the alphabet than s is. To determine 
the default variable that MATLAB differentiates with respect to, use the 
f indsym command: 

f indsym(f , 1 ) 
ans = 

t 

To calcúlate the second derivative of f with respect to t, enter 
diff(f,t,2) 

which returns 

ans = 
-sin(s*t)*s"2 

Note that diff ( f , 2 ) returns the same answer because t is the default 
variable. 
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More Examples 

To further illustrate the dif f command, define a, b, x, n, t, and theta in the 
MATLAB workspace by entering 

syms a b x n t theta 
The table below illustrates the results of entering dif f ( f ) . 



f diff(f) 


x"n x"n*n/x 


sin(a*t+b) 


cos(a*t+b) *a 


exp(i*theta) i*exp(i*theta) 



In the first example, MATLAB does not automatically simplify the answer. To 
simplify the answer, enter 

simplify(diff (x"n) ) 
ans = 
x"(n-1 )*n 

To differentiate the Bessel function of the first kind, besselj (nu,z), with 
respect to z, type 

syms nu z 

b = besselj (nu,z) ; 

db = diff(b) 

which returns 

db = 

-besselj (nu+1 ,z)+nu/z*besselj (nu,z) 

The dif f function can also take a symbolic matrix as its input. In this case, the 
differentiation is done element-by-element. Consider the example 

syms a x 

A = [cos(a*x) , sin(a*x) ; -sin(a*x) ,cos(a*x) ] 
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which returns 
A = 



[ 



cos(a*x) 
-sin(a*xl 



sin(a*x) ] 
cos(a*x) ] 



The command 

diff (A) 

returns 

ans = 

[ -sin(a*x) 

[ -cos(a*x) 



cos(a*x) 
■sin(a*x) 



You can also perform differentiation of a column vector with respect to a row 
vector. Consider the transformation from Euclidean (x, y, z) to spherical 
(r, X, (p) coordinates as given by x = rcos^coscp , y = rcos^sincp, and 
z = rsin^.. Note that X corresponds to elevation or latitude while (p denotes 
azimuth or longitude. 




To calcúlate the Jacobian matrix, J, of this transformation, use the j acobian 
function. The mathematical notation for J is 

j _ d(x,y,x) 
d(r, X, (p) 

For the purposes of toolbox syntax, use 1 for X and f for (p . The commands 
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syms r 1 f 

x = r*cos(l) *cos(f ) ; y = r*cos(l)*sin(f ) ; z = r*sin(l); 

J = jacobian([x; y; z], [r 1 f]) 

return the Jacobian 

J = 

[ cos(l) *cos(f ) , -r*sin(l)*cos(f ) , -r*cos(l) *sin(f ) ] 

[ cos(l) *sin(f ) , -r*sin(l)*sin(f ) , r*cos(l) *cos(f ) ] 

[ sin(l), r*cos(l), 0] 

and the command 

detJ = simple(det (J) ) 

returns 

detJ = 
-cos(l)*r"2 

Notice that the first argument of the j acobian function must be a column 
vector and the second argument a row vector. Moreover, since the determinant 
of the Jacobian is a rather complicated trigonometric exprés sion, you can use 
the simple command to make trigonometric substitutions and reductions 
(simplifications). See "Simplifications and Substitutions" on page 10-41 for 
more details on simplification. 

Atable summarizing diff and jacobian follows. 



Mathematical Operator 


MATLAB Command 


df 

dx 


diff (f ) or diff (f ,x) 


df 

da 


diff (f, a) 
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Mathematical Operator 



MATLAB Command 



d ¿ f 

db 2 


diff (f ,b,2) 


j _ d(r,í) 

d(u, v) 


J = jacobian( [ r;t ] , 


[u,v]) 



Límíts 

The fundamental idea in calculus is to make calculations on functions as a 
variable "gets cióse to" or approaches a certain valué. Recall that the definition 
of the derivative is given by a limit 

f Xx ) = lim f(x + h)-f(x) 
h->0 h 

provided this limit exists. The Symbolic Math Toolbox enables you to calcúlate 
the limits of functions directly. The commands 

syms h n x 

limit( (cos(x+h) - cos(x) ) /h,h,0 ) 

which return 

ans = 
-sin(x) 

and 

limit( (1 + x/n)"n,n,inf ) 

which returns 

ans = 
exp(x) 

illustrate two of the most important limits in mathematics: the derivative (in 
this case of eos x) and the exponential function. 



10-8 



Calculi 



One-Sided Limits 

You can also calcúlate one-sided limits with the Symbolic Math Toolbox. For 
example, you can calcúlate the limit of xl\x | , whose graph is shown in the 
following figure, as x approaches from the left or from the right. 



1 - 



x/abs(x) 



0.5- 



-0.5 



-0.5 0.5 

x 



To calcúlate the limit as x approaches from the left, 



x-»0~ W 

enter 

limit (x/abs(x) ,x,0, 'left ' 

This returns 

ans = 
-1 
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To calcúlate the limit as x approaches from the right, 

lim — = 1 

x^o + M 

enter 

limit (x/abs(x) , x,0, ' right ' ) 

This returns 

ans = 

1 

Since the limit from the left does not equal the limit from the right, the two- 
sided limit does not exist. In the case of undefined limits, MATLAB returns NaN 
(not a number). For example, 

limit (x/abs(x) ,x,0) 

returns 

ans = 
NaN 

Observe that the default case, limit ( f ) is the same as limit ( f , x , ) . Explore 
the options for the limit command in this table, where f is a function of the 
symbolic object x. 



Mathematical Operation 



MATLAB Command 



lim f(x) 

x^O 



lim f(x) 

x — ¥ a 

lim f(x) 



limit(f) 



limit (f,x,a) or 
limit (f, a) 

limit (f ,x,a, ' left ' ) 



lim f(x) 

x — > a+ 



limit (f ,x,a, ' right ' 
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Integratíon 



If f is a symbolic expression, then 
int(f) 

attempts to find another symbolic expression, F, so that dif f ( F) = f . That is, 
int (f ) returns the indefinite integral or antiderivative of f (provided one 
exists in closed form). Similar to differentiation, 

int(f,v) 

uses the symbolic object v as the variable of integra tion, rather than the 
variable determined by f indsym. See how int works by looking at this table. 



Mathematical Operation 



MATLAB Command 



g = eos (at + b) 
\g(t)dt = sin(at + b)/a 

jj 1 (z)dz = -J (z) 



re + 1 

\x dx = 

J n + 1 


int (x"n) or 
int (x"n,x) 


71/2 

f sin(2x)c¿x = 1 



int(sin(2*x) ,0,pi/2) or 
int(sin(2*x) ,x,0,pi/2) 



g = cos(a*t + b) 
int(g) or 
int(g,t) 

int (besselj (1 ,z) ) or 
int (besselj (1 ,z) ,z) 



In contrast to differentiation, symbolic integration is a more complicated task. 
A number of difficulties can arise in computing the integral: 



• The antiderivative, F, may not exist in closed form. 

• The antiderivative may define an unfamiliar function. 



The antiderivative may exist, but the software can't find the it. 

The software could find the antiderivative on a larger computer, but runs out 
of time or memory on the available machine. 
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Nevertheless, in many cases, MATLAB can perform symbolic integration 
successfully. For example, créate the symbolic variables 



syms a b theta x y n u z 

The following table illustrates integration of exprés sions containing those 
variables. 



f 




int(f) 




x"n 




x~(n+1 )/(n+1) 




y(-i) 




log(y) 




n"x 


1 /log(n) *n"x 




sin(a* 


theta+b) 


-1 /a*cos(a*theta+b) 


1/(1+u 


-2) 


atan(u) 


exp( -x 


~2) 


1/2*pi"(1/2)*erf (x) 



In the last example, exp( -x"2), there is no formula for the integral involving 
standard calculus expressions, such as trigonometric and exponential 
functions. In this case, MATLAB returns an answer in terms of the error 
function erf . 

If MATLAB is unable to find an answer to the integral of a function f , it just 
returns int (f ). 

Definite integration is also possible. The commands 

int (f ,a,b) 

and 

int (f , v,a,b) 

are used to find a symbolic expression for 



f f(x)dx and f f(v)dv 

a a 

respectively. 
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f 


a, 


b 


int(f,a,b) 






X"7 


o, 


1 


1/8 






1/x 


1, 


2 
1 


iog(2) 






log(x) *sqrt (x) 


o, 


-4/9 






exp( -x~2) 


o, 


inf 


1/2*pi"(1/2) 






besselj (1 ,z) "2 


o, 


1 


1 /12*hypergeom( [3/2, 
[2, 5/2, 3],-1) 


3/2], 





For the Bessel function (besselj ) example, it is possible to compute a 
numerical approximation to the valué of the integral, using the double 
function. The commands 

syms z 

a = int(besselj (1 ,z) "2,0,1 ) 

return 

a = 

1/12*hypergeom([3/2, 3/2], [2, 5/2, 3],-1) 

and the command 

a = double(a) 

returns 

a = 
0.0717 

Integration with Real Parameters 

One of the subtleties involved in symbolic integration is the "valué" of various 
parameters. For example, if a is any positive real number, the expression 
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is the positive, bell shaped curve that tends to as x tends to +°° . You can 
créate an example of this curve, for a = 1/2 , using the following commands: 



syms x 

a = sym(1 /2) ; 

f = exp( -a*x"2) ; 

ezplot(f) 



exp(-1/2x ) 




However, if you try to calcúlate the integral 



f — ax , 

e ax 



without assigning a valué to a, MATLAB assumes that a represents a complex 
number, and therefore returns a complex answer. If you are only interested in 
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the case when ais a positive real number, you can calcúlate the integral as 
follows: 

syms a positive; 

The argument positive in the syms command restricts a to have positive 
valúes. Now you can calcúlate the preceding integral using the commands 

syms x; 

f = exp( -a*x"2) ; 

int (f ,x, -inf ,inf ) 

This returns 

ans = 
1/(aP(1/2)*pi~(1/2) 

If you want to calcúlate the integral 

2 

e dx 

for any real number a, not necessarily positive, you can declare a to be real with 
the following commands: 

syms a real 

f =exp( -a*x"2) ; 

F = int(f, x, -inf, inf) 

MATLAB returns 

F = 

PIECEWISE([1/a~(1/2)*pi~(1/2) , signum(a) =1],[Inf, otherwise]) 
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You can put this in a more readable form by entering 
pretty (F) 



1/2 



Pi 



1/2 



Inf 



signum(a-) = 1 



otherwise 



The - after a is simply a reminder that a is real, and signum ( a- ) is the sign of 
a. So the integral is 

Ja 
when a is positive, just as in the preceding example, and » when a is negative. 

You can also declare a sequence of symbolic variables w, y, x, z to be real by 
entering 

syms w x y z real 

Integration with Complex Parameters 

To calcúlate the integral 



-oí , 

e ax 



for complex valúes of a, enter 

syms a x unreal % 

f = exp( -a*x"2) ; 

F = int(f, x, -inf, inf) 

Note that syms is used with the unreal option to clear the real property that 
was assigned to a in the preceding example — see "Clearing Variables in the 
Maple Workspace" on page 9-12. 
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The preceding commands produce the complex output 

F = 

PIECEWISE([1/a~(1/2)*pi~(1/2) , csgn(a) =1],[Inf, otherwise] 

You can make this output more readable by entering 

pretty (F) 

1/2 



Pi 

1/2 
a 

Inf 



csgn(a) 



otherwise 



The expression csgn (a) (complex sign of a) is defined by 

csgn(a) = í 1 ifRe(a)>0, orRe(a) = OandIm(a)>0 
1-1 if Re(a) < 0, or Re(a) = and Im(a) < 

The condition csgn(a) = 1 corresponds to the shaded región of the complex plañe 
shown in the following figure. 




Región corresponding to csgn(a) = 1 



The square root of a in the answer is the unique square root lying in the shaded 
región. 
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Symbolic Summatíon 

You can compute symbolic summations, when they exist, by using the symsum 
command. For example, the p-series 

2 2 3 2 
sums to n 2 /6 , while the geometric series 

1 2 

1 + X + X + ... 

sums to 1/(1 -x) , provided \x\ < 1 . Three summations are demonstrated 
below: 

syms x k 

s1 = symsum ( 1 /k"2,1 ,inf) 

s2 = symsum(x"k,k,0,inf ) 

s1 = 

1/6*pi~2 

s2 = 

-1/(x-1) 

Taylor Seríes 

The statements 

syms x 

f = 1/(5+4*cos(x)) 

T = taylor(f,8) 

return 

T = 
1/9+2/81*x~2+5/1458*x~4+49/131220*x"6 

which is all the terms up to, but not including, order eight in the Taylor series 
for f{x): 
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ni 



re = 

Technically, T is a Maclaurin series, since its basepoint is a = 0. 
The command 
pretty(T) 

prints T in a format resembling typeset mathematics: 

2 4 49 6 

1/9 + 2/81 x + 5/1458 x + x 

131220 

These commands 

syms x 

g = exp(x*sin(x) ) 

t = taylor(g,12,2) ; 

genérate the first 12 nonzero terms of the Taylor series for g about x = 2. 

Next, plot these functions together to see how well this Taylor approximation 
compares to the actual function g: 

xd = 1:0.05:3; yd = subs(g,x,xd) ; 

ezplot(t, [1,3]); hold on; 

plot (xd, yd, ' r- . ' ) 

title( 'Taylor approximation vs. actual function'); 

legend( 'Taylor' , 'Function' ) 
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Taylor approximation vs. actual function 




Special thanks to Professor Gunnar Backstr0m of UMEA in Sweden for this 
example. 

Calculus Example 

This section describes how to analyze a simple function to find its asymptotes, 
máximum, minimum, and inflection point. The section covers the following 
topics: 

• "Defining the Function" on page 10-21 

• "Finding the Asymptotes" on page 10-22 

• "Finding the Máximum and Minimum" on page 10-24 

• "Finding the Inflection Point" on page 10-26 



10-20 



Calculi 



Defining the Function 

The function in this example is 

Cl n 3x + 6x - 1 

X tX — O 

To créate the function, enter the following commands: 

syms x 

num = 3*x"2 + 6*x -1 ; 
denom = x"2 + x - 3; 
f = num/denom 

This returns 

f = 

(3*x"2+6*x-1 )/(x~2+x-3) 

You can plot the graph of f by entering 
ezplot(f) 
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This displays the following plot. 







(3 x 2 +6 x- 


1)/(x 2 +x-3) 


8 


- 




\ 


6 


- 




\ 


4 


: J 




^—~ 


2 




- 





- 




\ 


2 


- 




\ 


4 









Finding the Asymptotes 

To find the horizontal asymptote of the graph of f , take the limit of f as x 
approaches positive infinity: 

limit(f, inf) 

ans = 

3 

The limit as x approaches negative infinity is also 3. This tells you that the line 
y = 3 is a horizontal asymptote to the graph. 

To find the vertical asymptotes of f , set the denominator equal to and solve 
by entering the following command: 

roots = solve(denom) 
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2 
This returns to solutions to x + x - 3 = : 

roots = 

[ -1/2+1/2*13^(1/2)] 

[ -1/2-1/2*13^(1/2)] 

This tells you that vertical asymptotes are the lines 
-1 + VÍ3 

and 

_ -1-7Í3 
2 

You can plot the horizontal and vertical asymptotes with the following 
commands: 

ezplot(f) 

hold on % Keep the graph of f in the figure 
% Plot horizontal asymptote 
plot([-2*pi 2*pi], [3 3] , 'g' ) 
% Plot vertical asymptotes 
plot(double(roots(1))*[1 1], [-5 10],'r') 
plot(double(roots(2) )*[1 1], [-5 10],'r') 
title( ' Horizontal and Vertical Asymptotes') 
hold off 

Note that roots must be converted to double to use the plot command. 
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The preceding commands display the following figure. 



Horizontal and Vertical Asymptotes 



8 
6 

4 


l 




v_ 




_^y 


\\ 


- 


2 

2 
4 





To recover the graph of f without the asymptotes, enter 
ezplot(f) 

Finding the Máximum and Minimum 

You can see from the graph that f has a local máximum somewhere between 
the points x = 2 and x = 3, and might have a local minimum between x = -4 and 
x = -2. To find the x-coordinates of the máximum and minimum, first take the 
derivative of f : 

f1 = diff(f) 

This returns 

f1 = (6*x+6)/(x~2+x-3)-(3*x~2+6*x-1 ) / (x~2+x-3p2* (2*x+1 ) 

To simplify this expression, enter 
f1 = simplify(f1 ) 
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which returns 

f1 = -(3*x~2+16*x+17)/(x~2+x-3P2 

You can display f 1 in a more readable form by entering 
pretty(fl) 

which returns 

2 
3 x + 16 x + 17 



2 2 

¡x + x - 3) 



Next, set the derivative equal to and solve for the critical points: 
crit_pts = solve(fl) 

This returns 

ans = 

[ -8/3-1/3*13^(1/2)] 

[ -8/3+1/3*13^(1/2)] 

It is clear from the graph of f that it has a local minimum at 
_ -8-7Í3 

X i — "~ 



and a local máximum at 

_ -8 + VÍ3 

iA/Q — "~ 



Note MATLAB does not always return the roots to an equation in the same 
order. 



You can plot the máximum and minimum of f with the following commands: 

ezplot(f) 
hold on 
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plot (double(crit_pts) , double(subs(f ,crit_pts¡ 
title( 'Máximum and Minimum of f ' ) 
text( -5.5,3.2, ' Local mínimum') 
text( -2.5,2, ' Local máximum') 
hold off 

This displays the following figure. 

Máximum and Minimum of f 



ro 




Finding the Inflection Point 

To find the inflection point of f , set the second derivative equal to and solve. 

f2 = diff (f1) ; 
inflec_pt = solve(f2); 
double(inf lec_pt) 

This returns 



ans 



■5.2635 

■1 .3682 - 0.851U 

-1 .3682 + 0.8511Í 



10-26 



Calculi 



In this example, only the first entry is a real number, so this is the only 
inflection point. (Note that in other examples, the real solutions might not be 
the first entries of the answer.) Since you are only interested in the real 
solutions, you can discard the last two entries, which are complex numbers. 

inflec_pt = inflec_pt(1) 

To see the symbolic expression for the inflection point, enter 
pretty (simplify(inf lec_pt) ) 

This returns 

1/2 2/3 1/2 1/3 

(676 + 156 13 ) +52 + 16 (676 + 1 56 1 3 ) 

1/6 

1/2 1/3 
(676 + 156 13 ) 

To plot the inflection point, enter 

ezplot(f , [-9 6]) 

hold on 

plot (double(inf lec_pt) , double(subs(f ,inf lec_pt) ) , ' ro ' ) 

title(' Inflection Point of f) 

text ( -7,2, ' Inflection point') 

hold off 

The extra argument, [-9 6],inezplot extends the range of % valúes in the plot 
so that you see the inflection point more clearly, as shown in the following 
figure. 
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Inflection Point of f 



-0- 



Inflection point 



-8-6-4-2 2 4 6 

x 



Extended Calculus Example 

This section presents an extended example that illustrates how to find the 
máxima and mínima of a function. The section covers the following topics: 

• "Defining the Function" on page 10-28 

• "Finding the Zeros of f3" on page 10-30 

• "Finding the Máxima and Minima of f2" on page 10-33 

• "Integrating" on page 10-35 

Defining the Function 

The starting point for the example is the function 

ftx) = 



5 + 4cos(x) 
You can créate the function with the commands 

syms x 

f = 1/(5+4*cos(x)) 
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which return 

f = 

1/(5+4*COS(X)) 

The example shows how to find the máximum and mínimum of the second 
derivative of fix) . To compute the second derivative, enter 



f2 = diff(f,2) 

which returns 

f2 = 
32/(5+4*cos(x))~3*sin(x)"2+4/(5+4*cos(x))~2*cos(x) 

Equivalently, you can type f2 = diff (f ,x,2). The default scaling in ezplot 
cuts off part of the graph of f 2. You can set the axes limits manually to see the 
entire function: 

ezplot(f2) 

axis( [ -2*pi 2*pi -5 2]) 
title( 'Graph of f 2 ' ) 
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From the graph, it appears that the máximum valué of f"(x) is 1 and the 
mínimum valué is -4. As you will see, this is not quite true. To find the exact 
valúes of the máximum and minimum, you only need to find the máximum and 
minimum on the interval (-% %]. This is true because f"(x) is periodic with 
period 2k, so that the máxima and minima are simply repeated in each 
translation of this interval by an integer múltiple of 2% . The next two sections 
explain how to do find the máxima and minima. 



Finding the Zeros of f3 

The máxima and minima of f"{x) occur at the zeros of f'"{x) . The statements 

f3 = diff (f2) ; 
pretty(f3) 

compute f'"(x) and display it in a more readable form: 

3 

sin(x) cos(x) 



sin(x) 

384 - - - + 96 - 

4 3 

(5+4 COS(X) ) (5+4 COS(X) ) 



sin x 



[5+4 cos(x) 
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You can simplify this expression using the statements 

f3 = simple(f3) ; 
pretty(f3) 

2 2 
sin(x) (96 sin(x) + 80 cos(x) + 80 cos(x) - 25) 
4 --- 

4 
(5+4 COS(X) ) 

Now, to find the zeros of f"'(x) , enter 
zeros = solve(f3) 

This returns a 5-by-l symbolic matrix 

zeros = 

0] 

atan ((-255-60*1 9" (1/2)) "(1/2) , 1 0+3*19* (1 /2) ) ] 

atan(-( -255-60*1 9* (1/2))* (1/2) , 10+3*19* (1 /2) ) ] 

atan( (-255+60*19* (1/2))* (1/2)/ (10-3*19* (1/2)))+pi] 

-atan(( -255+60*19* (1/2))* (1/2)/ (10-3*19* (1/2))) -pi] 

each of whose entries is a zero of f'"(x) . The commands 

format; % Default format of 5 digits 
zerosd = double(zeros) 

convert the zeros to double form: 

zerosd = 


0+ 2.4381Í 
0- 2.4381Í 
2.4483 
-2.4483 
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So far, you have found three real zeros and two complex zeros. However, as the 
following graph of f3 shows, these are not all its zeros: 

ezplot(f3) 

hold on; 

plot (zerosd,0*zerosd, ' ro ' ) % Plot zeros 

plot([-2*pi,2*pi], [0,0] , 'g- . ' ) ; % Plot x-axis 

title( 'Graph of f 3 ' ) 



Graph of f3 






The red circles in the graph correspond to zerosd(1 ), zerosd(4), and 
zerosd (5 ) . As you can see in the graph, there are also zeros at +7t . The 
additional zeros occur because f'"(x) contains a factor of sin(x) , which is zero 
at integer múltiples of k . The function, solve ( sin ( x ) ) , however, only finds the 
zero atx = 0. 

A complete list of the zeros of f"\%) in the interval (-K n] is 

zerosd = [zerosd(1) zerosd(4) zerosd(5) pi] ; 
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You can display these zeros on the graph of f'"(x) with the following 
commands: 

ezplot(f3) 

hold on; 

plot(zerosd,0*zerosd, 'ro' ) 

plot( [-2*pi,2*pi] , [0,0] , 'g- . ' ) ; % Plot x-axis 

title( 'Zeros of f 3 ' ) 

hold off: 



Zeros of f3 



2- 




-2 




-<>- ~z¿*~ 



-O- i- -_^ 




Finding the Máxima and Minima of f2 

To find the máxima and minima of f"(x) , calcúlate the valué of f"(x) at each 
of the zeros oíf'"(x) . To do so, substitute zeros into f2 and display the result 
below zeros: 

[zerosd; subs(f2,zerosd) ] 
ans = 

2.4483 -2.4483 3.1416 
0.0494 1.0051 1.0051 -4.0000 
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This shows the following: 

• f"{x) has an absolute máximum at x = +2.4483 , whose valué is 1.0051. 



• f"(x) has an absolute mínimum at x = n , whose valué is -4. 

• f"{x) has a local minimum at x = n, whose valué is 0.0494. 

You can display the máxima and minima with the following commands: 

clf 

ezplot(f2) 

axis([-2*pi 2*pi -4.5 1 .5]) 

ylabel( ' f 2 ' ) ; 

title( 'Máxima and Minima of f 2 ' ) 

hold on 

plot(zeros, subs(f2,zeros) , ' ro ' ) 

text(-4, 1.25, 'Absolute máximum') 

text ( -1 , -0 .25, ' Local minimum') 

text(.9, 1.25, 'Absolute máximum') 

text(1.6, -4.25, 'Absolute minimum') 

hold off; 
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This displays the following figure. 



Máxima and Mínima of f2 



Absolute máximum Absolute máximum 




Absolute mínimum 



The preceding analysis shows that the actual range of f"(x) is [-4, 1.0051]. 

Integrating 

To see whether integrating f"(x) twice with respect to x recovers the original 
function f(x) = 1/(5 + 4cosx) , enter the command 

g = int(int(f2)) 
which returns 

9 = 
-8/(tan(1/2*x) A 2+9) 

This is certainly not the original exprés sion for f(x) . Now look at the difference 

f(x)-g(x). 

d = f - g 
pretty (d) 
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1 8 
- + 

5+4 cos(x) 2 

tan (1/2 x) + 9 

You can simplify this using simple (d) or simplif y(d). Either command 
produces 

ans = 
1 

This illustrates the concept that differentiating f{x) twice, then integrating the 
result twice, produces a function that may differ from f(x) by a linear function 
of x . 

Finally, intégrate f(x) once more: 

F = int(f) 

The result 

F = 
2/3*atan(1/3*tan(1/2*x)) 

involves the arctangent function. 

Note that F{x) is not an antiderivative of f{x) for all real numbers, since it is 
discontinuous at odd múltiples of n , where tanx is singular. You can see the 
gaps in F(x) in the following figure. 
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ezplot(F) 



o 

-0.2 
-0.4 



2/3atan(1/3tan(1/2x)) 



To change F(x) into a true antiderivative of f(x) that is differentiable 
everywhere, you can add a step function to F(x) . The height of the steps is the 
height of the gaps in the graph of F(x) . You can determine the height of the 
gaps by taking the limits of F(x) as x approaches n from the left and from the 
right. The limit from the left is 



limit(F, x, 
ans = 
1/3*pi 



Pi, 



'left' ) 



On the other hand, the limit from the right is 

limit(F, x, pi, 'right') 
ans =-1/3*pi 

The height of the gap is the distance between the left and right hand limits, 
which is 2k/3 , as shown in the following figure. 
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71/3 



Gapofheight 2 71/ 3 




-71/3 



You can créate the step function using the round function, which rounds 
numbers to the nearest integer, as follows: 

J = sym(2*pi/3)*sym( ' round(x/(2*pi) ) ' ) ; 

Each step has width 27t and the jump from one step to the next is 2n/3 , as 
shown in the following figure. 
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2/3 7iround(1/2x/ít) 
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Next, add the step function J(x) to F(x) with the following code: 

F1 = F+J 

F1 = 

2/3*atan(1/3*tan(1/2*x))+2/3*pi*round(1/2*x/pi) 

Adding the step function raises the section of the graph of F(x) on the interval 
[7i 3jt) up by 2k/3 , lowers the section on the interval (-37: -7t] down by 
2tc/3 , and so on, as shown in the following figure. 
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Moves down 
2tc/3 



Moves up 
27C/3 



When you plot the result by entering 
ezplot(F1 ) 

you see that this representation does have a continuous graph. 



2.5 : 

2 - 

1.5 - 

1 - 

0.5 

- 

-0.5 - 

-1 - 

-1.5 - 

-2 ; 
-2.5 : 



2/3 atan(1/3 tan(1/2 x))+2/3 n round(1/2 x/it) 
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Símplífícatíons and Substitutions 



There are several functions that simplify symbolic expressions and are used to 
perform symbolic substitutions: 

Símplífícatíons 

Here are three different symbolic expressions. 

syms x 

f = x~3-6*x~2+11*x-6 
g = (x-1)*(x-2)*(x-3) 
h = -6+(11+(-6+x)*x)*x 

Here are their prettyprinted forms, generated by 
pretty(f), pretty(g), pretty(h) 

3 2 

x - 6 x + 1 1 x - 6 

(x - 1) (x - 2) (x - 3) 

-6 + (11 + (-6 + x) x) x 

These expressions are three different representations of the same 
mathematical function, a cubic polynomial in x. 

Each of the three forms is preferable to the others in different situations. The 
first form, f , is the most commonly used representation of a polynomial. It is 
simply a linear combination of the powers of x. The second form, g, is the 
factored form. It displays the roots of the polynomial and is the most accurate 
for numerical evaluation near the roots. But, if a polynomial does not have such 
simple roots, its factored form may not be so convenient. The third form, h, is 
the Horner, or nested, representation. For numerical evaluation, it involves the 
fewest arithmetic operations and is the most accurate for some other ranges of 
x. 

The symbolic simplification problem involves the verification that these three 
expressions represent the same function. It also involves a less clearly defined 
objective — which of these representations is "the simplest"? 
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This toolbox provides several functions that apply various algebraic and 
trigonometric identities to transform one representation of a function into 
another, possibly simpler, representation. These functions are collect, 
expand, horner, factor, simplif y, and simple. 

collect 

The statement 

collect(f) 

views f as a polynomial in its symbolic variable, say x, and collects all the 
coefficients with the same power of x. A second argument can specify the 
variable in which to collect terms if there is more than one candidate. Here are 
a few examples. 



f 


eolleet(f) 


(x-1 )*(x-2)*(x-3) 


x"3-6*x"2+11*x-6 


x*(x*(x-6)+11 ) -6 


x~3-6*x~2+11*x-6 


(1+x)*t + x*t 


2*x*t+t 
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expand 

The statement 

expand(f ) 

distributes producís over sums and applies other identities involving functions 
of sums as shown in the examples below. 

f expand(f) 



a*(x + y 

(x-1)*(: 

x*(x*(x 

exp(a+b 

cos(x+y 



a*x + a*y 



(x-1 )*(x-2)*(x-3) 


x"3-6*x~2+11*x-6 


x*(x*(x-6)+11 ) -6 


x"3-6*x"2+11*x-6 


exp(a+b) 


exp(a)*exp(b) 



cos(3*acos(xl 



cos(x)*cos(y)-sin(x)*sin(y) 
4*x"3-3*x 



horner 

The statement 

horner(f ) 

transforms a symbolic polynomial f into its Horner, or nested, representation 
as shown in the following examples. 



f 


horner(f) 


x"3-6*x"2+11*x-6 


-6+(11+(-6+x)*x)*x 


1 .1+2.2*x+3.3*x"2 


11/10+(11/5+33/10*x)*x 
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factor 

If f is a polynomial with rational coefficients, the statement 

f actor(f ) 

exprés ses f as a product of polynomials of lower degree with rational 
coefficients. If f cannot be factored over the rational numbers, the result is f 
itself. Here are several examples. 



f 


factor(f) 


x"3-6*x"2+1 1*x-6 


(x-1 )*(x-2)*(x-3) 


x"3-6*x"2+1 1*x-5 


x"3-6*x"2+1 1*x-5 


x"6+1 


(x"2+1 )*(x"4-x"2+1 ) 



Here is another example involving factor. It factors polynomials of the form 
x " n + 1 . This code 

syms x; 
n = (1:9)'; 
p = x. "n + 1 ; 
f = factor(p) ; 
[P. f] 

returns a matrix with the polynomials in its first column and their factored 
forms in its second. 



x+1 






x+1 ] 


'2+1 






x"2+1 ] 


"3+1 






(x+1 )*(x"2-x+1 ) ] 


"4+1 






x"4+1 ] 


"5+1 






(x+1 )*(x"4-x"3+x"2-x+1 ) ] 


"6+1 






(x"2+1 )*(x"4-x"2+1 ) ] 


"7+1 


(x+1 


)*(1 


-x+x"2-x"3+x"4-x"5+x"6) ] 


"8+1 






x"8+1 ] 


"9+1 




(x+1 


)*(x"2-x+1 )*(x"6-x"3+1 ) ] 
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As an aside at this point, factor can also factor symbolic objects containing 
integers. This is an alternative to using the factor function in the MATLAB 
specf un directory. For example, the following code segment 

N = sym(1 ) ; 
for k = 2:11 

N(k) = 10*N(k-1)+1 ; 
end 
[N 1 factor(N')] 

displays the factors of symbolic integers consisting of ls: 



1 

11 

111 

1111 

11111 



1 ] 

(11)] 


(3)*(37)] 


(1 1 )*(101 ) ] 


(41)*(271)] 


3)*(7)*(11 )*(13)*(37)] 


(239)*(4649)] 


(11 )*(73)*(101 )*(137)] 


(3)~2*(37)*(333667)] 


11)* (41 )*(271)*(9091)] 


(513239)*(21649)] 
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simplify 

The simplify function is a powerful, general purpose tool that applies a 
number of algebraic identities involving sums, integral powers, square roots 
and other fractional powers, as well as a number of functional identities 
involving trig functions, exponential and log functions, Bessel functions, 
hypergeometric functions, and the gamma function. Here are some examples. 



f 


simplify(f) 




x*(x*(x-6)+1 1 ) -6 


x"3-6*x"2+1 1*x-6 




(1-x"2)/(1-x) 


x+1 





;i/a*3+6/a*2+12/a+8) A (1/3) 



((2*a+1)*3/a*3)*(1/3) 



syms x y positive 
log(x*y) 



exp(x) * exp(y) 
besselj(2,x) + besselj (0,x) 
gamma(x+1 ) -x*gamma(x) 
cos(x)"2 + sin(x)"2 



log(x)+log(y) 
exp(x+y) 

2/x*besselj (1 ,x) 




1 



simple 

The simple function has the unorthodox mathematical goal of finding a 
simplification of an expression that has the fewest number of characters. Of 
course, there is little mathematical justification for claiming that one 
expression is "simpler" than another just because its ASCII representation is 
shorter, but this often proves satisfactory in practice. 

The simple function achieves its goal by independently applying simplify, 
collect, factor, and other simplification functions to an expression and 
keeping track of the lengths of the results. The simple function then returns 
the shortest result. 

The simple function has several forms, each returning different output. The 
form 

simple(f ) 
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displays each trial simplification and the simplification function that produced 
it in the MATLAB Command Window. The simple function then returns the 
shortest result. For example, the command 

simple(cos(x) "2 + sin(x)~2) 

displays the following alternative simplifications in the MATLAB Command 
Window: 

simplify : 
1 

radsimp: 

cos(x) '2+sin(x) "2 

combine (trig) : 
1 

factor: 
cos(x)"2+sin(x) "2 

expand: 

cos(x) ~2+sin(x) "2 

combine: 

1 

convert (exp) : 
(1/2*exp(i*x)+1/2/exp(i*x)) A 2-1/4*(exp(i*x) -1/exp(i*x) ) "2 

convert (sincos) : 
cos(x) ~2+sin(x) "2 

convert (tan) : 

(1 -tan(1/2*xK2p2/(1+tan(1/2*xK2K2+ 

4*tan(1/2*x)~2/(1+tan(1/2*x)~2K2 

collect (x) : 
cos(x) "2+sin(x) "2 
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and returns 

ans = 
1 

This form is useful when you want to check, for example, whether the shortest 
form is indeed the simplest. If you are not interested in how simple achieves 
its result, use the form 

f = simple(f) 

This form simply returns the shortest expression found. For example, the 
statement 

f = simple(cos(x) ~2+sin(x) "2) 

returns 

f = 

1 

If you want to know which simplification returned the shortest result, use the 
múltiple output form: 

[F, how] = simple(f) 

This form returns the shortest result in the first variable and the simplification 
method used to achieve the result in the second variable. For example, the 
statement 

[f, how] = simple(cos(x) ~2+sin(x) ~2) 

returns 

f = 

1 

how = 
combine 

The simple function sometimes improves on the result returned by simplif y, 
one of the simplifications that it tries. For example, when applied to the 
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examples given for simplif y, simple returns a simpler (or at least shorter) 
result in two cases. 



simplify(f) 



simple(f) 



;i/a A 3+6/a"2+12/a+8P(1/3) ( (2*a+1 ) ~3/a~3) A (1 /3) 



(2*a+1 )/a 



syms x y positive 
log(x*y) 



log(x)+log(y) 



log(x*y) 



In some cases, it is advantageous to apply simple twice to obtain the effect of 
two different simplification functions. For example, the statements 

f = (1/a A 3+6/a"2+12/a+8K(1/3) ; 
simple(simple(f ) ) 

return 

2+1 /a 

The first application, simple (f), uses radsimp to produce (2*a+1 ) /a; the 
second application uses combine (trig) to transform this to 1 /a+2. 

The simple function is particularly effective on expressions involving 
trigonometric functions. Here are some examples. 



f 


simple(f) 




cos(x) ~2+sin(x) "2 


1 


2*cos(x)"2-sin(x)"2 


3*cos(x)"2-1 


cos(x) ~2-sin(x) '2 


cos(2*x) 


cos(x) + (-sin(x)-2P(1/2) 


cos(x)+i*sin(x) 


cos(x)+i*sin(x) 


exp(i*x) 




cos(3*acos(x) ) 


4*x"3-3*x 
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Substítutíons 

There are two functions for symbolic substitution: subexpr and subs. 

subexpr 

These commands 

syms a x 

s = solve(x*3+a*x+1 ) 

solve the equation x"3+a*x+1 = for x: 

s = 

[ 1/6*(-108+12*(12*a*3+81K(1/2)K(1/3) -2*a/ 

(-108+12*(12*a~3+81 ) " (1 /2) ) " (1 /3) ] 
[ -1/12*(-108+12*(12*a"3+81 ) " (1 /2) ) ' (1 /3)+a/ 

(-108+12*(12*a"3+81 ) ~ (1 /2) ) ~ (1 /3)+1 /2*i*3~ (1 /2)*(1 / 

6*(-108+12*(12*a~3+81 ) " (1 /2) ) ~ (1 /3)+2*a/ 

(-108+12*(12*a~3+81 ) * (1 /2) ) " (1 /3) ) ] 
[ -1/12*(-108+12*(12*a"3+81 ) " (1 /2) ) ' (1 /3)+a/ 

(-108+12*(12*a"3+81 ) * (1 /2) ) ~ (1 /3) -1 /2*i*3~ (1 /2)*(1 / 

6*(-108+12*(12*a~3+81 ) " (1 /2) ) ~ (1 /3)+2*a/ 

(-108+12*(12*a"3+81 ) * (1 /2) ) " (1 /3) ) ] 
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Use the pretty function to display s in a more readable form: 

pretty(s) 



1/3 a 

1/6 %1 - 2 

1/3 
%1 



] 



1/3 a 1/2 / 1/3 a V 

1/12 %1 + + 1/2 i 3 |1/6%1 +2 |; 

1/3 | 1/3|; 

%1 \ %1 /; 



1/3 
1/12 %1 + 



1/3 



1/2 / 1/3 

1/2 i 3 11/6 %1 +2 



\: 



%i 



\ 



i/3i: 

%i /: 



3 1/2 

%1 := -108 + 12 (12 a + 81 ) 

The pretty command inherits the %n (n, an integer) notation from Maple to 
denote subexpressions that occur múltiple times in the symbolic object. The 
subexpr function allows you to save these common subexpressions as well as 
the symbolic object rewritten in terms of the subexpressions. The 
subexpressions are saved in a column vector called sigma. 

Continuing with the example 

r = subexpr(s) 

returns 

sigma = 

-108+12*(12*a"3+81)"(1/2) 
r = 

[ 1/6*sigma"(1/3) -2*a/ sigma" (1 /3) ; 

[ -1/12*sigma"(1/3)+a/ sigma" (1/3)+1/2*i*3"(1/2)*(1/6*sigma~ 
(1/3)+2*a/sigma~(1/3))] 
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[ -1 /12*sigma'(1 /3)+a/sigma~ (1 /3) -1 /2*i*3~(1 /2)*(1 /6*sigma~ 
(1/3)+2*a/sigma~(1/3))] 

Notice that subexpr creates the variable sigma in the MATLAB workspace. 
You can verify this by typing whos , or the command 

sigma 

which returns 

sigma = 
-108+12*(12*a"3+81 P(1/2) 

subs 

The following code finds the eigenvalues and eigenvectors of a circulant matrix 
A: 

syms abe 

A= [abe; bea; cab]; 

[v,E] = eig(A) 



[ -(a+(b~2-b*a-c*b-c*a+a~2+c~2K(1/2)-b)/(a-c) , 

-(a-(b~2-b*a-c*b-c*a+a~2+c~2P(1/2)-b)/(a-c) , 1 ] 

[ -(b-c-(br2-b*a-c*b-c*a+a A 2+c*2) A (1/2))/(a-c) , 

-(b-c+(b~2-b*a-c*b-c*a+a~2+c~2P(1/2) )/(a-c) , 1 ] 

[ 1, 

1, 1] 



[ (b~2-b*a-c*b- 

c*a+a~2+c~2K(1/2) , 0, 0] 

[ 0, -(b~2-b*a-c*b- 

c*a+a~2+c~2)"(1/2) , 0] 

[ 0, 0, b+c+a] 
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Note MATLAB might return the eigenvalues that appear on the diagonal of 
E in a different order. In this case, the corresponding eigenvectors, which are 
the columns of v, will also appear in a different order. 



Suppose you want to replace the rather lengthy exprés sion 

(tT2-b*a-c*b-c*a+a~2+c~2P(1/2) 

throughout v and E. First, use subexpr 
v = subexpr(v, ' S' ) 

which returns 

S = 
(b~2-b*a-c*b-c*a+a~2+c~2)~(1/2) 

v = 

[ -(a+S-b)/(a-c), - (a-S-b) / (a-c) , 1] 

[ -(b-c-S)/(a-c), -(b-c+S)/(a-c), 1] 

[ 1, 1, 1] 

Next, substitute the symbol S into E with 
E = subs(E,S, 'S' ) 

E = 

[ S, 0, 0] 

[ 0, -S, 0] 

[ 0, 0, b+c+a] 

Now suppose you want to evalúate v at a = 10. You can do this using the subs 
command: 

subs(v,a, 10) 

This replaces all occurrences of a in v with 10. 

[ -(10+S-b)/(10-c), -(10-S-b)/(10-c) , 1] 

[ -(b-c-S)/(10-c) , -(b-c+S)/(10-c) , 1] 

[ 1, 1, 1] 
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Notice, however, that the symbolic expression that S represents is unaffected 
by this substitution. That is, the symbol a in S is not replaced by 10. The subs 
command is also a useful function for substituting in a variety of valúes for 
several variables in a particular expression. For example, suppose that in 
addition to substituting a = 1 in S, you also want to substitute the valúes for 
2 and 10 for b and c, respectively. The way to do this is to set valúes for a, b, 
and c in the workspace. Then subs evaluates its input using the existing 
symbolic and double variables in the current workspace. In the example, you 
first set 

a = 10; b = 2; c = 10; 
subs(S) 

ans = 
8 



To look at the contents of the workspace, type whos, which gives 
Ñame Size Bytes Class 



A 


3x3 


878 


sym object 


E 


3x3 


888 


sym object 


S 


1x1 


186 


sym object 


a 


1x1 


8 


double array 


ans 


1x1 


140 


sym object 


b 


1x1 


8 


double array 


c 


1x1 


8 


double array 


V 


3x3 


982 


sym object 



a, b, and c are now variables of class double while A, E, S, and v remain symbolic 
expressions (class sym). 

If you want to preserve a, b, and c as symbolic variables, but still alter their 
valué within S, use this procedure. 

syms abe 
subs(S,{a,b,c}, {10,2,10}) 

ans = 
8 

Typing whos reveáis that a, b, and c remain 1-by-l sym objeets. 
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The subs command can be combined with double to evalúate a symbolic 
expression numerically. Suppose you have the following expressions 

syms t 

M = (1 -t"2)*exp(-1/2*-r2) ; 

P = (1 -t"2)*sech(t) ; 

and want to see how M and P differ graphically. 
One approach is to type 

ezplot(M); hold on; ezplot(P); hold off; 

but this plot does not readily help us identify the curves. 

(1-t 2 ) sech(t) 




10-55 



I w Using the Symbolic Math Toolbox 



Instead, combine subs, double, and plot 

T = -6:0.05:6; 

MT = double ( subs (M,t,T) ) ; 

PT = double(subs(P,t,T) ) ; 

plot(T,MT, 'b' ,T,PT, r-. ' ) 

titleC ') 

legend( ' M 1 , 'P' ) 

xlabel( ' t ' ) ; grid 

to produce a multicolored graph that indicates the difference between M and P. 
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Finally the use of subs with strings greatly facilitates the solution of problems 
involving the Fourier, Laplace, orz-transforms. 
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Variable-Precisión Arithmetic 



Overview 

There are three different kinds of arithmetic operations in this toolbox: 
Numeric MATLAB floating-point arithmetic 

Rational Maple's exact symbolic arithmetic 

VPA Maple's variable-precisión arithmetic 

For example, the MATLAB statements 

format long 
1/2+1/3 

use numeric computation to produce 
0.83333333333333 

With the Symbolic Math Toolbox, the statement 
sym(1/2)+1/3 

uses symbolic computation to yield 
5/6 

And, also with the toolbox, the statements 

digits(25) 
vpa( '1/2+1/3' ) 

use variable-precisión arithmetic to return 
. 8333333333333333333333333 

The floating-point operations used by numeric arithmetic are the fastest of the 
three, and require the least computer memory, but the results are not exact. 
The number of digits in the printed output of MATLAB double quantities is 
controlled by the format statement, but the internal representation is always 
the eight-byte floating-point representation provided by the particular 
computer hardware. 
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In the computation of the numeric result above, there are actually three 
roundoff error s, one in the división of 1 by 3, one in the addition of 1/2 to the 
result of the división, and one in the binary to decimal conversión for the 
printed output. On computers that use IEEE floating-point standard 
arithmetic, the resulting internal valué is the binary expansión of 5/6, 
truncated to 53 bits. This is approximately 16 decimal digits. But, in this 
particular case, the printed output shows only 15 digits. 

The symbolic operations used by rational arithmetic are potentially the most 
expensive of the three, in terms of both computer time and memory. The results 
are exact, as long as enough time and memory are available to complete the 
computations. 

Variable-precisión arithmetic falls in between the other two in terms of both 
cost and accuracy. A global parameter, set by the function digits, controls the 
number of significant decimal digits. Increasing the number of digits increases 
the accuracy, but also increases both the time and memory requirements. The 
default valué of digits is 32, corresponding roughly to floating-point accuracy. 

The Maple documentation uses the term "hardware floating-point" for what 
you are calling "numeric" or "floating-point" and uses the term "floating-point 
arithmetic" for what you are calling "variable-precisión arithmetic." 



Example: Using the Dífferent Kínds of Arithmetic 

Rational Arithmetic 

By default, the Symbolic Math Toolbox uses rational arithmetic operations, i.e., 
Maple's exact symbolic arithmetic. Rational arithmetic is invoked when you 
créate symbolic variables using the sym function. 

The sym function converts a double matrix to its symbolic form. For example, if 
the double matrix is 



A 


= 






1 . 


,1000 


1 .2000 


1 .3000 


2. 


,1000 


2.2000 


2.3000 


3. 


,1000 


3.2000 


3.3000 
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its symbolic form, S = sym(A), is 

S = 

[11/10, 6/5, 13/10] 

[21/10, 11/5, 23/10] 

[31/10, 16/5, 33/10] 

For this matrix A, it is possible to discover that the elements are the ratios of 
small integers, so the symbolic representation is formed from those integers. 
On the other hand, the statement 

E = [exp(1) (1+sqrt(5))/2; log(3) rand] 

returns a matrix 

E = 

2.71828182845905 1 .61803398874989 
1 .09861228866811 0.76209683302739 

whose elements are not the ratios of small integers, so sym(E) reproduces the 
floating-point representation in a symbolic form: 

ans = 

[ 6121026514868074*2" ( -51 ) , 7286977268806824*2" ( -52) ] 

[ 4947709893870346*2" ( -52) , 6864358026484820*2" ( -53) ] 

Variable-Precisión Numbers 

Variable-precisión numbers are distinguished from the exact rational 
representation by the presence of a decimal point. A power of 10 scale factor, 
denoted by ' e ' , is allowed. To use variable-precisión instead of rational 
arithmetic, créate your variables using the vpa function. 

For matrices with purely double entries, the vpa function generates the 
representation that is used with variable-precisión arithmetic. For example, if 
you apply vpa to the matrix S defined in the preceding section, with digits (4 ) , 
by entering 

vpa(S) 
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MATLAB returns the output 

S = 

[1 .100, 1 .200, 1 .300] 

[2.100, 2.200, 2.300] 

[3.100, 3.200, 3.300] 

Applying vpa to the matrix E defined in the preceding section, with digit s ( 25 ) , 
by entering 

digits(25) 
F = vpa(E) 

returns 

F = 

[2.718281828459045534884808, 1 .414213562373094923430017] 

[1 .098612288668110004152823, .2189591863280899719512718] 

Converting to Floating-Point 

To convert a rational or variable-precisión number to its MATLAB 
floating-point representation, use the double function. 

In the example, both double(sym(E) ) and double (vpa (E) ) return E. 

Another Example 

The next example is perhaps more interesting. Start with the symbolic 
exprés sion 

f = sym( 'exp(pi*sqrt(163) ) ' ) 

The statement 
double(f ) 

produces the printed floating-point valué 

ans = 

2 . 6253741 26407687e+01 7 
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Using the second argument of vpa to specify the number of digits, 

vpa(f,18) 

returns 

262537412640768744. 

whereas 

vpa(f,25) 

returns 

262537412640768744.0000000 

you suspect that f might actually have an integer valué. This suspicion is 
reinforced by the 30 digit valué, vpa (f , 30 ) 

26253741 2640768743 . 999999999999 

Finally, the 40 digit valué, vpa (f ,40) 

26253741 2640768743 . 9999999999992500725944 

shows that f is very cióse to, but not exactly equal to, an integer. 
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Basic Algebraic Operations 

Basic algebraic operations on symbolic objects are the same as operations on 
MATLAB objects of class double. This is illustrated in the following example. 

The Givens transformation produces a plañe rotation through the angle t. The 
statements 

syms t; 

G = [cos(t) sin(t); -sin(t) cos(t)] 

créate this transformation matrix. 

G = 

[ cos(t) , sin(t) ] 
[ -sin(t) , cos(t) ] 

Applying the Givens transformation twice should simply be a rotation through 
twice the angle. The corresponding matrix can be computed by multiplying G 
by itself or by raising G to the second power. Both 

A = G*G 

and 

A = G~2 

produce 

A = 

[cos(t)"2-sin(t)"2, 2*cos(t)*sin(t)] 
[ -2*cos(t)*sin(t) , cos(t) A 2-sin(t) A 2] 

The simple function 
A = simple (A) 

uses a trigonometric identity to return the expected form by trying several 
different identities and picking the one that produces the shortest 
representation. 
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A = 

[ cos(2*t), sin(2*t)] 

[-sin(2*t), cos(2*t)] 

The Givens rotation is an orthogonal matrix, so its transpose is its inverse. 
Confirming this by 

I = G. ' *G 



which produces 




I = 
[cos(t)"2+sin(tp2, 

[ o, 


0] 
cos(t)"2+sin(t)"2] 


and then 




I = simple (I) 
I = 
[1, 0] 
[0, 1] 





Linear Algebraíc Operatíons 

The following examples show to do several basic linear algebraic operations 
using the Symbolic Math Toolbox. 

The command 

H = hilb(3) 

generates the 3-by-3 Hilbert matrix. With f ormat short, MATLAB prints 



H 


= 










1 


.0000 





.5000 





.3333 





.5000 





.3333 





.2500 





.3333 





.2500 





.2000 



The computed elements of H are floating-point numbers that are the ratios of 
small integers. Indeed, H is a MATLAB array of class double. Converting H to 
a symbolic matrix 

H = sym(H) 
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gives 

[ 1, 1/2, 1/3] 
[1/2, 1/3, 1/4] 
[1/3, 1/4, 1/5] 

This allows subsequent symbolic operations on H to produce results that 
correspond to the infinitely precise Hilbert matrix, sym(hilb(3) ), not its 
floating-point approximation, hilb(3). Therefore, 

inv(H) 

produces 

[ 9, -36, 30] 
[-36, 192, -180] 
[ 30, -180, 180] 

and 

det(H) 

yields 
1/2160 

You can use the backslash operator to solve a system of simultaneous linear 
equations. For example, the commands 

b = [1 1 1 ] ' 

x = H\b % Solve Hx = b 

produce the solution 

[ 3] 
[-24] 

[ 30] 

All three of these results, the inverse, the determinant, and the solution to the 
linear system, are the exact results corresponding to the infinitely precise, 
rational, Hilbert matrix. On the other hand, using digits ( 16), the command 

V = vpa(hilb(3)) 
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returns 

[ 1., .5000000000000000, .3333333333333333] 

[.5000000000000000, .3333333333333333, .2500000000000000] 
[.3333333333333333, .2500000000000000, .2000000000000000] 

The decimal points in the representation of the individual elements are the 
signal to use variable-precisión arithmetic. The result of each arithmetic 
operation is rounded to 16 significant decimal digits. When inverting the 
matrix, these errors are magnified by the matrix condition number, which for 
hilb(3) is about 500. Consequently, 

inv(V) 

which returns 

ans = 

[ 9.000000000000179, -36.00000000000080, 30.00000000000067] 

[ -36.00000000000080, 192.0000000000042, -180.0000000000040] 

[ 30.00000000000067, -180.0000000000040, 180.0000000000038] 

shows the loss of two digits. So does 

det(V) 

which gives 

.462962962962953e-3 

and 
V\b 

which is 

[ 3.000000000000041] 
[-24.00000000000021 ] 
[ 30.00000000000019] 

Since H is nonsingular, calculating the nuil space of H with the command 
null(H) 

returns an empty matrix, and calculating the column space of H with 
colspace(H) 
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returns a permutation of the identity matrix. A more interesting example, 
which the following code shows, is to find a valué s for H ( 1 , 1 ) that makes H 
singular. The commands 

syms s 
H(1,1) = s 
Z = det(H) 
sol = solve(Z) 

produce 

H = 

[ s, 1/2, 1/3] 

[1/2, 1/3, 1/4] 

[1/3, 1/4, 1/5] 

Z = 
1/240*s-1/270 

sol = 
8/9 

Then 

H = subs(H,s,sol) 

substitutes the computed valué of sol for s in H to give 

H = 

[8/9, 1/2, 1/3] 

[1/2, 1/3, 1/4] 

[1/3, 1/4, 1/5] 

Now, the command 
det(H) 

returns 

ans = 



and 

inv(H) 
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produces an error message 

??? error using ==> inv 

Error, (in inverse) singular matrix 

because H is singular. For this matrix, Z = null(H) and C = colspace(H) are 
nontrivial: 

Z = 
[ 1] 
[ -4] 
[10/3] 

C = 

[ 1, 0] 

[ 0, 1] 

[ -3/10, 6/5] 

It should be pointed out that even though H is singular, vpa ( H ) is not. For any 
integer valué d, setting 

digits(d) 

and then computing 

det(vpa(H)) 

inv(vpa(H) ) 

results in a determinant of size 1 0" ( -d ) and an inverse with elements on the 
order of 10~d. 

Eígenvalues 

The symbolic eigenvalues of a square matrix A or the symbolic eigenvalues and 
eigenvectors of A are computed, respectively, using the commands 

E = eig(A) 
[V,E] = eig(A) 

The variable-precisión counterparts are 

E = eig(vpa(A)) 
[V,E] = eig(vpa(A)) 
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The eigenvalues of A are the zeros of the characteristic polynomial of A, 
det (A-x*I ), which is computed by 

poly(A) 

The matrix H from the last section provides the first example: 

H = 

[8/9, 1/2, 1/3] 

[1/2, 1/3, 1/4] 

[1/3, 1/4, 1/5] 

The matrix is singular, so one of its eigenvalues must be zero. The statement 
[T,E] = eig(H) 

produces the matrices T and E. The columns of T are the eigenvectors of H: 

T = 

[ 1, 28/153+2/153*12589"(1/2) , 28/153-2/153*12589^(12)] 

[ "4, 1, 1] 

[ 10/3, 92/255-1/255*12589"(1/2) , 292/255+1/255*12589^(12)] 

Similarly, the diagonal elements of E are the eigenvalues of H: 

E = 

[0, 0, 0] 

[0, 32/45+1/180*12589^(1/2), 0] 

[0, 0, 32/45-1/180*12589"(1/2) ] 

It may be easier to understand the structure of the matrices of eigenvectors, T, 
and eigenvalues, E, if you convert T and E to decimal notation. To do so, proceed 
as follows. The commands 



Td 


= double(T) 




Ed 


= double(E) 




-un. 
Td 


L 






1.0000 1.6497 


-1 .2837 




-4.0000 1.0000 


1 .0000 




3.3333 0.7051 


1 .5851 
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Ed 


= 






















1 .3344 













0.0878 



The first eigenvalue is zero. The corresponding eigenvector (the first column of 
Td) is the same as the basis for the nuil space found in the last section. The 
other two eigenvalues are the result of applying the quadratic formula to 

x~2-64/45*x+253/2160 

which is the quadratic factor in f actor (poly(H) ). 

syms x 

g = simple(f actor(poly(H) ) /x) ; 

solve(g) 

ans = 

[ 32/45+1/180*12589"(1/2) ] 

[ 32/45-1/180*12589^(1/2)] 

Closed form symbolic expressions for the eigenvalues are possible only when 
the characteristic polynomial can be exprés sed as a product of rational 
polynomials of degree four or less. The Rosser matrix is a classic numerical 
analysis test matrix that illustrates this requirement. The statement 



R = sym(gallery( ' rosser ' 


)) 








generates 
















R = 
[ 611 


196 


-192 


407 


-8 


-52 


-49 


29 


[ 196 


899 


113 


-192 


-71 


-43 


-8 


-44 


[-192 


113 


899 


196 


61 


49 


8 


52 


[ 407 


-192 


196 


611 


8 


44 


59 


-23 


[ -8 


-71 


61 


8 


411 


-599 


208 


208 


[ -52 


-43 


49 


44 


-599 


411 


208 


208 


[ -49 


-8 


8 


59 


208 


208 


99 


-911 


[ 29 


-44 


52 


-23 


208 


208 


-911 


99 


The commands 














p = po 


iy(R); 














pretty (factor 


(P)) 
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produce 



1020) (x 



1020 x + 100) 



1040500) 



1000) 



The characteristic polynomial (of degree 8) factors nicely into the product of two 
linear terms and three quadratic terms. You can see immediately that four of 
the eigenvalues are 0, 1020, and a double root at 1000. The other four roots are 
obtained from the remaining quadratics. Use 

eig(R) 
to find all these valúes 



ans 



10*10405' 

-10*10405' 

510+100*26' 

510-100*26' 



0] 

1020] 

:i/2)] 

:i/2)] 
:i/2)] 
:i/2)] 

1000] 
1000] 



The Rosser matrix is not a typical example; it is rare for a full 8-by-8 matrix to 
have a characteristic polynomial that factors into such simple form. If you 
change the two "córner" elements of R from 29 to 30 with the commands 

S = R; S(1 ,8) = 30; S(8,1 ) = 30; 

and then try 
p = poly(S) 

you find 

P = 

x"8-4040*x"7+5079941 *x"6+82706090*x"5 - 5327831 91 8568*x~4+ 
428783291 271 9760*x~3 - 1 08269938841 1 1 66000*x~2+51 264008540948000*x 
+40250968213600000 

You also find that f actor(p) is p itself. That is, the characteristic polynomial 
cannot be factored over the rationals. 
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For this modified Ros ser matrix 

F = eig(S) 

returns 

F = 

.21803980548301606860857564424981 ] 
999 . 94691 786044276755320289228602 ] 
1 000 . 1 206982933841 33571 281 7075454 ] 
1019. 524355263201 6358324933278291 ] 
1019. 9935501 291 629257348091 8081 73 ] 
1 020 . 4201 88201 50472781 85457498840 ] 

■ .17053529728768998575200874607757] 

■ 1 020 . 05321 4255891 51 65931 894252600 ] 

Notice that tríese valúes are cióse to the eigenvalues of the original Rosser 
matrix. Further, the numerical valúes of F are a result of Maple's floating-point 
arithmetic. Consequently, different settings of digits do not alter the number 
of digits to the right of the decimal place. 

It is also possible to try to compute eigenvalues of symbolic matrices, but closed 
form solutions are rare. The Givens transformation is generated as the matrix 
exponential of the elementary matrix 

1 
-1 



The Symbolic Math Toolbox commands 

syms t 

A = sym([0 1 ; -1 0]) ; 

G = expm(t*A) 

return 

[ cos(t) , sin(t) ] 
[ -sin(t) , cos(t) ] 

Next, the command 
9 = eig(G) 
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produces 

9 = 

[ cos(t) + (cos(tK2-ip(1/2)] 

[ cos(t)-(cos(t)"2-1)-(1/2)] 

you can use simple to simplify this form of g. Indeed, a repeated application of 

simple 

for j =1:4 

[g,how] = simple(g) 
end 

produces the best result: 



[ cos(t) + (-sin(tp2r(1/2) 
[ cos(t)-(-sin(tr2P(1/2) 

how = 
mwcos2sin 

g = 

[ cos(t)+i*sin(t) ] 
[ cos(t) -i*sin(t) ] 

how = 
radsimp 

g = 

[ exp(i*t)] 
[ 1/exp(i*t)] 

how = 
convert (exp) 

g = 

[ exp(i*t)] 
[ exp(-i*t)] 

how = 
simplify 
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Notice the first application of simple uses mwcos2sin to produce a sum of sines 
and cosines. Next, simple invokes radsimp to produce eos (t) + i*sin(t) for 
the first eigenvector. The third application of simple uses convert (exp) to 
change the sines and cosines to complex exponentials. The last application of 
simple uses simplif y to obtain the final form. 

Jordán Canonical Form 

The Jordán canonical form results from attempts to diagonalize a matrix by a 
similarity transformation. For a given matrix A, find a nonsingular matrix V, 
so that inv(V) *A*V, or, more succinctly, J = V\A*V, is "as cióse to diagonal as 
possible." For almost all matrices, the Jordán canonical form is the diagonal 
matrix of eigenvalues and the columns of the transformation matrix are the 
eigenvectors. This always happens if the matrix is symmetric or if it has 
distinct eigenvalues. Some nonsymmetric matrices with múltiple eigenvalues 
cannot be diagonalized. The Jordán form has the eigenvalues on its diagonal, 
but some of the superdiagonal elements are one, instead of zero. The statement 

J = jordan(A) 

computes the Jordán canonical form of A. The statement 
[V,J] = jordan(A) 

also computes the similarity transformation. The columns of V are the 
generalized eigenvectors of A. 

The Jordán form is extremely sensitive to perturbations. Almost any change in 
A causes its Jordán form to be diagonal. This makes it very difficult to compute 
the Jordán form reliably with floating-point arithmetic. It also implies that A 
must be known exactly (i.e., without round-off error, etc.). Its elements must be 
integers, or ratios of small integers. In particular, the variable-precisión 
calculation, jordan(vpa(A) ), is not allowed. 

For example, let 

A = sym([12,32,66,116;-25,-76,-164,-294; 
21 ,66, 143, 256; -6, -19, -41 , -73]) 



12, 


32, 


66, 


116] 


-25, 


-76, 


-164, 


-294] 


21, 


66, 


143, 


256] 


-6, 


-19, 


-41, 


-73] 
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Then 






[V,J] = 


jordan(A) 




produces 






V = 






[ 4, 


-2, 4, 


3 


[ "6, 


8, -11, 


-8 


[ 4, 


-7, 10, 


7 


[ -1, 


2, -3, 


-2 


J = 






[ 1, 1, 


0, 0] 




[ 0, 1, 


0, 0] 




[ o, 0, 


2, 1] 




[ o, 0, 


0, 2] 





Therefore A has a double eigenvalue at 1, with a single Jordán block, and a 
double eigenvalue at 2, also with a single Jordán block. The matrix has only 
two eigenvectors, V ( : , 1 ) and V ( : , 3 ) . They satisfy 

A*V(:,1) = 1*V(:,1) 
A*V( : ,3) = 2*V( : ,3) 

The other two columns of V are generalized eigenvectors of grade 2. They 
satisfy 

A*V(:,2) = 1*V(:,2) + V(:,1) 
A*V( : ,4) = 2*V( : ,4) + V(: ,3) 

In mathematical notation, with v¡ = v ( : , j ) , the columns of V and eigenvalues 
satisfy the relationships 

(A - X 1 I)v 2 = v 1 
{A-X 2 I)v A = v 3 

Singular Valué Decomposítíon 

Only the variable-precisión numeric computation of the complete singular 
vector decomposition is available in the toolbox. One reason for this is that the 
formulas that result from symbolic computation are usually too long and 
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complicated to be of much use. If A is a symbolic matrix of floating-point or 
variable-precisión numbers, then 

S = svd(A) 

computes the singular valúes of A to an accuracy determined by the current 
setting of digits. And 

[U,S,V] = svd(A); 

produces two orthogonal matrices, U and V, and a diagonal matrix, S, so that 
A = U*S*V ; 

Consider the n-by-n matrix A with elements defined by 
A(i,j) = 1/(i-j+1/2) 

For n = 5, the matrix is 

[ 2 -2 -2/3 -2/5 -2/7] 

[2/3 2 -2 -2/3 -2/5] 

[2/5 2/3 2 -2 -2/3] 

[2/7 2/5 2/3 2 -2] 

[2/9 2/7 2/5 2/3 2] 

It turns out many of the singular valúes of these matrices are cióse to n . 

The most obvious way of generating this matrix is 

fon i=1 :n 

for j =1 : n 

A(i,j) = sym(1/(i-i+1/2)); 
end 
end 

The most efficient way to genérate the matrix is 

[J , I] = meshgrid(1 :n) ; 
A = syra(1./(I - J+1/2)); 

Since the elements of A are the ratios of small integers, vpa(A) produces a 
variable-precisión representation, which is accurate to digits precisión. Henee 

S = svd(vpa(A) ) 
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computes the desired singular valúes to full accuracy. With n 
digits(30), the result is 



16 and 



1 .20968137605668985332455685357 
2 . 691 621 58686066606774782763594 
3 . 07790297231 1 1 9748658424727354 
3 . 1 3504054399744654843898901 261 
3.141 0604466347006380521 8371 924 
3.141 5575435991 8083691 050658260 
3.141 590754586058487289825771 1 9 
3.141 592569254923064702848631 02 
3.141 5926505265488081 556947961 3 
3.141 59265349961 0531 43856838564 
3.141 59265358767361 71 239261 2384 
3.141 59265358975439206849907220 
3.141 59265358979270342635559051 
3.141 592653589793233252901 42781 
3.141 5926535897932384306684671 2 
3.141 59265358979323846255035974 



There are two ways to compare S with pi, the floating-point representation of 
n . In the vector below, the first element is computed by subtraction with 
variable-precisión arithmetic and then converted toadouble. The second 
element is computed with floating-point arithmetic: 



format short e 
[double(pi*ones(16,1 ) -S) 



pi-double(S) ] 



10-76 



Linear Algebra 



The results are 






1 .9319e+00 


1 


,9319e+00 


4.4997e-01 


4 


,4997e-01 


6.3690e-02 


6 


.3690e-02 


6.5521e-03 


6 


.5521e-03 


5.3221e-04 


5 


.3221e-04 


3.5110e-05 


3 


.5110e-05 


1 .8990e-06 


1 


.8990e-06 


8.4335e-08 


8 


,4335e-08 


3.0632e-09 


3 


.0632e-09 


9.0183e-11 


9 


.0183e-11 


2.1196e-12 


2 


.1196e-12 


3.8846e-14 


3 


,8636e-14 


5.3504e-16 


4 


.4409e-16 


5.2097e-18 







3.1975e-20 







9.3024e-23 








Since the relativeaccuracyofpiis pi*eps, which is6.9757e-16, either column 
confirms the suspicion that four of the singular valúes of the 1 6-by-1 6 example 
equal n to floating-point accuracy. 

Eígenvalue Trajectoríes 

This example applies several numeric, symbolic, and graphic techniques to 
study the behavior of matrix eigenvalues as a parameter in the matrix is 
varied. This particular setting involves numerical analysis and perturbation 
theory, but the techniques illustrated are more widely applicable. 

In this example, you consider a 3-by-3 matrix A whose eigenvalues are 1, 2, 3. 
First, you perturb A by another matrix E and parameter t: A -> A + tE . As t 
increases from to 10" 6 , the eigenvalues X l = 1 , X 2 = 2 , X 3 = 3 change to 
V - 1.5596 + 0.2726¿ , X 2 ' - 1.5596 - 0.2726¿ , X 3 ' - 2.8808 . 



10-77 



I w Using the Symbolic Math Toolbox 





XV) 


































Mi) 


1(2) 
Ó 


X(3) 






fc(3) 


















- 






X 

V(2) 






- 



This, in turn, means that for some valué of t = x, < x < 10" 6 , the perturbed 
matrix A(t) = A + tE has a double eigenvalue X 1 = X 2 . The example shows 
how to find the valué of t, called x , where this happens. 

The starting point is a MATLAB test example, known as galle ry (3) . 



A = gall 


ery(3) 




A = 






-149 


-50 


-154 


537 


180 


546 


-27 


-9 


-25 



This is an example of a matrix whose eigenvalues are sensitive to the effects of 
roundoff errors introduced during their computation. The actual computed 
eigenvalues may vary from one machine to another, but on a typical 
workstation, the statements 
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format long 
e = eig(A) 

produce 

e = 

1 .00000000001122 
1 .99999999999162 
2.99999999999700 

Of course, the example was created so that its eigenvalues are actually 1, 2, 
and 3. Note that three or four digits have been lost to roundoff. This can be 
easily verified with the toolbox. The statements 

B = sym(A) ; 
e = eig(B) ' 
p = poly(B) 
f = factor(p) 

produce 

e = 

[1, 2, 3] 

P = 
x~3-6*x~2+11*x-6 

f = 

(x-1 )*(x-2)*(x-3) 

Are the eigenvalues sensitive to the perturbations caused by roundoff error 
because they are "cióse together"? Ordinarily, the valúes 1, 2, and 3 would be 
regarded as "well separated." But, in this case, the separation should be viewed 
on the scale of the original matrix. If A were replaced by A/1 000, the 
eigenvalues, which would be .001, .002, .003, would "seem" to be closer 
together. 

But eigenvalue sensitivity is more subtle than just "closeness." With a carefully 
chosen perturbation of the matrix, it is possible to make two of its eigenvalues 
coalesce into an actual double root that is extremely sensitive to roundoff and 
other errors. 
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One good perturbation direction can be obtained from the outer product of the 
left and right eigenvectors associated with the most sensitive eigenvalue. The 
following statement creates 



E = [130, -390, 0;43, -129, 0;133, -399,0] 
the perturbation matrix 



E = 






130 


-390 





43 


-129 





133 


-399 






The perturbation can now be expressed in terms of a single, scalar parameter 
t. The statements 

syms x t 
A = A+t*E 

replace A with the symbolic representation of its perturbation: 

A = 

[-149+130*1:, -50-390*t, -154] 

[ 537+43*t, 180-129*+., 546] 

[ -27+133*t, -9-399*+., -25] 

Computing the characteristic polynomial of this new A 
p = poly(A) 

gives 

P = 

x~3-6*x~2+11*x-t*x~2+492512*t*x -6- 1221271*+. 

Prettyprinting 

pretty (collect (p,x) ) 

shows more clearly that p is a cubic in x whose coefficients vary linearly with t. 

3 2 

x + (- t - 6) x + (492512 t + 11) x - 6 - 1221271 t 
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It turns out that when t is varied over a very small interval, from to 1.0e-6, 
the desired double root appears. This can best be seen graphically. The first 
figure shows plots of p, considered as a function of x, for three different valúes 
of t : t = 0, t = 0.5e-6, and t = 1.0e-6. For each valué, the eigenvalues are 
computed numerically and also plotted: 

x = .8: .01 :3.2; 
for k = 0:2 

c = sym2poly (subs(p,t,k*0 .5e-6) ) ; 

y = polyval(c,x) ; 

lambda = eig(double(subs(A,t , k*0 .5e-6) ) ) ; 

subplot(3,1 ,3-k) 

plot(x,y, ' - ' ,x,0*x, ' : ' , lambda, 0* lambda, ' o ' ) 

axis([ .8 3.2 -.5 .5]) 

text(2.25, .35, [ 't = ' num2str( k*0.5e-6 )]); 
end 
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The bottom subplot shows the unperturbed polynomial, with its three roots at 
1, 2, and 3. The middle subplot shows the first two roots approaching each 
other. In the top subplot, these two roots have become complex and only one 
real root remains. 

The next statements compute and display the actual eigenvalues 

e = eig(A) ; 
pretty (e) 

showing that e ( 2 ) and e ( 3 ) form a complex conjúgate pair: 

2 
1/3 492508/3 t - 1/3 - 1/9 t 

;i/3 %1 - 3 + 1/3 t + 2 

1/3 
%1 



1/3 492508/3 t - 1/3 - 1/9 t 

1/6 %1 + 3/2 + 1/3 t + 2 

1/3 
%1 



/ 2\] 

1/2 | 1/3 492508/3 t - 1/3 - 1/9 t |] 

+ 1/2 13 |1/3 %1 +3 - |] 

I 1/3 |] 

\ %1 /] 



1/3 492508/3 t - 1/3 - 1/9 t 

1/6 %1 + 3/2 + 1/3 t + 2 

1/3 
%1 



/ 2\] 

1/2 | 1/3 492508/3 t - 1/3 - 1/9 t |] 

1/2 I 3 |1/3 %1 +3 - |] 

I 1/3 |] 

\ %1 /] 
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2 3 3 

%1 := -2216286 t + 3189393 t + t +3 (358392752910068940 t 

2 4 1/2 

- 1052829647418 t - 181922388795 t + 4432572 t - 3) 

Next, the symbolic representations of the three eigenvalues are evaluated at 
many valúes of t 

tvals = (2:- .02:0) ' * 1 .e-6; 
r = size(tvals,1 ) ; 
c = size(e,1 ) ; 
lambda = zeros(r,c) ; 
f or k = 1 :c 

lambda(:,k) = double(subs(e(k) ,t , tvals) ) ; 
end 

plot ( lambda, tvals) 
xlabel( ' \lambda' ) ; ylabel('t'); 
title( ' Eigenvalue Transition') 



to produce a plot of their trajectories. 
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Eigenvalue Transition 




Above t = . 8e , the graphs of two of the eigenvalues intersect, while below 
t = 0.8e" 6 , two real roots become a complex conjúgate pair. What is the precise 
valué of t that marks this transition? Let x denote this valué of t. 

One way to find x is based on the fact that, at a double root, both the function 
and its derivative must vanish. This results in two polynomial equations to be 
solved for two unknowns. The statement 

sol = solve(p,diff (p, ' x ' ) ) 

solves the pair of algebraic equations p = and dp/dx = and produces 
sol = 





t: 


[4x1 


sym] 






x: 


[4x1 


sym] 




ind x 


now by 






f ormat 


short 






tau 


= 


double 


(sol 


t(2) 
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which reveáis that the second element of sol . t is the desired valué of x : 

tau = 

7.8379e-07 

Therefore, the second element of sol . x 
sigma = double(sol.x(2) ) 

is the double eigenvalue 

sigma = 
1 .5476 

To verify that this valué of x does indeed produce a double eigenvalue at 

o = 1.5476 , substitute x for t in the perturbed matrix A(t) = A + tE and find 

the eigenvalues of A(t) . That is, 

e = eig(double(subs(A,t,tau) ) ) 

e = 

1 .5476 + O.OOOOi 
1 .5476 - O.OOOOi 
2.9048 

confirms that o = 1.5476 is a double eigenvalue of A(t) for t = 7.8379e-07. 
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Solvíng Algebraíc Equations 

If S is a symbolic expression, 
solve(S) 

attempts to find valúes of the symbolic variable in S (as determined by 
f indsym) for which S is zero. For example, 

syms a b c x 

S = a*x"2 + b*x + c; 

solve(S) 

uses the familiar quadratic formula to produce 

ans = 

[1/2/a*(-b+(b"2-4*a*c)"(1/2))] 

[1/2/a*(-b-(b~2-4*a*c)~(1/2))] 

This is a symbolic vector whose elements are the two solutions. 

If you want to solve for a specific variable, you must specify that variable as an 
additional argument. For example, if you want to solve S for b, use the 
command 

b = solve(S,b) 

which returns 

b = 

- (a*x"2+c) /x 

Note that these examples assume equations of the form f{x) = . If you need 
to solve equations of the form f{x) = q{x) , you must use quoted strings. In 
particular, the command 

s = solve( ' cos(2*x)+sin(x)=1 ' ) 
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returns a vector with four solutions. 

s = 

[ 0] 
[ Pi] 
[ 1/6*pi] 
[ 5/6*pi] 

Several Algebraíc Equations 

This section explains how to solve systems of equations using the Symbolic 
Math Toolbox. As an example, suppose you have the system 

2 2 
x y = 

x-2- = a 
2 

and you want to solve for x and y. First créate the necessary symbolic objects. 
syms x y alpha 

There are several ways to address the output of solve. One is to use a 
two-output cali 

[x,y] = solve(x~2*y~2, x-y/2-alpha) 

which returns 

x = 

[ 0] 
[ 0] 
[ alpha] 
[ alpha] 

y = 

[ -2*alpha] 
[ -2*alpha] 
[ 0] 
[ 0] 

Consequently, the solution vector 
v = [x, y] 
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appears to have redundant components. This is due to the first equation 
2 2 = > which has two solutions in x and y: x = ±0 , y = ±0 . Changing the 



x y 
equations to 

eqsl = 'x"2*y~2=1, x-y/2-alpha 
[x,y] = solve(eqs1 ) 

produces four distinct solutions: 

x = 

[ 1/2*alpha+1/2*(alpha A 2+2P(1/2) ] 
[ 1/2*alpha-1/2*(alpha A 2+2P(1/2) ] 
[ 1/2*alpha+1/2*(alpha A 2-2P(1/2) ] 
[ 1/2*alpha-1/2*(alpha~2-2P(1/2) ] 



y = 

[ -alpha+(alpha*2+2)*(1/2) ] 
[ -alpha-(alpha"2+2)"(1/2)] 
[ -alpha+(alpha*2-2)*(1/2)] 
[ -alpha-(alpha"2-2)"(1/2)] 

Since you did not specify the dependent variables, solve uses f indsym to 
determine the variables. 

This way of assigning output from solve is quite successful for "small" systems. 
Plainly, if you had, say, a 10-by-lO system of equations, typing 

[x1 ,x2,x3,x4,x5,x6,x7,x8,x9,x10] = solve(...) 

is both awkward and time consuming. To circumvent this difficulty, solve can 
return a structure whose fields are the solutions. In particular, consider the 
system u~2-v"2 = a"2, u + v = 1,a"2-2*a = 3. The command 

S = solve( 'iT2-v"2 = a A 2','u + v = 1','a*2-2*a = 3') 

returns 

S = 

a: [2x1 sym] 

u: [2x1 sym] 

v: [2x1 sym] 
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The solutions for a reside in the "a-field" of S. That is, 
S.a 

produces 

ans = 
[ 3] 
[ -1] 

Similar comments apply to the solutions for u and v. The structure S can now 
be manipulated by field and index to access a particular portion of the solution. 
For example, if you want to examine the second solution, you can use the 
following statement 

s2 = [S.a(2) , S.u(2) , S.v(2) ] 

to extract the second component of each field. 

s2 = 

[ -1, 1, 0] 

The following statement 
M = [S.a, S.u, S.v] 

creates the solution matrix M 

M = 

[ 3, 5, -4] 

[ -1, 1, 0] 

whose rows comprise the distinct solutions of the system. 

Linear systems of simultaneous equations can also be solved using matrix 
división. For example, 

clear u v x y 

syms u v x y 

S = solve(x+2*y-u, 4*x+5*y-v) ; 

sol = [S.x;S.y] 
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and 



A = 


[1 2; 


4 


5]; 


b = 


[u; v; 


l; 




z = 


A\b 






result i 


in 






sol 


= 







[ -5/3*u+2/3*v] 
[ 4/3*u-1/3*v] 

z = 

[ -5/3*u+2/3*v] 

[ 4/3*u-1/3*v] 

Thus s and z produce the same solution, although the results are assigned to 
different variables. 

Single Dífferentíal Equatíon 

The function dsolve computes symbolic solutions to ordinary differential 
equations. The equations are specified by symbolic expressions containing the 
letter D to denote differentiation. The symbols D2, D3, ... DN, correspond to the 
second, third, ..., Nth derivative, respectively. Thus, D2y is the Symbolic Math 
Toolbox equivalent of d y/dt . The dependent variables are those preceded by 
D and the default independent variable is t. Note that ñames of symbolic 
variables should not contain D. The independent variable can be changed from 
t to some other symbolic variable by including that variable as the last input 
argument. 

Initial conditions can be specified by additional equations. If initial conditions 
are not specified, the solutions contain constants of integration, C1 , C2, etc. 

The output from dsolve parallels the output from solve. That is, you can cali 
dsolve with the number of output variables equal to the number of dependent 
variables or place the output in a structure whose fields contain the solutions 
of the differential equations. 
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Example 1 

The following cali to dsolve 

dsolve( 'Dy=1+y"2' ) 

uses y as the dependent variable and t as the default independent variable. 

The output of this command is 

ans = 
tan(t+C1 ) 

To specify an initial condition, use 
y = dsolve( ' Dy=1+y"2' , ' y (0 ) =1 ' ) 

This produces 

y = 
tan(t+1/4*pi) 

Notice that y is in the MATLAB workspace, but the independent variable t is 
not. Thus, the command diff (y,t) returns an error. To place t in the 
workspace, type syms t. 

Example 2 

Nonlinear equations may have múltiple solutions, even when initial conditions 
are given: 

x = dsolve( ' (Dx)*2+x*2=1 ' , 'x(0)=0' ) 

results in 

x = 

[ sin(t)] 

[ -sin(t)] 

Example 3 

Here is a second order differential equation with two initial conditions. The 
commands 

y = dsolve( ' D2y=cos(2*x) -y ' , ' y (0 ) =1 ' , ' Dy(0)=0 ' , 'x'); 
simplify (y) 
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produce 

ans = 

4/3*cos(x) -2/3*cos(x) "2+1/3 

The key issues in this example are the order of the equation and the initial 
conditions. To solve the ordinary differential equation 

3 
d u 

dx 

w(0) = 1, ií'(O) = -1, w"(0) = 7i 
simply type 

u = dsolve( 'D3u=u' , 'u(0)=1 ' , ' Du(0)=-1 ' , 'D2u(0) = pi ' , ' x ' ) 

3 3 

Use D3u to represent d u/dx and D2u(0) for u"(0) . 

Severa I Differential Equations 

The function dsolve can also handle several ordinary differential equations in 
several variables, with or without initial conditions. For example, here is a pair 
of linear, first-order equations. 

S = dsolve('Df = 3*f+4*g', ' Dg = -4*f+3*g') 

The computed solutions are returned in the structure S. You can determine the 
valúes of f and g by typing 

f = S.f 
f = 

exp(3*t)*(C1*sin(4*t)+C2*cos(4*t)) 

9 = S.g 

g = 

exp(3*t)*(C1*cos(4*t) -C2*sin(4*t)) 

If you prefer to recover f and g directly as well as include initial conditions, 
type 

[f,g] = dsolve( 'Df=3*f+4*g, Dg =-4*f+3*g', 'f(O) = 0, g(0) = 1') 
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exp(3*t)*sin(4*t) 



exp(3*t)*cos(4*t) 



This table details some examples and Symbolic Math Toolbox syntax. Note that 
the final entry in the table is the Airy differential equation whose solution is 
referred to as the Airy function. 



Differential Equation 



MATLAB Command 



^ + 4y(t)=e- t 
dt 



3/(0) = 1 



d y A , , -2x 

dx 



y(0) = 0,y(n)= 



d y , % 

— ^ = xy(x) 

dx 



y(0) = 0,y(3)= -K^S) 



(The Airy equation) 



y = dsolve( ' Dy+4*y = exp(-t) 
'y(O) =1') 



y = dsolve( ' D2y+4*y = exp(-2*x)', 
'y(0)=0', 'y(pi) = O 1 , 'x') 



y = dsolve('D2y = x*y','y(0) = O 1 , 
' y (3) = besselk(1/3, 2*sqrt (3) ) /pi ' , 
'x' ) 



The Airy function plays an important role in the mathematical modeling of the 
dispersión of water waves. 
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